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Introduction 



Strong interactions continue to be a challenge in theoretical physics. They involve all the phe- 
nomena regarding hadrons and actually the best suggestion for a model to describe them is given 
by Quantum Chromodynamics (QCD), which is a gauge quantum field theory based on a local 
SU(3) symmetry which satisfies the general requirements of renormalizability and unitarity of 
the scattering 5*-matrix. The elementary degrees of freedom of QCD are quarks and gluons and 
its complexity has prevented the solution of the theory i.e. finding the spectrum and correlation 
functions. On one side interesting work is devoted to analyzing numerically the euclidean version 
of the theory on a lattice, even if until now the available computational power really restrict the 
dimension of the lattice involved and thus the significance of obtainable results. 

On the other hand analytical efforts are done in two directions. Due to the asymptotic 
freedom property related to the decrease of the interaction strength with the increase of the 
energy of the particles involved, perturbation theory in the small coupling regime is applied. 
A perturbative approach has been also considered for other "small" parameters, such as 1/N C 
Jl|, ||, [| where N c is the number of colours which defines the special unitary gauge group. Typical 
non-perturbative methods concern condensates of quarks and gluons and topological objects such 
as instantons; however we are still far from the understanding the nature of composite bound 
objects in QCD. 

One of these objects seems to be the pomeron, originating from the works in Regge theory 
before a consistent field theory approach to strong interactions was established with QCD. The 
pomeron is related, in a hadronic scattering process, to the rightmost singularity (preferably a 
pole) of the partial wave amplitude in the complex angular momentum plane. The mathematical 
structure of the amplitude is commonly associated to the presence of an exchanged object with 
vacuum quantum numbers which can lead to a power-law-increase of the total cross section with 
the center of mass energy, in such a case in contrast with the unitarity bounds. In general the 
process is deeply related to non-perturbative phenomena and the only possibility to understand 
as much as possible of it with perturbative tools is to consider processes involving the scattering 
of highly virtual photons or onia-particles where the perturbative analysis in terms of elementary 
quarks and gluons make sense. 

At each order of perturbation theory one finds a leading behaviour proportional to a proper 
power of the logarithm of the center of mass energy (y/s). The approach which resums these 
dominant contributions is called leading log(s) approximation (LLA). One of its principal lim- 
itations is related to the fact that the coupling constant is considered fixed. The works of 
Lipatov et al. ||, |B| were devoted to these calculations for a SU (N c ) gauge theory. According to 
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previous ideas of Regge theory they use the gluon reggeization property, computed considering 
the dominant virtual corrections and which can be associated to a bootstrap condition related 
to unitarity, to construct the BFKL equation which describes the interactions in LLA between 
two reggeized gluons; the contributions to these interactions may be represented by ladder- type 
diagrams, resummed in a Bethe-Salpeter-like equation. The vacuum quantum number in the 
i-channel corresponds to the BFKL hard pomeron. Its singularity is not a pole but a branch-cut 
in the complex angular momentum plane. 

In order to go beyond the fixed coupling case a next-to-leading order (NLO) analysis has 
been carried on in the last years ||, [7|, ^. A quite complicated generalized BFKL equation 
was found but still it is not clear if such a model is really applicable to the HERA kinematical 
region or one may need to go beyond the NLO approximation to give satisfactory predictions 
and interpretations of the experimental data. 

All the theoretical predictions may be tested in various experiments in the kinematical Regge 
region. In the last years particular interest was devoted to the deep inelastic scattering (DIS) 
processes at HERA where the small x region is probed. Recent obtained measures [10| have 
shown a sharp rise of the F2 proton structure function at low x and this fact may be interpreted 
as a manifestation of the hard pomeron. Much attention was also given to the observed large 
rate of rapidity-gap events at small x, another possible signature of the pomeron interaction. 

The basic ingredients in the LLA analysis and a brief derivation of the BFKL equation are 
reported in the first chapter of this thesis. A well known problem related to the consistency 
of the approach is evidenced by noting that the diffusion-like process described by the BFKL 
equation may spoil the applicability of perturbation theory, since a huge contribution at large 
energies may involve small transverse momenta, deeply related to a non-perturbative regime in 
QCD. The other main problem in the BFKL theory is the violation of unitarity due to the high 
value of its intercept. The remarkable property of global conformal symmetry of the BFKL 
equation in the transverse two dimensional space is also shown and this permits one to also have 
a general explicit solution for the system. 

An alternative approach in perturbative QCD to describe the scattering of two onia-particles 



is given by the colour dipole picture Jllj, 12], which is based on the computation of the corrections 



to the light-cone wave functions of onia or virtual photons due to the emission of soft gluons 
strongly ordered in their scaling variables. This implies a leading logarithmic approach which 
turns out to be equivalent to the BFKL one. In the large N c limit it is possible to write 
an equation for the generating functional of the n-dipole distributions [11] which should be 
considered in general to have unitarity corrections. A brief introduction to the basic equations 
of this model is given in the last section of the first chapter, since they will be needed in the 
analysis performed in the third chapter. 

The problems related to the infrared sector one finds in the BFKL theory, which is con- 
structed on a small fixed coupling analysis thus far away from the non-perturbative region 
which instead can be reached in the diffusion precess, has motivated many phenomenological 
works, beyond the NLO analysis, devoted to the modification of the infrared sector of the BFKL 
equation. 

In the second chapter we try to analyze some properties of a phenomenological model which 
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have the merit of preserving the gluon reggeization property [|T3| , |l4| with the insertion of the 
running coupling in the generalized BFKL equation. The bootstrap condition is the crucial point 
since in this way unitarity is satisfied to leading order for the production amplitude with one 
gluon exchange in the i-channel while in all the other phenomenological attempts it is lost. The 
second input for the model is the compatibility with the predictions of the renormalization group 
equations for the gluon density (DGLAP) in the double leading log approximation (DLLA) The 
model depends at least on a mass parameter, related to infrared sector, on which the trajectories 
of all the singularities depend. 

The insertion of the running coupling provides a deep change in the properties of the pomeron 
equation [jl5| . A part of the spectrum is found to be discrete. The values are found numerically 
while the corresponding asymptotic behaviour of the eigenfunction can be derived analytically. 
The parameter dependence can be fixed, for example, by forcing the value of the slope (computed 
perturbatively) of the leading singularity. For some reasonable values one finds two supercritical 
pomeron states (corresponding to negative "energy" for the two reggeized gluons bound states in 
the Schrodinger-like picture). A rough attempt to estimate some effects of unitarity corrections 
has been made by considering the multipomeron exchange in the scattering of two virtual photons 
in the large N c limit, a case for which a corresponding eikonal formulation can be given. As a 
result, for the parameter choice considered, one finds that these corrections occur at energies 
much higher than the ones now-a-days reached. However unitarity requires one to include also 
the effects of the triple pomeron interaction, as is also shown in the third chapter. Thus unitarity 
effects can still be significant at lower energies. 

In the framework of the standard BFKL approach the process of jet production has been 
extensively studied showing the importance of mini-jet production and the logarithmic rise of 
the multiplicity. The fixed coupling analysis requires the introduction of an infrared parameter 
and the inclusions of some unitarity and coherence effects to restore the correct asymptotic 
behaviour of the differential cross section for large jet transverse momentum. In the model with 
running coupling considered one does not have, of course, any problem in the infrared sector 
and finds, due to the asymptotic behaviour of the pomeron states, a correct tail for the cross 
section at large jet transverse momenta. The calculations have been carried out on taking into 
account the two supercritical pomeron states. 

In order to obtain a better description of the physical processes at present energies a first 
improvement consists in considering all the states present in the spectrum of the pomeron 
equation, and not only the leading contributions of the two supercritical states. This program 
flTH can be achieved by considering the associated evolution equation in 1/x (with x the Bjorken 
variable) for the pomeron wave function. To this end we are forced to address the problem of 
the coupling of a virtual photon to the gluon in the running coupling framework, since otherwise 
there is no link, for example, between the gluon distribution in the proton and its structure 
function. The bootstrap relation, which concerns only the gluon interaction, is of no help. The 
idea is to find a corrected dipole density distribution for the photon in order to take into account 
a minimal running coupling effect. We adopt here a phenomenological approach and guess a 
possible form of the dipole density from the matching in the DLLA of the quark distribution 
(which one has after coupling the gluons to the dipole) with the one obtained by the DGLAP 
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evolution. Numerical calculations show a still steep increase of the gluon distribution and the 
singlet proton structure function, but smaller than the one implied by the standard BFKL 
theory. 

In the second part of the thesis we deal with the unitarity problem from a theoretical point 
of view. There is a limitation to the validity of the BFKL (fixed coupling) theory for large 
s or equivalently large 1/x and to extend the kinematical region of applicability of the model 
one has to go beyond the LLA. This problem is shared also by the DGLAP formalism and to 
solve it one needs to find a systematic way to include unitarity corrections. Continuing the 
restriction to the onia-particles or virtual photons scattering there is the hope to reach some 
new insights by means of perturbative tools even if the complete understanding surely requires 
a true non-perturbative approach (the idea of grasping non-perturbative effects with the help of 
a perturbative approach can be seen for example in the papers about renormalons). 

In order to understand which kind of corrections one has to deal with, it is very useful to 
consider the evolution of the parton densities in the rapidity. For example it is well know for DIS 
processes that at lower x the photon probes in the proton a higher density of partons, of typical 
size 1/Q at virtuality Q 2 . The resulting cross section then increases; the unitarity corrections 
correspond in this picture to the fusion of partons which diminishes the rate of increase of the 
parton density. A first quantitative model for such a process was given by Gribov, Levin and 
Ryskin |l8|] by constructing a modified non linear BFKL equation generating a fan diagrams 
structure (based on triple ladder vertices). The key point is the presence of multi-parton states 
in the t-channel and eventually also multi-scattering processes between the partons of the two 
scattering onia have to be considered as suggested by the dipole picture. 

To this end a systematic approach to take into account a minimal set of corrections necessary 
to restore unitarity has been studied by Bartels fl9[| . Higher order amplitudes, which contain 
any number of (reggeized) gluons in the t-channel, are constructed using unitarity and dispersion 
relations starting from the lower order ones and one has a coupled system of equations for the n 
reggeized gluon state amplitudes. The main feature is the presence of transitions vertices which 
change the number of reggeons in the t-channel and thus leads to a reggeon field theory. 

These corrections have been studied by considering a system of up to four reggeized gluons 
in the t-channel p0|. Among the building blocks an effective vertex for the transition between 



2 and 4 gluons states in the t-channel appears. In [21| its conformal invariance property was 
proven and in the appropriate color subspace it is related to the triple pomeron vertex. Another 
element of some interest is the 4-gluon interacting state (and its generalizations to n gluons) 
which, as a system with a fixed number of reggeized gluons in the t-channel, can be studied by 
means of the so called BKP equation [p2| , 23]. The related amplitude has been studied in the 
large N c limit and the system has been proven to be complete integrable and equivalent to an 



XXX Heisenberg zero spin model p4j , 25, gq]. These elements have suggested an interpretations 
in terms of correlation functions of a 2D conformal field theory (CFT) |27], ^] and one hopes to 
have an effective formulation in terms of a (2+1) dimensional CFT for QCD in the Regge limit. 

In the third chapter some efforts are addressed to obtaining a clearer picture of the interac- 
tions present in the system of 4 reggeized gluons in the t-channel in the large N c limit [^, |3(| , 
in particular on having some understanding of the relation between multipomeron exchange and 
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triple pomeron interactions and making a connections between this reggeon analysis based on 
the s-channel unitarity and the colour dipole picture for QCD. 

We shall study the formal general system of equations for the 4 gluon system at leading and 
next-to-leading order in 1/N C The double pomeron amplitude at NLO is studied to understand 
the role of the direct double pomeron exchange (DPE) and the triple pomeron interaction (TPI). 
The TPI term is analyzed in terms of 3 point function building blocks (a generalization to the 
non- forward case of a function previously introduced in |2C| ]) to describe in simpler terms the 
conformal invariance properties of the system. It is later shown in detail the fact, already noted 
for N c = 3 in f20|] , that the DPE contribution can be substituted by a completely equivalent TPI 
term while the inverse is not true and there is a complete agreement with the vertex form in |2(| 
in the limit N c — > do. On coupling to pomerons and working in the coordinate representation one 
can compare the results obtained from s-channel unitarity with the ones found by R. Peschanski 
in [31] starting from the colour dipole picture. A substantial agreement is found, apart from an 
asymmetry factor in the triple pomeron vertex and some single pomeron contributions absent in 
the dipole approach. A final section is devoted to studying some properties of higher order dipole 
densities in the A.H. Mueller approach which are found to be represented by a set of pomeron fan 
diagrams with only a triple pomeron coupling. It is expected that the s-channel unitarity analy- 
sis for a larger number of reggeized gluons in the t-channel might tell us the form of higher order 
pomeron vertices and thus also if the dipole picture has to be improved, as already suggested [31]. 



The thesis is structured in the following way. 

In the first chapter we review some important results of QCD in LLA which are at the basis 
of the developments in the second and third chapters. In section 1.1 the BFKL equation is 
briefly derived and in the following one the pomeron solution with some properties presented. 
Section 1.3 is devoted to show in some details the conformal invariance property of the BFKL 
equation and finally in section 1.4 a minimal review of the colour dipole approach in the A.H. 
Mueller formulation is given. In the second chapter we analyze a phenomenological model for 
considering the running coupling in the hard pomeron. In section 2.1 the definition of the 
model is given and the properties of the pomeron studied fll5f |; in particular in section 2.1.1 
and 2.1.2 the problem is formulated to make possible a numerical analysis, performed in section 
2.1.3, while multi-scattering effects are considered in section 2.1.4 In section 2.2 the inclusive jet 
production process is considered fi"q| . Section 2.2.1 is devoted to set up the formalism and to 
make asymptotic estimates while section 2.2.2 collects some numerical results. The evolution of 
the gluon density |l7]] is considered in section 2.3 and in particular in section 2.3.1 we address 
the problem of the coupling to the virtual photon, in section 2.3.2 we look at the starting 
distribution and in section 2.3.3 some numerical calculations are shown. Finally in the third 
chapter the triple pomeron vertex is analyzed in the large N c limit |30[| . In section 3.1 we 
introduce some basic elements of the unitarization program and in section 3.2 we present the 
analysis for the 4 gluons system in the large iV c limit, in particular at the leading order in 3.2.1 
and at the NLO in 3.2.2 where we establish the form of the DPE and TPI terms. In section 
3.3 the properties of conformal invariance are re-analyzed in terms of a simple function and in 
section 3.4 the interplay between the DPE and the TPI terms is analyzed in detail. In section 3.5 
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the coupling to pomerons and the relation to the dipole approach is studied while section 3.6 is 
devoted to study the solutions for general higher order dipole distributions in the standard dipole 
picture and analyze the triple pomeron interactions structure resulted. In the appendix one can 
find some more technical parts. In A.l, related to the phenomenological model considered in 
the third chapter, the asymptotic form of the bound state solution of the pomeron equation is 
derived. In A. 2 the dipole density in the non-forward case for the coupling of two gluons to a 
quark loop is given. In A. 3 the coordinate representation of the function G, needed in chapter 
4, is calculated and its infrared stability explicitly shown. In A. 4 some relations, used in the 
fourth chapter, are derived in a graphical notation. 



Chapter 1 



Pomeron and the Regge limit in 
perturbative QCD 

To understand the concept of Pomeron it is necessary to go back to the old Regge theory 
and Pomeron field theory. The Pomeron is the rightmost singularity in the complex angular 
momentum plane of the partial wave signatured amplitude of a hadronic scattering process 
characterized by the quantum number of the vacuum exchanged in the i-channel. The singularity 
can be associated to a rise of the cross-section with the center of mass energy, thus leading to 
a violation of unitarity which can be restored only if more complicated processes than a single 
pomeron exchange are taken into account. The understanding of the strong interactions in the 
Regge limit is therefore deeply connected to unitarity from a general point of view. 

The study of strong interactions as described by QCD theory is well known to be an incredibly 
difficult problem and due to the asymptotic freedom which characterizes QCD many efforts have 
been directed in perturbative analysis, hoping to not lose some important aspects, which could 
serve as a guideline to understand the features of the complete real QCD. 

In particular the concept of the Pomeron has been intensively studied in perturbative QCD 
in the Regge limit for a long time H. At the basis of the approach is the property of Reggeization 
of the gluon. In general the poles in the complex angular momentum plane of the partial wave 
amplitude depend on the exchanged t squared momenta so that they may describe a family of 
exchanged particles of integer even or odd angular momenta, thus rising with spin and mass. An 
elementary particle, when exchanged in the t channel, gives a non-analytic contribution to the 
partial amplitude dependence in the angular momentum, and the modern idea of reggeization 
is related to the fact that summing all the contributions to the amplitude not containing poles 
in the t-channel leads to a partial wave amplitude with a pole in the t-channel which is located 
at the spin of the elementary particles if t is taken to be the mass squared of the particle itself. 
This was verified for the gluon (massive after considering the Higgs mechanism) up to the 8-th 
order (and partially up to 10-th) in perturbation theory by Cheng and Lo |32|| , after the work 
of Grisaru |3^] where a first confirmation that the gluon reggeizes was given. 

Another method which has been used to check gluon reggeization consists in taking a 
gluon ladder diagram, constructed with reggeized gluons on the basis of elastic unitarity (the on 
shell condition turns the reggeized gluons into elementary gluons). On summing these ladder 
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diagrams one obtains a partial wave amplitude which, in the colour channel, has the same 
reggeized gluon Regge pole. This is called a "bootstrap" condition and it is at the basis of the 
BFKL theory and all its extensions. 

The analysis of the processes in the Regge limit in perturbation theory reduces to summing 
the leading log terms. In fact for a small fixed coupling g and high c.m.s. energy s such 
that s — ► oo, g — > and g 2 \ns ~ 1 the Regge pole contribution has a behaviour s a ^ where 
a(t) = l + u(t) is the trajectory and u[t) is of second order in the coupling constant g. Therefore 
a perturbative expansion leads one to consider amplitudes of the order g 2 s{g 2 lns) n . 

Apart from considering a fixed coupling, its smallness is generally ensured in the Lipatov 
perturbative approach by giving to the gluon a mass M >> A using the Higgs-Brout-Englert 
mechanism [35fl . Thus one has also an infrared regularization but we should remember that 
in this way one is not really working with true QCD since it is not possible to safely perform 
the m — > limit without avoiding the non-perturbative region. 

The LLA analysis of the amplitude for gluon-gluon scattering in the vacuum channel was 
started in (|, 0]. Multiparticle production and unitarity were the ingredients used to construct 
a Bethe-Salpeter equation for ladder contributions to the partial wave amplitude with reggeized 
gluons in the t-channel. The resulting singularity was found to be a branch cut and not a pole in 
the complex angular momentum and the BFKL pomeron to be a "bound" state of two reggeized 
gluons. In the colour octet channel instead the Regge pole ansatz for the gluon applied to the 
BFKL equation gave directly the self-consistency bootstrap condition. 

The BFKL equation can be seen as a Schrodinger-like equation in the transverse space 
with the evolution flow in the rapidity variable. We just remember, beyond the violation of 
unitarity, an intrinsic problem in the BFKL perturbative approach due to the diffusion of the 
gluon momenta in the ladders. For a highly virtual projectile (Q 2 ) the scattering is characterized 
by a distribution which tends to be Gaussian in In k 2 with an average value In Q 2 and a width 
proportional to vrn s, thus with a support which can be strongly present in the non-perturbative 
region. 

The opportunity of applying the results found in this perturbative approach to analyze some 
real processes is based on the possibility of separating in them the contribution of short distance 
interactions (hard scales), where perturbation theory can be applied, from the long distance 
ones (soft scales) where perturbation theory is meaningless. 

This is the meaning of the factorization theorems. A first kind used to study many hard 
processes characterized by a hard scale Q 2 at which the coupling is small is the collinear factor- 
ization theorem ||36|| . An object such as, for example, a structure function in DIS can be written 
as a convolution between a cross section for an hard parton subprocess (perturbatively calculable 
and depending on the process virtuality) and a distribution for the parton probed in the hadron 
involved. The latter distribution depends on the virtuality chosen to probe the hadron by means 
of a renormalization group equations (DGLAP [ |37|| ) so that a non-perturbative input at a lower 
scale Qq is necessary to provide the initial conditions for these equations (which are constructed 
on all the collinear logarithmic contributions). The DGLAP equations can be associated, in 
the physical gauge, with summing all the contributions of ladder graphs with strongly ordered 
parton momenta but, in the small x (high energy) limit, higher order contributions, beyond the 
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leading logs, are found to be significant. This fact has led to consider the perturbative resum- 
mation of contributions of the order (a s lnl/x) n which correspond to the ladder graphs with 
reggeized gluons in the i-channel and no ordering in the transverse momenta. In this framework 
a new k-factorization theorem has been given pg|]. 

We shall introduce in the next section the BFKL equation and show some of its properties 
and the solutions in the subsequent one. In the third section the proof of its conformal invariance 
is given, a feature which turns out to be very interesting and also characterizes the symmetry of 
the three pomeron vertex. We shall return to these properties in the last chapter. The Colour 



dipole approach in the A. H. Mueller formulation [Tl], 12] is briefly reviewed in the last section 
since its use is pertinent in the high energy (and a large number of colours) limit and it has 
been shown to give an alternative formulation to the hard pomeron phenomena. In chapter 3 
we shall try to compare the pomeron interactions described in the two formulations to see the 
degree of equivalence of the two approaches. 



1.1 The BFKL equation 

The main results about the LLA are briefly derived in the following || |3^]. We consider the 
absorptive part A s (s, t) = ImA(s, t) of the dominant non-helicity-flip amplitude for the scattering 
of two gluons in the limit s — ► oo, t finite and a s = <7 2 /(47r) << 1 in the t-channel colour state 
R. Using unitarity one can compute this quantity by means of the amplitude for the production 
of (n + 2) gluons from two: 

p(R) r 

ImA (R \s,t) = / dT n+2 A*(n + 2,p' 1 ,p' 2 )A(n + 2, Pl ,p 2 ) (1.1) 

where = is the colour projector, p\ = (p+,0, 0j_) and p 2 = (0,p_,0j_) are the momenta of 
the incoming gluons in the c.m.s. in the light cone representation, so that s = 2p\ ■ p 2 = 2p"^_ 
and d,T n+2 is the (n + 2) particle phase measure given by 

n+l 

dr n+2 = J] (d 3 k i /(27r) 3 2k i+ )(2TT) 4 5 4 (J2 h-pi- P2) (1.2) 

i=0 

The momenta k% of the intermediate particles satisfy an on-shell condition. It is convenient to 
use the Sudakov variables on and Pi such that k{ = aip\ + flip 2 + ki±, oiiPi = rnfj_/s with m?j_ = 
m? — kf ± . Other useful variables are the transferred momenta qi = p\ — S}=o % (^» = 1i~ 
and the partial energies carried by two adiacent intermediate real particles s^j+i = (ki + hi + i) 2 . 

In the Regge limit the dominant leading log(s) contribution comes from the region in the 
phase space characterized by 

OLi » CKi+1 , Pi « Pi+l (1.3) 

and thus one has approximatively the equalities Sj,j+i = (ai/ai + i)mj +1 ± and qf = qf ± . Chang- 
ing the variables which represent the measure in phase space to % and sa+i and performing an 
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integration in the tranverse space and the integrations w.r.t. «o an d Oi n +i one finally may write 

dr n+2 = 7 n ^y^ ^ii s ^+ i - s n m ?±) ( l4 ) 

Let us now consider the amplitude for the production of (n + 2) gluons. In particular we 
shall denote by Aj, aj and A^, a! i the helicities and colours of the two incoming gluons and of 
the two outgoing gluons with nearly the same momenta as the initial ones and by and di the 
polarizations and colours of the remaining n gluons. 

In the Regge limit the virtual corrections to the tree diagrams may be analyzed by looking 
at the one-loop corrections of the scattering process gg — ► gg with both the diagrams related to 
the s-channel and u-channel physical region. One finds a correction of the type a(t) Ins and it 
is natural to make the reggeization ansatz for a correction (s/so) a ^ in all orders in a s in the 
LLA. The trajectory found for the gluon is a(t) = 1 + to(t) with 

^= n ^I ^iZ-i? ; t= ~ q2 (L5) 

Alternatively one may use a general reggeization ansatz for the gluons in the i-channel and find 
from the boostrap condition the trajectory, for a case of the two gluons amplitude in the colour 
adjoint representation; this confirms the fact that the gluon reggeizes. 

The real corrections in the LLA leads one to consider a three gluons reggeon-physical-reggeon 
vertex (the reggeons in the t-channel) whose explicit expression is 

F(q i+ i,qi) = -q i± - ft+i^+Pi + r^) ~Pz( + ~, (1-6) 

PlP2 hPl VW2 kP2 

It is orthogonal (F(q^i,qi)ki = 0) to the emitted gluon momentum. Prom the analysis of the 
two physical-reggeon-physical gluons vertices (external in the ladder) the dominant contribution 
is helicity conserving and proportional to s so that the amplitude A2—>2+n can be written as 

A 2 ^ +n = 2sg5 Xl yT^ ai s^h^gT^ 

4j(t„ + l) -1 j. rpln + 1 / 1 7 \ 

- s n,n+l t n+l9 6 X 2 X' 2 1 a ' 2 a2 ^ ' > 

where ti = —qf In order to obtaining an expression for ( |1.1| ) it is necessary to sum over the 
polarizations and the colours of the intermediate real gluons and consider the projection onto 
the desired colour subspace (we shall be mainly interested in the vacuum state). 

The sum over polarizations leads to an effective interaction vertex between the two reggeized 
gluons. Starting from Ylxt ^^xj^^i ili-* Qi+l)^ (lij Qi+i) an d performing the sum over the po- 
larizations (actually using the replacement e Ai e ^ ~~ * —g^ iUi ) one obtains 

F-F' = -2K q (q i ,q / i ,q i+ i,q' i+1 ) ; K q (q h q[, q i+1 , q' i+l ) = q 2 - (qf ± q'l± + q%±q'i±)/ k ± ( L8 ) 

Regarding the colour structure in the t-channel, the gluon-gluon amplitude has to be de- 
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composed into irreducible representations of the (N% — 1) (8> (N% — 1); by means of a projection 
operator one has 

4* W,t) =T,Cia ia ^ R) (s,t) (1.9) 
(R) 

For each real gluon emission one has a colour factor which depends on the representation R 
according to the decomposition of the product of the two corresponding colour matrices 



di 

so that the interation term is given by 

= \ R K q (i.ii) 

Also for the first and last intermediate gluons one has the same colour factor. In particular one 
has for the vacuum channel Ai = Aq = N c and for the adjoint representation A^v2„ 1 = N c /2. 



(R) 



The absorptive amplitude (|l.l| ) can thus be written in the form 

J z^zir) sq tit i 



IJ 2K^( qi ,q l+1 )5{\{ s i)i+1 - s J] (1-12) 

1=1 i=0 1=1 

and has the signature fixed by the R representation (positive for symmetric and negative for 
antisymmetric representations). Taking the Mellin transform one obtains for it an asymptotic 
contribution (which thus corresponds to the Mellin transform nearby the rightmost singularity) 

. n+l 



n (R)(n) M _ o„2(n+2) \2 / TT " q 



(t)-2 5 2 (^)Ai/n iwu*iV-i- 

f[\k^r l 2K^\ qi ,q i+1 ) (1.13) 

i=i 

We note that in the validity region (j close to 1) the factor may be omitted. Changing 

the notation to the 2-dimensional transverse vectors In = qi±, In = —Qi± such that In + l2i = q, 
U = —I 2 and t = —q 2 one may write in an operatorial form 

a™ n \t) = 2g 2 X 2 R I ^^(hi\D(V^Dr\h, n+1 )^^ (1.14) 



167T 3 tiit21 ' 167T 3 



where we have defined the operators 



= : y^)^^i (1 . u) 
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It is now very easy to sum the contributions of all the ladder diagrams and one has 



,2, A*) 



with 

72,/ 



= J J^^ 1 - V (R) D)-%) (1.17) 

The Lipatov equation is a Bethe-Salpeter equation for 4>) q \h) which can be obtained by taking 
the matrix element of the operatorial identity D~ 1 T^ = l + V^T^ w.r.t. the states and 
where = D(l - V^D)' 1 and integrating w.r.t. l[. It reads 

\j-l- u (t 1 )- U (t 2 )]4f(h) = l + 2a a X R J ^V q {hA)^\l'i) (1-18) 
with the interaction potential given by 

v < M = 4 + Pwk?-$7 (1 ' 19) 

This equation in the gluon channel gives the bootstrap condition if one requires that the solution 
of this equation must have a Regge pole of a reggeized gluon, i.e. at j = 1 + uj(t). Denoting by 
Cq(li) the residue of the solution at this pole, the bootstrap condition can be written for N c = 3 
as f 

[u(t) - w(ti) - u{t 2 )]c q {h) = 3a s J ^V q {hA)cS'i) (1-20) 



The gluon trajectory found by perturbative calculations corresponds to a residue c q independent 

31 



of l\. Hence in the gluon channel the whole solution for 4>f„ A \h) is determined by this pole. 



1.2 Solution for the pomeron amplitude 

The Lipatov equation can be written as an inhomogeneous Schrodinger equation (H — E)(f> = f 
with an "energy" E = 1 — j and an Hamiltonian 

H = -u(h)-u(l 2 )-2a s \ R V q (1.21) 

In the vacuum case one has \r = N c and one may check that for q / 0, as will be shown also 
in the next section, the infrared infrared divergences cancels between the gluon trajectories and 
the interaction potential term. For q = the divergences at l[ = 1' 2 = remain but actually 
disappear for the coupling of the pomeron to a colourless physical object, since the function 
4>j q (l — > 0, q) vanishes. 

Let us consider the forward scattering case (q = 0). It is convenient to consider a new 
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function ip(li) = for which one can write the equation 

1 f d 2 l' ibil' ) 

[j - 1 - 2u{t x )]ij{h) = ^ + 4a s X R J ^ ^ (1.22) 

The interaction term is of the Coulomb type and both the trajectory and the potential can be 
studied in the limit of zero regularizing mass. Using the expressions given in the next section 
( Jl.35 ) and (1.37) one may write the Hamiltonian for the vacuum colour channel (Ao = N c = 3) 

H = (3a s /ir)h ; h = \nl 2 + hip 2 - 2(ln2 + ^(1)) ; E = 1 - j = (3a s /vr)e (1.23) 

where the dependence on the gluon mass has disappeared and is referred to the logarithmic 
derivative of the T function. To find the eigenvalues and eigenfunction of this Hamiltonian one 
can take advantage of the scaling symmetry p —* cp , Z — >• Z/c which is manifest. The symmetry 
between the coordinate and the momentum space is complete and one may write easily the 
proper functions for the system both in coordinate and momentum space 

^n(p) = a vnP - 1+iv d n t> ; i> vn {l) = K n r l - iv j M (1-24) 

for any real v and integer n. It is also easy to show that these functions are orthonormalized if 
\a vn \ 2 = \b un \ 2 = 1/47T 2 . Defining \i = (1/2)(1 + n + iv) the relation between the two normalizing 
factors is for n > 

b un = a vn i n T w ^ (1.25) 

while the n < case is obtained by complex conjugation. The eigenvalues are computed through 
the definition h Vw = £un ipun \ in particular one has 

e un = 2ReV( 1 + |n 2 l+U/ ) - 2^(1) (1.26) 

The minimal value eo = —4 In 2 is reached at v = and n = and gives for the BFKL pomeron 
the intercept 

j = l + ^41 n2 (1.27) 

From this value the violation of the Froissart bound |4j|] derived from unitarity is evident. 

The pomeron Green function in the forward direction can be easily constructed considering the 
inhomogeneous term of the equation for tp, f = (2tt) 2 5(1 — l')l~ 2 l'~ 2 , and in momentum space is 

To see which behaviour the solution has in the high momentum (short range) limit one 
should study the singularities of the Green function in the v plane. Since we are interested in 
the dominant contribution, which one has for e close to eo, we consider the expansion up to the 
second order 

e„o = eo + au 2 (1.29) 
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where a = (7/2)£(3). Defining w = eo — e = (n/3a s )(j — 1) — |eo| the singularity is at w = and 
the asymptotic behaviour for large I is calculated by taking the residue at v = iyjw/a. Thus 
one obtains 

= V e** f °° d./r^fcj^ * (^ e -Vf (1.30) 
^ J-oo £yn - e 21 Jaw 



n=—oo 



and it is evident that the pomeron singularity is branch cut of type — 1 — |eo|. The 

leading behavior at high energies s can be found by performing an inverse Mellin transform for 
the positive signatured amplitude 



i>(s,l) = - - WiW ■ (1-31) 
4 j sin 7rj 



giving the dominant contribution 

nr 1+ 3a, ko i3a s (tpoo\f) -tt '"^ , n QrA 

1P(S,1) ~ S ^ T 1 01 - P. 12a s alns (1-32) 

4/ 7r \/3a s alns 

Thus the average value of I increases with s since the distribution is Gaussian in In/ with a 
width proportional to Vm s. This fact is characteristic of a diffusion process in the transverse 
momenta; the Lipatov equation can in fact be seen as an diffusion equation in the rapidity flow. 
From this fact we address an important problem in the BFKL theory. Such a diffusion process 
may spoil the applicability of perturbation theory since at sufficient large energy the distribution 
in the transverse momenta can have support in the small values region corresponding to large 
distances and typically regarding a non-perturbative regime. 



1.3 Conformal Invar iance 



The BFKL equation manifests two important properties: the holomorphic separability and con- 



formal invariance. An elegant way to prove them [40, ^7], ^] can be achieved by rewriting the 
Hamiltonian of the two interacting reggeized gluons in the complex coordinate space represen- 
tation for the transverse plane. In the following we will fix the notation: given two vectors 
a = (a x , a y ) and b = (b x , by) we have 



a = a x + ia y ; ,a.-b = -(ab* + a*b) (1.33) 



The complex vectors for the impact parameter r{ will be conjugate to the momenta Zj = —i-gp- 
It is also convenient to introduce the relative momentum I = — h) and relative distance 
p = ri — T2 (vector) and p = \{r* — r|) (complex number) so that p,l and p*,l* forms, each one, 
an independent pair of canonical variables and \p\ 2 = App* . It should be clear from the context 
which is, each time, used. Further if in the configuration space we are referring to a vector or 
to a complex number should be clear from the context. Moreover we define k = l\ — l[. 
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The interaction potential ( |1,19| ) can be rewritten as 

V g (h, Zl) = hll^l'xl'l)- 1 + c.c. = h^WipWlr 1 + c.c. (1.34) 

where W(p) is a Coulombian potential in the coordinate representation. The explicit form of 
W, after having introduced an infrared regulator m, is 

W(p) = *p» „ J. (1.35) 

27T m—>0 27T L 2 J 

and consequently one obtains 

1 r 



/iln / 9/f 1 +/ 2 lnp^ 1 + lnm 2 -2V'(l) + c.c. (1.36) 



The regularized form of the gluon Regge trajectory reads 

"(h) « —^(lnh-lnm + c.c.) (1.37) 

so that the Hamiltonian is holomorphic separable, i.e. H = H(p,li) +c.c, and, for the R colour 
channel, 

^,^) = ^fln/ 1 +ln/ 2 + ^(Z 1 lnp^ 1 +^lnp/ 2 - 1 ) + (^-l)lnm 2 -2^V'(l)l (1-38) 

In the vacuum channel the infrared divergences (regularized by m) in the interaction and the 
trajectory terms cancel and on rescaling the wave function <fi in the form <fi = lil^ip, the new 
Hamiltonian for if) will be given by H' = ^^-(^12 + ^12) with 

hi2 = In Zi + In l 2 + li 1 In p l x + Z^ 1 In p Z 2 - 2^(1) (1-39) 

The holomorphic separability of the Hamiltonian is an important property and as a consequence 
the eigenvalues will be given by the sum of two contributions, one the complex conjugate of the 
other, and the corresponding eigenfunctions will be factorized in two parts, the first depending 
only on ri,r 2 and the second part on the conjugate quantities. 

Let us now consider a global conformal transformation in two dimension, represented by the 
Mobius mapping 

, ar + b j , 

r^r'= 1.40 

cr + d v ' 

where a,b,c,d are complex number such that ad — be = 1. The generators of the Mobius 
group are M° = rl, M~ = /, M + = r 2 l where I = Mo generates the rotation and scale 

transformations and M~ the translations. The Casimir operator is given by M 2 = M° 2 — 
^{M + ,M~}. For a system of n gluons the generators are given by the direct sum of the 
generators acting on the different gluon spaces. 

The Hamiltonian , which describes the two gluon system, is evidently invariant under the 
action of the Mq and M~ generators (also /ii 2 given in ( |1.39| ) is invariant). To prove the 
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conformal invariance it is convenient |4(], [39| to consider the finite transformation denned by 
r — > — - which leads to 

P^— ; k^rfh (1.41) 
r x r 2 

and rewrite the Hamiltonian in a different form. Let us therefore consider the following operator 
identities 

d l =T- 1 (l + z)p- 1 T{\ + z) -> In ft = -r- 1 (l + z)ln / or(l+z) 

p 2 di = T(z)pT~ 1 (z) -> lnp 2 di = r(z)lnpT- 1 (z) (1.42) 

where d\ = and z = pd\. They are similarity transformations between p~ l and d\ and 
between p and p 2 d\. By noting that [z, \np] = 1 one can use the representation lnp = — J| and 
derive from the relations ( |1.42| ) 

ln<9i = — In p + ^(z)-\ — 
lnp 2 <9i = In p + ip(z) (1.43) 

Introducing these relations in ( 1.39| ) and also substituting the operator identity l^ 1 In pl\ = 
In p — i the holomorphic part h\ 2 can be written as 

h 12 = lnp 2 9i+ lnp 2 9 2 -2 lnp -2^(1) (1.44) 

In this last expression we made use of the cancellation of imaginary terms between the holomor- 
phic and antiholomorphic part. 

It is now easy to prove the invariance of the Hamiltonian with respect to the transformation 



Ti — * --ir A 



1,2. In fact applying the transformations ( [L.41 ) one has 



^ ^ ; P 2 d 2 ^ P -° 2 
and the conformal invariance is verified by direct substitution 



Once the conformal invariance is known to be a property of the system, one can take advan- 
tage of it by expanding the general solution of the BFKL equation in terms of the conformal 
basis functions, which are eigenfunctions of the Casimir operator of the Mobius group and are 
then related to its irreducible representations. 

For a system of two gluon the Casimir operator for the holomorphic part is given by M 2 = 
—pl\l 2 . Using the notation p = r\ 2 = T\ — r 2 the conformal basis is defined by 

M 2 M 2 E^ n ,r 2 ) = ^rf 2 d 2 d 2 E^ ri ,r 2 ) = ^p(p - l)p(p - l)E^(n,r 2 ) (1.46) 

where p = + iv and p = + iu are the conformal weights which label the representation; 
n, the conformal spin, is an integer and v is a real number. The basis functions (in complex 
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notation) are explicitly given by 

E,(r u r 2 ) = E m (r 1 ,r 2 ) = (^f-) 2 2 (1.47) 

V ^10^20 ' Vr 10 r 20 7 

One notes the presence of an additional quantum number, the coordinate tq, in some way 
related to translational invariance. The corresponding eigenvalues of the Casimir operator can 
be written conveniently to obtain 

M 2 M 2 E^n,r 2 ) = - E^ n ,r 2 ) (1.48) 

^Q"n J r\,vQ"n— l,i> 



where we have used the standard notation 



7T 4 1 



1.49) 



2 v 2 + n 2 /4 
The conformal basis forms a complete system: 

rf 2 5 2 (r u ,)5 2 (r 22 ,)=Y,Mr 1 ,r 2 )E;(r[y 2 ) (1.50) 
The notation used regarding the sum over the quantum numbers is the following 

5> £ \ dv ^-\ d2r » (L51) 



l-l 



Another important property is the orthogonality relation which contains an additional term due 
to the equivalence of the representations which correspond to the (n, v) and (— n, —v) quantum 
numbers. It is given by 



f d 2 nd 2 r 2 

/ — i — En^r^n^^E^yyin,^) = 

a n , v 8 nn ,8(y - v')8 2 {r w ) + b n ,J n ^ n ,5(v + z/)|r 0(y r 2 ~ 4 ^(^) n (1-52) 

r 00' 

where 

= 3 2 ^ r(-^ + (l + |n|)/2) r^ + Jnj/2) 

-M/ + M/2 r(w/ + (l + |n|)/2) r(-ii/+|n|/2) l ' J 



The proof of these relations can be found in|40[]. The spectrum of the BFKL Hamiltonian for 
the general non-forward case (q 7^ 0) is the same as the forward case, as a consequence of the 
global conformal invariance: 

JV,a. r , , s ,l + |n| , , ,1 + lnl ,i 
v» = u n ,u = 2^(1) - ^i^- 1 + iv) - ^{^r 1 - iv) (1-54) 

The solution of the BFKL equation with an inhomogeneous term ^/q can then be expanded in 
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terms of the conformal basis functions: in the Mellin representation one has 

*j(ri,r 2 ) = 2^ . _ _ (1.55) 

and moving to the rapidity space representations, where the BFKL equation is an evolution 
equation in the rapidity variable, one finds 

*(n,r 2 ;y) =^e^^(r 1 ,r 2 )( / u|^ ) (1-56) 
It 

Here we have defined 



<A*l*o> = / -^^^(ri.raJ^oCri,^) (1.57) 



1.4 The colour dipole approach in the large N c limit 

The light-cone wave function of colourless objects, such as onia or virtual photons, has been 



studied by many authors [11, 12, Q to compute the corrections to the wave function given by the 
emission of soft gluons with scaling variables strongly ordered. The kinematical region >> 
Zk thus relates the results found to the one obtained in the BFKL approach. Incidentally the 
derivation is quite different and simpler and avoids some suspicious assumptions taken in the 
construction of the BFKL theory. 

The colour dipole approach has also the merit of giving, without further effort, some insight 
in the multipomeron exchange and interactions as we shall see in chapter 3. In this section we 
shall sketch a brief derivation of the model in the approach given by A.H. Mueller 0, g|, || 
and show its relation to the BFKL pomeron. 

The first object introduced is the wave function of the colourless object (onium or virtual 
photon) with momentum (p + ,0,0) in the lowest order in the coupling, with no soft gluons 
emitted, so that we may think of it as split in a quark-antiquark pair. Let k\ and —k± denote 
the transverse momenta of the antiquark and of the quark respectively; the scaling variable for 
the antiquark is defined by z\ = k\ + /p + . We shall mainly use the coordinate representation in 
which the wave function reads 

where a and (3 are indices for degrees of freedom such as spin and colour. The incoherent sum 
over spin and colour of the square moduli of the wave function is 

^°\r 1 ,z 1 )=J2\^(r 1 ,z 1 )\ 2 (1.59) 

a/3 

and the normalization will be given by 
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For a hadron N is finite and can be chosen to be 1 while for a photon it would diverge. 



Let us consider the contribution to order a s due to the emission of one soft gluon. One will 
search for a component of the wave function of the form V , i/j a {t\ , Z\ ; r 2 , z 2 ) where r 2 , Z2 and a are 
the transverse coordinate, the scaling and the colour variables of the soft gluon (z 2 « Z\, 1 — Z\) 
respectively. Summing the contributions due to the emission of a soft gluon from the quark and 
the antiquark one finds 

^ a (n,z l ;r 2 ,Z2) = -i{g/ir)t a ^l{r 1 ,z 1 )[{r 2 e)/rl - (r 21 e)/r 2 21 ] (1.61) 

where e is the polarization vector of the soft gluon and rj,- = r% — rj. Introducing the quark 
coordinate r$ = 0, taking the square modulus of the wave function and the sum over colour and 
spin one finds the distribution 

2 

$W( ri , Zi;r2)22 ) = $(o)( ri)Zl ) c J^_ (i. 6 2) 

r 20 r 21 

where c = 4a s C F /TT, C F = (N 2 - 1)/(2JV C ). 



In order to consider the contributions due to the multiple emissions of soft gluons some 
simplifications have been taken into account so that it is possible to find an integral equation 
which resums all the orders. Precisely the large -/V c limit has been considered, under which the 
gluon lines can be represented in terms of qq pairs and only diagrams with a planar topology 
survive. 



For a second soft gluon emitted (characterized by r^,z^) one finally obtains 

^2 ^2 



^ 2 \r 1 ,z 1 ;r 2 ,z 2 ;r 3 ,z 3 ) = $W(n, z\\ r 2 , z 2 )c(-^\- + -^V) 

r 03 r 32 r 13 r 32 

= * (0 Hn,z 1 )c 2 rl(- T ^ + - T -^) (1.63) 

'02'23'31 '03'32'21 

and the formula for an emission of n soft gluons in the kinematical region z n+ \ « z n « ...z 2 
is easily inferred 

$ (n) (ri,zi;r 2 ,z 2 ;...r n+ i,z n+ i) = ^°\ ri , Zx)c n r 2 J 2 2^ 2 V perm, of 2,3, ...n + 1 

r 12 r 23-" r n+l,0 

(1.64) 

In the strongly ordered region the integration in the scaling variables of the emitted soft glu- 
ons can evidently be done and in the leading order one obtains a trivial logarithmical factor 
(lnzio) n /n! where z\q = Z\/zq. It is possible to introduce [11] a generating functional for the 
distribution of the transverse coordinates of the emitted gluons: 

n+1 s 

& n Xr 1 ,z 1 -r 2 ,...r n+1 -z ) = &°Xr 1 ,z 1 )(l/nV) ]J ——Z(r 1 ,r ,z 10 ) u=0 (1.65) 

2 OU(Vi) 
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The generating functional satisfies the following nonlinear equation 

Z(ri,ro,zi ) = 1 + cr\ [ — [ d ''I' 2 „ u(r 2 )Z(r 2 , r , z 20 )Z(r 2 , n, z 20 ) (1.66) 



One may be interested in a functional generating not the distribution in terms of the transverse 
positions of the emitted gluons but instead in terms of r^, already present in ( 1.66j ), which 



correspond to the dimensions of the colour dipoles. The proper generating functional in this 
case is 

D{no,z 10 ;u) = u(r 10 ) + crf / — - / - — 2 -^ L 2 -D{r 12 , z 20 ; u)D(r 2 i, z 20 ; u) (1.67) 

Jz a z 2 J 47rr 12 r 2 o 

The generating functionals given above are very useful because using them one can easily calcu- 
late the multiple inclusive distributions for gluons or dipoles. In fact after taking the functional 



derivatives with respect to u one has just to set its value to 1 instead of as in (1.65). On 
setting u = 1 without taking any derivative one has the probability to find any number of gluons 
or dipoles emitted with a scaling variable z in the range Zq < z < z\. These quantities for 
both ( 1.66j ) and (1.67) are evidently not normalized and divergent, the divergence due to the 
integration in the region of small inter-gluon distances in the transverse plane (rij = 0). 

In order to cure these ultraviolet divergences one has to include also the virtual contributions 
given by loop diagrams; in fact the previous computation was performed only at the tree level. 
One can even compute the contribution of the virtual corrections by imposing that the divergent 
part of real and virtual gluon emissions compensate each other to give a normalized probability 
and checking that the contribution is correct at the lowest loop order. This procedure has been 



followed for example in [11] and leads to the following equations for the functionals Z and D: 

r 02 ' 2 ! 



rrl \ 2w(r 10 ) . 2 f Zl dz 2 f d r 2 2u{r 10 ) s 

(1.68) 

n / \ i \ Mrw) , 2 f Zl dz 2 f d 2 r 2 2w(ri ) n / \ 

D(r 10 ,z 10 ;u) = u{r 10 )z 10 K '+ cr w /- — 3— 3-^2 ' D(r u , z 20 ;u)D(r 20 , z 20 ;u) 

Jza z 2 J 47rr 19 rr, n 



(1.69) 

where the quantity u{r) has the same functional form of the gluon Regge trajectory respect to 
its momentum argument, i.e. 



»M = "( c / 2 >oi / = (2a s C F /ny 0l [ (1.70) 

J * 7r '02'21 J ^ 7r '02'21 



The virtual corrections to the functional equations are related to the property of gluon reggeiza- 
tion of the BFKL approach. 



Let us now look at the connection between the colour dipole model and the BFKL pomeron. 
One first obtains an equation for the inclusive one dipole distribution taking a functional deriva- 
tive with respect to u of relation ( 1.69| ) and setting u = 1. The resulting equation in term of the 
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rapidity y = In zio reads 

n(n , r, y) = e 2 ^ r ^5 2 (r - r w ) + 2cr2 Q f dy > e 2(y-y'Mno) f _ n ( r21 ,r, y') (1.71) 

and taking a derivative with respect to y one has the evolution equation 

S-n(r w , r, y) = 2u(r 10 )n(r 10 ,r, y) + 2ois< ^ F f <j 2 r 2 2 10 2 n(r 2 i,r, y) (1.72) 
dy tH J rf 2 r| 

which describes a dual BFKL diffusion in the transverse coordinate space instead of in the 
transverse momentum space. This equation is also infrared and ultraviolet finite but one may 
need to regularize it with an ultraviolet cutoff; in it the role of the coordinates and momenta is 
interchanged. 
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Chapter 2 

Hard pomeron with a running 
coupling constant: a 
phenomenological approach 



The basic elements in the BFKL approach for the multi-Regge kinematics of Yang-Mills theories 
are the effective vertices (real contribution) and the Regge trajectory (virtual contribution) 
whose separate bad behaviour in the infrared sector cancels in the total cross section. In the 
LLA approximation all the contributions, which are present in the BFKL equation for the 
pomeron, are based on calculations in the first non-trivial order in perturbation theory. The 
effective QCD coupling constant (taken to be small since a s In s ~ 1) is considered fixed and one 
does not know the region of applicability of LLA and in particular the energies and momenta 
fixing the scale for a s . This is at the basis of the weakness of the LLA approximation, since the 
quantitative results can be strongly modified by changing the scale of virtuality. 

For this reason a program to compute the next-to-leading order corrections to the BFKL 
kernel of the Lipatov equation has been carried out in the last years |6| with the hope to study 
soon the so called quasi-Regge kinematics for gluon production. The one loop corrections to the 
effective vertices, the two loop corrections to the gluon trajectory and terms for the productions 
of two gluons and quark-antiquark with fixed invariant mass have been computed [f7j. The nice 
feature is that also at the NLO order the pomeron can be seen as a compound state of two 
reggeized gluons. The generalized BFKL equation however is quite complicated [||] and there is 
no hope to obtain an explicit analytical solution but the study of such a kernel should be very 
useful to understand by means of numerical calculations the validity of the predictions in the 
Regge kinematics and corrections in the quasi-Regge limit. Moreover it is still not clear if such 
results are applicable in the HERA kinematical region or one should go even beyond the NLO 
analysis ||. 

A part of this well established program some attempts have been made in order to take into 
account the running of the QCD coupling with a minimal added work. The first idea 27] has 
been to introduce directly in the BFKL equation for the pomeron the running coupling with 
some assumptions about the momentum and distance scales at which the gluons interact. Also 
in other attempts 47, 4^, 48] the BFKL equation has been modified in the infrared sector 
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in order to model some non-perturbative effects and at least try to avoid the strong inconsistency 
one has using the BFKL evolution which does not separate perturbative and non-perturbative 
scales. This fact is evident on considering the diffusion mechanism of the gluon distribution 
in the transverse momentum space during the evolution in the rapidity leading to a typical 
normal distribution in Ink 2 . In all these phenomenological models a better behaviour also due 
to some decrease of the pomeron intercept to more reasonable values was found, but one of 
the fundamental properties shown by the QCD in the Regge kinematical limit, that the gluon 
reggeizes, is lost. 

The idea to introduce a running coupling in a manner consistent with the gluon reggeization 
has been proposed by M.Braun |l3|, [14]] and in this way one is able to guarantee that the 
production amplitude in the one reggeized gluon exchange approximation, which serves as an 
input in the BFKL theory, satisfy unitarity in the leading order p^] . Thereby the whole scheme 
becomes self-consistent: otherwise one should add to the input amplitudes corrections following 
from the unitarity. Since the bootstrap condition, as discussed in section 1.x, is a crucial 
element for the reggeization of the gluon and for the theory of the reggeized gluons as a whole 
the only way to introduce the running coupling constant in a manner compatible with the gluon 
reggeization is to preserve the bootstrap. A second basic ingredient is the requirement that the 
outcoming gluon distribution has the same asymptotic behavior as the one predicted by the 
DGLAP equations (describing the renormalization group flow) in the DLLA. 



In section 1 we present the basic equations in the form suitable for numerical analysis 15] 
for the cases q = (forward scattering) and q ^ 0. After that we describe the method of the 
solution and present the numerical results for the intercept, slope and the wave function at q = 0. 
The results for the intercept and slope, on the whole, agree with those found in [[0], 14] by the 
variational approach. An interesting new result is the existence of a second pomeron with the 
intercept minus one, roughly speaking, two time less than for the leading one, but still positive. 
We end this section applying these results to the study of the asymptotical behaviour of the 
cross-section for the 7*7* scattering, including some unitarization effects due to multipomeron 
scattering by mean of an eikonal formulation valid in the large iV c limit ]j5C[] , and try also to 
extrapolate to the hadron case. 

Therefore the results of our study show that with the running coupling included the pomeron 
equation possesses bound state solutions which have negative energy and thus intercepts greater 
than unity. These solutions correspond to supercritical pomerons in the old sense, that is, they 
represent simple poles in the complex angular momentum plane. However the sub-dominant 
pomeron does not seem to play an important role in describing the asymptotical behaviour of 
the amplitudes. The intercept of the leading pomeron singularity depend weakly on the infrared 
regulator parameter and stay in the region 0.35-0.5 for its physically reasonable values. The 
introduction of the running coupling and thus a scale provides for a nontrivial slope for the 
pomeron, which is responsible for the physically reasonable behaviour of the cross-sections at 
very high energies. 

For realistic photonic cross-sections and, with a rather crude approximation, for the hadronic 
ones the estimated unitarization effects begin to be felt at extraordinary high energies, of the 
order 100 — 1000 TeV (or equivalently x < 10~~ 10 — 10~ 12 ). Until these energies a single pomeron 
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exchange remains a very good approximation to the asymptotic amplitude. By the way we 
remember that no effects due to triple pomeron interaction are here taken into account and 
these are seen to be required as shown by the unitarization program (see the third chapter). 
Assuming that the picture is not completely changed by the last consideration, a comparison to 
the experimental cross-sections and structure functions at the highest energy (lowest x) achieved 
seems to confirm the widespread opinion that we are still rather far from the asymptotical regime 
and that other states, different from the supercritical pomerons, give the dominant contribution. 

In section 2 we apply our model with the two found supercritical pomerons to jet production 



16]. This process has been extensively studied in the framework of the standard BFKL approach 



with a fixed coupling [pl| , 45]. As is well-known, this analysis has lead to some far-reaching 



conclusions as to the importance of mini-jet production at high energies and the logarithmic rise 
of the multiplicity. However in the fixed coupling approach, various ad hoc modifications of the 
basic BFKL model had to be introduced to cut off the spectrum at low k± and also to correctly 
reproduce the high k± tail of the spectrum. Both these problems are naturally resolved by the 
introduction of a running coupling in our model, in which no new parameters appear in contrast 
to the fixed coupling approach. We shall study the asymptotic behaviour of the jet production 
cross section using the asymptotic form of the pomeron wave functions shown in appendix A.l 
and present some numerical calculations. 

Since we feel that to describe the present experimental data it is necessary to take into 
account all the states from the the spectrum of pomeron equation, this one has been converted 
into an evolution equation in 1/x and solved with an initial condition at some (presumably 
small) value x = xq. In such an approach [17], taking a non-perturbative input at x = xq 
adjusted to the experimental data, also the problem of coupling the pomeron to the hadronic 
target is solved in an effective way. This approach is shown in details in section 3 where we 
state the basic equations and we try to pass from the gluon density to the observable structure 
function. 



2.1 Pomeron equation with running coupling 

Let us consider the bootstrap condition for the reggeization of the gluon: 

[co(q) -co( qi ) - wfeMfi) = J ^K^q'^qi) (2.1) 

This is in general a very complicated functional integral equation which relates u, K and <f>. 
Looking for a particular class of solutions of this equation (we look for u> and K, limiting 
ourselves to constant </>) we can start to write them making use of more than one arbitrary 
function. In fact we can write the reggeized gluon trajectory 



(2.2) 
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and the gluon pair interaction kernel in the gluon channel 

u( /n _ N c f jA(qi) ./UfeK 1 fA(q) \ ( , 

One must note that there are some constrains on the kernel K which must be satisfied. The 
first one is that K should be symmetric in the two gluons (indices 1 and 2) so we must require 
that fs = fc from the last term in the kernel. The second requirement is that K should be 
symmetric in the initial and final gluons (quantities with momentum not primed and primed), 
after separating the two "gluon propagators" ( f B ( q l )f B ( q 1 ) )- From the first term in the kernel 
we must thus require this symmetry for the expression 

fA{qi)fB{q' 2 )+fA{q2)fB{q' 1 ) , , 

fc(qi -q[) { ' j 

Matching the expression with the corresponding one where the primed quantities are exchanged 
with the not primed we finally get = Jb- So we arrive at /a = /b = /c = V- The standard 
BFKL theory is characterized by r\ = q 2 /(2a s ) and in general the gluon trajectory and the 
interaction kernel are correlated in terms of the rj function. 



Let us write the equation for the amplitude in the vacuum channel in the form from which 
the typical cancellation of the infrared divergences in the pure BFKL theory is manifest 

(j - mqi ) = N c J wy2 [^-^ ( m - v{q[) + r]{qi _ q[) ) + (1 - 2) - ^^r^l 

(2.5) 

The idea in 14 1 was to change 77(g) so that it corresponds to a running rather then to a 
fixed coupling. For the running coupling some conclusions about the form of rj(q) can be made 
considering the vacuum channel equation in the limiting case of very large q. In this case one has, 
assuming that r](q) grows with q, the approximated equation in the forward direction (q = 0) 

and in the leading approximation one gets, after an inverse Mellin transformation, an evolution 
equation for the fully amputated wave function 

" 3 -(x, qi ) = ^^-<P(x, qi ) (2.7) 



91n^91nl/x 2ir r](qi) 

which, on comparing with the DGLAP evolution equation in the leading order in lnl/x (that 
is, in the double leading log approximation (DLLA)) for the gluon distribution, 

i} 3 - xg ( x , qi ) = N c ^®xg(x, qi ) (2.8) 



d In qf d In 1/ x tt 
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allows one to find the asymptotic form of 77(g) : 

* (?) *2^)' q ^°° (2 ' 9) 

and a proportionality relation between the amputated pomeron wave function (j)(x, k 2 ) and the 
gluon distribution xg(x,k 2 ). It evidently differs from the fixed coupling case by changing the 
fixed coupling constant a s to a running one a s (q 2 ). As a result, with a running coupling, both 
the gluon trajectory and its interaction have to be changed simultaneously in an interrelated 
manner, so that the resulting equation is different from the BFKL one already in the leading 
order. 

The behaviour of rj(q) at small q, comparable or even smaller than the QCD parameter A, 
cannot be established from any theoretical calculation, since this domain is non-perturbative. 
To take into account these confinement effect we choose rj(q) at finite q in a simple form: 

r ] (q) = (b/2n)f(q) ; f(q) = (q 2 + m\) ln((q 2 + m 2 )/A 2 ) (2.10) 

with b = (1/4)(11 — (2/3)Np) and m > A, which agrees with ( p.Sj ) for large q and remains finite 
up to q = 0. It allows for the freezing of the coupling and the confinement proper to occur at 
somewhat different scales (m and m\ respectively). However, on physical grounds, one feels that 
they should be of the same order. 

A preliminary study of the properties of the pomeron with 77(g) given by ( |2.10| ) for m = m\ 



was performed by the variational technique in 13, 14]. It was found that the intercept depended 



on the ratio m/A quite weakly: as m/A changes from 1.5 to 5.0 the intercept (minus one) A 
falls from 0.4 to 0.25. On the other hand, the slope depends on this ratio very strongly. This 
allows to fix the ratio m/A to values in the interval 3.0-^4.0. 

This variational study, although very simple, cannot however give values for the intercept 
and especially for the slope with some precision. Still less can be found by this method about the 
properties of the pomeron wave function essential for the high-energy behaviour of the physical 



amplitudes, for some aspects studied in ||50| . Finally, one does not receive any knowledge about 
the existence of other solutions with a positive intercept. All these reasons give us a motivation 
to undertake a numerical study of the two-gluon vacuum channel equation with in general the 
gluon trajectory and interaction given by 

u(q) = f ^ "® 



2 J (2tt) 2 77(^)77(5; 



2) 



K (a n>\ TT (M®-) I ^ 2 \ 1 

\ vWi) »?wi-9i) n{Qim%)J 

where T is the gluon colour vector (T1T2 = — N c in the vacuum channel) and the function r](q) 



satisfying (|2jj|) . 

The parametrization ( |2.10| ), with b fixed, does not evidently take into account the change 
in Np which occurs on flavour thresholds. In principle this can easily be remedied by assuming 
Np and thus b depending on q 2 and varying accordingly. Our theory possesses a scale (actually 
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several scales, m,mi and A), so that the asymptotic behaviour (5) used to determine r/ refers to 
q 2 » m 2 , m 2 ,A 2 . Since m ~ mi is of the order of 1 GeV (the actual choice we make in following 
is m = mi = 0.82 GeV), the only thresholds to be accounted for are those of the heavier quarks, 
b and t. However one should take into account that our theory is much better behaved than the 
scale-less BFKL model. In particular, in our case the pomeron possesses a normalizable wave 
function. So, in contrast to the BFKL case, the characteristic momenta which form the pomeron 
are determined by the same scale m ~ mi and are therefore of the order of several GeV jc. For 
this reason we expect to find the pomeron properties with sufficient precision fixing Nf = 4 and 
forgetting about b and t thresholds. Accordingly, the calculations presented are performed with 
a constant Np = 4. Clearly the found momentum wave function is strictly speaking valid up to 
momenta squared of the order 100 (GeV/c) 2 , when the beauty threshold should be taken into 
account to slightly lower the value of our parameter b (by 8%). In applications to the structure 
functions this restricts the domain of applicability of the calculated pomeron wave function to 
Q 2 < 100 (GeV/c) 2 . 



2.1.1 Basic equations. Pomeron at q = 

We consider the physical case N c = 3. The units are chosen to have A = 1. In relating to 
observable quantities we take A = 0.2 GeV. 

The pomeron equation is the eigenvalue equation 

(-u J (q 1 )-u{q 2 ))<t>(q 1 ) + / ^K^(qi, q ' 1 )<p( q [) = E( q )^ q i) (2.12) 

where the "energy" eigenvalue E(q) is related to the pomeron trajectory via 

a(q) = 1 - E(q) ~ 1 + A — aq 2 (2.13) 



The last equation, well satisfied for small q, defines the intercept A and the slope a'. In ( p. 12 ) 



the trajectories u and the kernel K vac are given by the eqs. ( p. 11 ) with T1T2 = —3 and the 



function rj given by ( p. 10 ). To symmetrize the kernel we pass to the function 



HQi) = <K<Zi)/vWM92) (2-14) 

We also take out the common numerical factor 6/((ll — 2/3Nf)tt) and express all terms via the 
function f(q) also defined in ( p. 10 ). Then the equation for if) takes the form 



A q (qi)^(qi) + J d 2 q[L q ( qi ,q[)^(q[) = e{q)^{qi) (2.15) 
Here the "kinetic energy" is 

Mqi) = (V2) J fWi)f(qi-q'i) + (V2) J (2 - 16) 



The interaction kernel consists of two parts, a quasi-local and a separable one: 
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They are given by 



L q 



/fa) 1 Ml IM 1 r9 17 n 

/(<&) f(Qi - Q'l) V V f(Qi) f(Q2 ~ Q' 2 ) V /(«£) j 1 J 



and 

4 Sep) (^?i) = / = (2-18) 



Both parts are evidently symmetric in q\ and q[. The scaled energy e is related to the initial 
one by 

E = j( u-(W) e (2 ' 19) 



eq. ( |2. 15 ) simplifies in the case when the total momentum of the two gluons is equal to zero. 



With q = the two parts of the kinetic energies become equal and the square roots in (2.17) 
turn to unity. So at q = the equation retains its form ( p. 15 ) with 



f(q[)f(qi-q[) 
and the interaction where now 

4*Wi) = - ,, 2 .s (2-21) 

j{qi - <?i) 

has really become local and 

This is the equation which we shall solve numerically. 

To reduce to a one-dimensional problem we introduce the angular momentum of the gluons 
n and choose the solution in the form 

i>{q) = ^n{q 2 ) exp in<f> (2.23) 

where (ft is the azimuthal angle. Integrating over it in the equation, we obtain an one-dimensional 
integral equation for the radial function ip n (q 2 )'- 

A Q (q)Mq 2 ) + J dqlL n {q\ql)^{ql) = e^ n (q 2 ) (2-24) 
with the kernel now given by 

L n (q 2 , q\) = -B n (q\ql) + <5 n0 7T /(0) (2.25) 

f(q)f(qi) 

where 

B n {q\ql)= d<ft f 2.26 

Jo f{q 2 +q(-2qq l cos (ft) 
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Note that Aq can be expressed via Bq: 

Mi) = (1/2) / d q fB ( q 2 ,ql)M^ (2.27) 

Evidently eq. Q2.24|) is very similar to a Schrodinger equation with an attractive interaction 
provided by the local term and a positive kinetic energy described by A, which however grows 
very slowly at high momenta (as lnlng). Therefore the attraction becomes smaller with growing 
n and we expect to find negative energies, corresponding to intercepts larger than unity, only 
for small n. Remember that for the BFKL pomeron only the isotropic state with n = has a 
negative energy. Our calculations reveal that the introduction of the running coupling following 
does not change this situation: states with |n| > all have positive energies. So in the following 
we consider the case n = 0. 



2.1.2 Pomeron at q 0: the slope 

With q 7^ the pomeron equation becomes essentially two dimensional. Rather than attempt to 
solve it numerically at all q we limit ourselves to small values of q and determine not the whole 
trajectory a(q) but only the slope a' defined by ( p. 13 ). This can be done in a much simpler 



manner using a perturbative approach. We present "the Hamiltonian" in ( |2.15| ) H q = A q + L q 
in the form H q = Hq + W(q) and calculate analytically W(q) up to terms of the second order 
in q. Then for small q the value of the energy e(q) will be given by the standard perturbation 
formula 

e(q) = e(0)+ < W(q) > (2.28) 

where < > means taking the average with the wave function at q = 0, determined from the 
numerical solution of the equation discussed in the previous section. Thus we avoid solving the 
two-dimensional problem, but, of course, cannot determine more than the slope. Fortunately it is 
practically all we need to study the high-energy asymptotics (although, of course, the knowledge 
of the trajectory as a whole might be of some interest). 

In order to derive an expression for W(q) we pass to the relative momenta I and I' q\(2) = 
(l/2)q + (— )l and similarly for the primed momenta. Up to the second order in q we have 

f(qi) = /(0[1 +ai(q-l) + ^q 2 + y (q ■ I) 2 ] (2.29) 

where 2 2 

ai = Ol (0 = [1 + ln(/ 2 + m 2 ) - ^ l^ JtT) (2 - 30) 

1 to 2 — m\ 1 

02 = a2{l) = 75 o + 77T5 ovjItTTT (2-31) 

w L P+m 2 (I 2 + m 2 ) 2 f{l) K ' 

The expansion for fijq^) differs by changing the sign of I (or of q); for f(q[) and /(q^) it suffices 
to replace I by I' in the expressions for f(q%) and /(<?2)- We use the notation a[ = ai(l') and 
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0-2 = a 2(l')- We also need the expansion for f(q): 
f(q) = f(0)(l+a 3 q 2 )+O(q 4 ) 



03 



mj 1 2 \ 1 
-rr + mm 



/(0) 



(2.32) 



The perturbation W(<?), up to second order in q, can be expressed via the introduced functions 
a\ 2 and 2 and the constant 03. After some calculations we find a part of W coming from the 
kinetic term in Hamiltonian in the form 



wi(0 



f(i')^f(i-n /( 

aia'^q ■ l)(q ■ l') aia[(q ■ l)(q ■ I') 



fii-v) 



+ 



n-i - n 



}+Ao(Z)[f 9 2 + f(,-0 5 



The part of coming from the quasi-local part of the interaction can be written as 

W2(l,l') = \[a l (q-l)-a' 1 (q-l')} 2 L^\l,l') 
and the one coming from the separable part as 

W 3 (l,l') = [(a 3 - a -±±A)<? - h a2 - a 2 )(q • if - \(a' 2 - a'\)(q ■ l')*]4"*\l,l') 



(2.33) 



(2.34) 



(2.35) 



As mentioned, only isotropic solutions have the intercept larger than one and are of interest. 
Then the expression for W(q) = Yli=i 2 3 nas to be integrated over the azimuthal angles. Thus 
integrated values will be denoted Wj, i = 1,2,3. Using ( |2.25D - ([2.27 ), they can be conveniently 
expressed via the kernel B n (eq. ( p. 26 )): 



Wx 
2^2 



\ J dl,2 j^{i-Y + ( a 'i ~ Y^ Bo ^ " <*<W B & 0} + + «2/ 2 ) (2.36) 



w 2 

2nq 2 ~ 
/(0) 



a^ 2 + a'll' 2 )B (l, I') + oioi«'Bi(l, /') 



W3 



2V /(/)/(/') 



ai + a'i . 1 



27TI „ ( - — -p-L. , - ^/ 2 (a 2 - a 2 ) - V(a 2 - a' 2 ) 



(2.37) 



(2.38) 



The slope is given by the momentum average of the sum of these expressions, taken with a given 
isotropic wave function: 



(1/2) j di 2 di ,2 ^ii)^{i')(w 2 + w 3 ) + j d/vco 2 ^! 

2V/d/ 2 V' 2 (/) 



a 



(2.39) 



2.1.3 Numerical procedure and results 



Eq. ( [2.24 ) was first changed to the variable t = lnq 2 whereupon the wave function and the 
kernel transform according to 

HQ 2 ) - m = q^(q 2 ) (2.40) 
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Figure 2.1: .Pomeron intercepts as a function of the infrared regulator mass m = mi; 

A = 0.2 GeV (a) and for different values of the confinement parameter mi 
and the coupling freezing parameter m (b). 



and 

L(q 2 , q\) - L(t, h) = qq X L{q 2 , q\) (2.41) 

Performing a discretization on a grid the equation was reduced to a finite system of linear 
equations by approximating the integral by a sum 

/oo n 

dtF{t)~Y. W i F ^) ( 2 - 42 ) 

-°° i=l 

with the grid points ti and weights it?j depending on the chosen approximation scheme. The 
final equation is thus 

n 

By-ay = ea^i, i = 1, ■ (2.43) 

where 



and 



(2.44) 

By = + y/nHwJLit^tj) (2.45) 

The vector space for x has been chosen in order to obtain a symmetric matrix B. After de- 
termining the lowest eigenvalues e and the corresponding eigenvectors x% the wave function in 
the momentum space is directly given by fl2~4C| ) and (fOl at points <? 2 = exptj. It should be 
normalized according to 

l^WI 2 = 1 (2-46) 

Note that this wave function is a partially amputated one (see eq. ( [2.14D ). The full (non- 
amputated) wave function is given by $((/) = ip(q)/rj(q) (in the forward case). It is this function 
that appears in the physical amplitudes. 

The results for the lowest (and negative) eigenvalues of energy for the case n = (isotropic 
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Figure 2.2: Pomeron slopes as a function of the infrared regulator mass m — mi; 
A = 0.2 GeV. 



pomeron) are presented in Fig. 2A. Actually energies with an opposite sign are shown, which 
according to fl2,13j ) give precisely the intercepts (minus one). As mentioned, the QCD scale here 
and in the following is taken to be A = 0.2 GeV . In Fig. |2.1|.a the intercepts are shown for the 



case when the two scales m and m-i in ( |2.10| ) are equal. Fig. ^.b illustrates the dependence 



of the intercepts on the ratio m/m\. The most interesting observation which follows from these 
figures at once is that in all cases one observes two positive intercepts, which correspond to two 
different supercritical pomerons, the leading and subleading ones. The single cut, characterizing 
the BFKL spectrum, is deformed with the introduction of the running coupling and a discrete 
part of the spectrum, giving some poles, appears. The intercept of the leading pomeron is found 
to be in accordance with the earlier calculations pjl, 14], performed by the variational method 



(it is slightly larger, which was to be expected). For physically realistic values of m and m\ in 
the interval 0.5 -j- 1.0 GeV it takes on values in the region of 0.5 -j- 0.3 falling with the masses m 
and mi. The same trend is seen for the intercepts of the subleading pomeron, which lie in the 
interval 0.25 4- 0.15. 



The slopes of the two found pomerons are shown in Fig. 2.2 as a function of m for the case 
m = mi. The slopes depend very strongly on the value of the regulator mass. The physically 
reasonable slopes for the dominant pomeron of the order of a' ~ 0.25 (GeV/c)~ 2 restrict the 
values of m to the region 0.7 — 0.9 GeV. So finally we choose 

m = 0.82 GeV (2.47) 

which leads to the following parameters of the leading (0) and subleading (1) pomerons 

A = 0.384, a' = 0.250 (GeV/c)- 2 ; Ai = 0.191, a[ = 0.124 {GeV/cy 2 (2.48) 



In Fig. 



2.2 



we show the coordinate space wave functions $(r) of these two pomerons. 
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Figure 2.3: Coordinate space wave functions for the leading ($ (r)) and subleading 
($i(r)) pomerons. Both r and 5> are in units 1/A ~ 1 fm. 



2.1.4 Pomerons and the high-energy scattering. 

To apply the found results to the actual physical processes one has to couple the pomerons to the 
external sources corresponding to the colliding particles. The only way to do it in a more or less 
rigorous manner is to assume that both the projectile and target are highly virtual photons, or 
heavy onia, with momenta q and p respectively, — q 2 = Q 2 >> A 2 and —p 2 = P 2 » A 2 . Then 
the non-perturbative effects inside the target and projectile can safely be neglected. We shall 
address to the very complicated problem of coupling the gluons to the virtual photons in the 
running coupling framework in section 2.3 and consider here the simple lowest order approach, 
sufficient for the BFKL pomeron, where the coupling is small and fixed. This corresponds to 
taking the contribution of a pure quark loop into which the incoming photon goes. We shall then 
need of the colour densities of the virtual photons whose explicitly form, in the forward case, 



was found in [46] for both transverse and longitudinal photons. In the appendix A. 2 we present 
a generalization for the non-forward case, both in momentum and coordinate representation, 
since we shall need of it, at least formally, in the third chapter, where we return to the fixed 
coupling formalism. 

We shall try in the following to understand the role of the two found pomeron states and 
analyze the situation in which more than one pomeron is exchanged. The last situation can be 
seen as an attempt to understand some aspects which are related to higher order corrections 
needed to restore the unitarity. Incidentally the true situation is more complicated, as is shown 
in chapter 3, since, for a two pomeron amplitude, one has to take into account the triple pomeron 
interaction, which can take a form which absorbs completely the double pomeron exchange. 

Nevertheless the contribution of multipomeron scattering to some processes is of interest 



and we shall try to analyze it as done in [5C] where the scattering amplitude in the large colour 



number limit takes an eikonal form for fixed transverse dimensions of the projectile and target 
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Figure 2.4: Structure functions of a virtual photon (P = 1 GeV/c) at Q 2 = 4 and 
100 (GeV/c) 2 as a function of x (a) and of the proton at Q 2 = 4 and 
100 (Gel//c) 2 as a function of a; (b) (solid curves). Dashed curves show the 
contribution of a single pomeron exchange. 



and leads to a cross-section 

a = 2 J d 2 Rd 2 rd 2 r'p q {r)p p {r') (l - e -*<rAr,^ (2.49) 

where 

z(u,R,r,r') = (1/8) / ^ gfo g, gi, g!)^* (1 - e^) (1 - e^') (2.50) 



i=l,2 



is essentially the Fourier transform of the (non-amputated) Green function of eq. (2.12), 
G{y, g, gi , q[ ) , considered as a function of the energetic variable u = pq and with g = gi + g2 = 
q[ + g2- The functions p q and /3 p correspond to the colour densities of the projectile and target 
photons, respectively. 

The found supercritical pomerons represent a part of the total pomeron spectrum which 
contributes to the Green function, in the high energy limit, with a term given by 

Gp(^g,gi,gi) = E ^ l(<?) " 1 ^(9i,g 2 )^(gi,g^) (2.51) 

i=0,l 

where and $j are the trajectories and wave functions of the leading (0) and subleading (1) 
pomerons. At high v we can neglect the dependence on the total momentum g of the wave 
functions, taking them at g = 0, and approximate the trajectories according to ( p.13 ). Then 



all the quantities in ( |2.49| ) become determined, so that we can calculate the cross-sections for 
both the transversal and longitudinal projectile photon and thus find the structure function of 
the virtual photon target. We have taken for the latter a transversal photon with the lowest 
momentum admissible of P = 1 GeV/c. The resulting structure functions are presented in Fig. 



2. 4. a for the interval of small x which we extended to extraordinary small values to clearly see 



the unitarization effects. 
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To move closer to reality one has to consider hadronic target and projectiles. The confinement 
effects make any rigorous treatment of such a case impossible. Rather than to introduce arbitrary 
parameters (in fact, functions) we extend the formula (2.49) to hadronic target and projectile 
substituting the photonic colour densities by hadronic ones. For the latter we take a Gaussian 
form and a normalization which follows from the comparison to the electromagnetic densities 
with only the simplest quark diagrams taken into account. In particular for the proton we take 
the Gaussian p, with the observed electromagnetic proton radius and normalized to three. Such 
a treatment, in all probability, somewhat underestimates the density, since it does not include 
coupling to gluons. 

The proton structure functions and the proton-proton total cross-sections which follow from 
this approximation for the densities are shown in Figs. ^J.b and [T^ respectively. To see the 
unitarization effects we had again to consider extraordinary high values of 1/x and energies, well 
beyond our present experimental possibilities. 

In discussing these results, we have first to note that their overall normalization is somewhat 
undetermined, since the exact scale at which v enters into lnz^ factors is unknown. A second 
point to note is that the subleading pomeron contribution is always very small: it amounts to a 
few percent at smallest values of 1 jx and s considered and naturally gets still smaller at higher 
1/x or s. 

As one observes from Figs. ^4-2Jj, the structure functions and cross-sections monotonously 
rise with 1/x, s and Q 2 . Studying the asymptotics of the solutions of eq. ( [2.12 ) at high q and of 
eq. ( 2.49| ) one can show that this rise is logarithmic. In particular, the structure function of the 
virtual photon rises as ln 4 (l/x) and as ln^(q 2 ) with f3 ~ 2.5. The proton-proton cross-section 
eventually rise as In 2 s, as expected. Comparison to the Froissart bound (dash-dotted line in 
Fig. 7) shows however that it remains far from being saturated. 



The most interesting result that follows from Figs. |2.4]2.5| is that the unitarization effects 
become visible only at exceedingly very small values of x or very large values of s, well outside 
the range of the present experiment. They appear earlier at lower Q 2 . Still at the smallest value 
Q = 2 GeV/c considered, the exchange of more than one pomeron achieves only 15% of the 
total for the proton structure function at x = 10~ 10 . Likewise the relative contribution of many 
pomerons to the proton-proton cross-section rises to 23% only at yfs ~ 10 5 GeV . 

Comparing the calculated proton structure functions and the cross-sections with the ex- 
perimental results at highest 1/x and s available we observe that our results are essentially 
smaller than the observed ones. Experimental value of F2 P (Q 2 ,x) at Q 2 = 8.5 (GeV/c) 2 and 
x = 0.000178 is 1.19 ± 0.05 ± 0.16 [j52j| . Our calculations only give a value 0.17. The pp cross- 
section at ^fs = 1800 GeV is around 80 mbn 53], whereas our result is 18.5 mbn. Of course, 
having in mind the uncertainties in the overall normalization and a very crude picture for the 
pomeron coupling to the proton assumed, one cannot ascribe too much importance to this fact. 
However one is tempted to explain this underestimation of the experimental values by the sim- 
ple reason that we are too far from the pure asymptotic regime yet and that other solutions of 
eq. ( |2.12| ) different from the found supercritical pomerons and having their intercepts around 
unity give the bulk of the contribution at present energies. This may also explain the notorious 
discrepancy between a high value of the hard pomeron intercept, of the order 0.35-0.5, and the 
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Figure 2.5: Proton-proton total cross-sections as a function of cm. energy yfs (the 
solid curve) . The dashed curve shows the contribution of a single pomeron 
exchange. The dash-dotted curve marks the Froissart bound. 



observed slow growth of the experimental cross-section, well described by the "soft pomeron" 
with an intercept around 0.08. 

Forgetting for the moment the triple pomeron interactions, if this picture is correct then we 
may expect that with the growth of energy the cross-sections will grow faster and faster, until 
at \fs ~ 10 TeV they will become well described by a pure hard pomeron with the intercept 
0.35-0.5. This power growth will continue until energies of an order 1000 TeV when finally the 
unitarity corrections set in to moderate the growth in accordance with the Froissart bound. 



2.2 Inclusive jet production. 
2.2.1 General formalism. 

We shall consider in the following the scattering of two highly virtual photons (or heavy " onia" ) 
to make the derivation more rigorous. To obtain the formula for the inclusive jet production let 
us recall the total cross section for one pomeron exchange pOfl: 



o = \j <Pr<Pr'p q {r)p p {r') J f^f G (u,0,qi,q[) JJ (1 - e^)(l - e^') (2.52) 

Here q(p) is the momentum of the projectile (target); v = qp = (l/2)s; G(v,q,qi,q[) is the 
non-amputated pomeron Green function in the momentum space, q being the total momentum 
of the two gluons and q\ (q[) being the initial (final) momentum of the first gluon; p q [ p ) is the 
dipole colour density of the projectile (target). For the virtual photons we take the colour dipole 
densities from [Hj. 

In the high energy limit the dominant term of the Green function comes from the two 
mentioned supercritical pomeron states ( |2.51[ ) where on and <3?j are the trajectories and wave 
functions of the leading (0) and subleading (1) pomerons. 

Mini-jets appear as intermediate gluon states in the Green function ( |2.51| ). They possess 
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=> 



Figure 2.6: The substitution for the Green function used to calculate the inclusive jet 
production cross section. 



arbitrary k± subject to condition k\ « s. The inclusive cross-section for their production can 
be calculated by splitting the Green function in ( |2.52 ) as indicated in Fig. 2J: whihc means that 
we perform the substitution 



G(u,0,q 1 ,q' 1 ) 



d 2 k 1 d 2 k[ 



G(y u 0, qi,ki)V(kiX)G{v2,0, k{,^)S^ (h +k±- k[) (2.53) 



where, for the running coupling case, 

Remembering that q\ + qi = q[ + q' 2 = q = we obtain 
d 3 a 



(2.54) 



dyd 2 k A 



d 2 rd 2 r' Pq (r)p p (r') [ ^f/ 1 (1 - e^)(l - e*^ f 



(2tt)' 



(2k)* 



5 (2) {kl + k ± - fci)G(yi,0, ?1> fc!)G(z/ 2) 0, k[,q[)[ ^ )v{ f) - 3*7(0)] (2.55) 



The appropriate kinematical variables are the rapidity and transverse momentum of the 



observed gluon (jet), y = ^ In ^ and k±, respectively ; we also have Si = 2v, L and su 2 ) 
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kj_^/se It is convenient to use a mixed (momentum-coordinate) representation for the 

Green functions: 



G(u,0,r,k) = J ^ e ^ r G(u,0, gi ,k) 
Defining A$i(r) = $j(r) — 3>i(0) we can write 

G(i/,0,r,&) - 0, 0, Jfe) = ^ ^ Al A^(r)$*(A:) 



(2.56) 



(2.57) 



i=0,l 



We also introduce 

= (A*i) ff(p) = | dV^^A^r) (2.58) 
In this notation and also using both the semi-amputated and full wave functions we find 

I(y,k ± ) = £ fl?J#^ /" d " kl \M^)^x + k ± ) 

i,i=0,l 



(2tt)^ 



r?(/cj 



fyCO^CfcO^-Cfci + fcj.) (2.59) 



Since the found supercritical pomerons are isotropic in the transverse space the inclusive 
cross-section I(y,k±) also turns out to be isotropic and therefore we can integrate over az- 
imuthal angles in ( p.59[ ). Defining for the full and semi-amputated wave functions the integrated 
quantities 



Hh,k ± ) 



2- 



2- 



da$(ki + k±) 

dar](ki + kx)<S>(ki + k±) 



(2.60) 



where a is the angle between k\ and k±, we finally obtain the inclusive cross section 



<Fa 



dydk 



£ 

i,j=0 



4 



d/cf r 2 ^(fci)^j(fei,fcj 



(27T) 



T/(0)*i(fci)*,-(A!i,fei.) 



(2.61) 



The results of numerical calculations of the cross-section ( |2.6l| ) and also its generalization to 
more interesting cases of hadronic targets or/and projectiles will be discussed in the next sec- 
tion. In the rest of this section we shall study the asymptotic behaviour of the found inclusive 
jet production cross-section at very small and very large transverse momenta and also its y 
-dependence. 

As to the latter, all y-dependence in ( p. 61 ) comes from the factor exp(— y(Aj — Aj)) which has 
its origin in the existence of two different supercritical pomerons. Evidently in the limit s — > oo 
this dependence dies out, since the relative contribution of the sub-dominant pomeron becomes 
negligible. The model thus predicts an asymptotically flat y plateau at very high energies. 

At small k± the cross-section ( [2.61 ) evidently goes down as k^ 1 , since all other factors are 
finite in this limit. However one should remember that ( |2.51| ) gives the dominant contribution 
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only while (kj_s) continues to be large. At too small k±, when the above quantity becomes finite, 
all other states from the spectrum of the two-gluon equation ( p. 12 ), hitherto neglected, begin 
to give comparable or even dominant contribution, so that the found k^ 1 behaviour ceases to 
be valid. 

To find the asymptotic behaviour of the inclusive cross section for k± — * oo we need to know 
the behaviour of the pomeron wave functions in the momentum space at q — * oo and in the 
ordinary space at r — > 0. In the appendix A.l we show that 



ip{q) 
ip{r) 



q — >oo 



i 



~ (In-) 

r^O r 



13+1 



(2.62) 



" so that in the forward scattering case depends just on the intercept of 



bE 



where (3 = — 

the corresponding pomeron state. Let us study now the behaviour of (2.59) for k± — > oo. Using 
&(q) = ip(q)/'n( ( l) ~ ^(hi^ 2 )^ 1 we find / d 2 q<fr(q) < oo, so that for the second term in the 
integrand of ( |2.59| ) we get 



[<Pk 1 $i(ki)$ j (k 1 +k x ) — > $j(k ± ) fcPh^ih) ~ i(lnfc 

J k i —too J K I 



(2.63) 



To analyze the first term we use the identity 
d 2 L 



(27T) 



^rlf>i(kl)ll)j(kl + k A 



2tt 



rdrJ (k±r)i()i(r)ijjj(r) 



and the relation 



which is valid for naturally behaved f, such that rf{r)J\{k±r) 



poo 2 f°° 

I rdrJo(k±r)f(r) = / rdrJi(k±r)f'(r) 

o ~ k± Jo 



r=oo 
r=0 



0. 



Putting /(r) ~ (lni)^ l+ ^ j+2 we obtain the leading behaviour 



d 2 k^ i (k 1 )ifj j (ki + fcjj ~ -^(lnk 2 ^ 



So at k± — > oo we find for the inclusive cross-section ( [2. 59 ) 



I(y,k ± )~ R l R % 
j,i=o,i 



kjs. A i+ A o — V (A-—A -) 1 A ,2 



(2.64) 
(2.65) 



(2.66) 



(2.67) 



This asymptotics corresponds to the standard quark-counting rules behaviour (~ l/k]_), modi- 
fied by a power factor due the pomeron energetic dependence and a logarithmic factor coming 
from the pomeron wave function, i.e. from the running of the coupling. Note that in ( 2.6 7] ) the 
k± and s dependencies are separated. As a result we find that the average (k±) is finite and 
independent of s and ~ s A , since it formally diverges for ( |2.67| ) and one has to restrict 
k 2 , < s. 
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For the the multiplicity 



1 / A . . a , 



the standard asymptotic behaviour (n) = a Ins + b is obtained. Indeed integrating ( |2,59 ) we 



get / dyd 2 k± d yj± k± = Bs^° Ins + Cs A °. Since the total cross-section has the form a = As^° at 
large s, we get the mentioned asymptotical expression for the multiplicity. Of course this result 
is valid only in the extreme limit s — » oo. At large but finite s the existence of two different 
pomerons leads to some additional non-trivial s-dependence. 

It is instructive to compare our asymptotic results with those obtained in the BFKL fixed 
coupling model. In the latter case the inclusive cross section is of course badly behaved in the 
k± — > limit due to scale invariance. It is also very different in the high k± limit. At very large 
k± such that In k± ~ y/lns one finds 

/ d 2 a \ (k 2 lS ) A e -^ 2 kl/aHy) 

(zr^r) ~ a(y) [ -^^- — _ (2.69) 

\dydk\)BFKL fc ± -oo k\ ln/c^vOn^ 

where a 2 ~ (Ins — 4y 2 /lns). Thus for large k± the BFKL cross section goes down faster than 
any power. Also one obtains that (lnfcj_) ~ y/\n s so that both k± and k\ grow with s. 

To bring these predictions in better correspondence with the physical reality, as mentioned, 
various modifications of this orthodox BFKL approach have been introduced. In particular 
in pl| , [45| the fusion of gluons via "fan" diagrams was assumed to take place at high gluonic 
densities, which was considered as a way to partially restore the s-channel unitarity. Then, 
under some additional assumptions, a asymptotic behaviour similar to ( 2.67] ) was found. 



In our model such a behaviour naturally follows without introducing additional assumptions or 
imposing the unitarity restrictions. 



2.2.2 Numerical results. 



Taking the wave function evaluated numerically in [15] we have computed the cross section 



In Fig. 2/7. a we present d 2 a/dydk± for the process 7*7* (in units c = 1). We have chosen 
the projectile photon to have virtuality Q = 5GeV/c and the target one to have virtuality 
P = lGeV/c. The center of mass energy is = 540GeV. 

Of course, processes involving hadronic targets or /and projectiles are much more interesting 
from the practical point of view. However these require some non-perturbative input for the 
colour densities of the colliding hadrons. A possible way to introduce it is evidently to convert 
eq. ( p. 12 ) into an evolution equation in 1/x and take initial conditions for it from the exist- 



ing experimental data. Postponing this complicated procedure for future studies, we use here 
a simpler approach, taking for the hadron (proton) a Gaussian colour density with a radius 
corresponding to the observed electromagnetic one. Such an approximation, in all probability, 
somewhat underestimates the coupling of the hadron to the pomeron, since the coupling of the 
latter to constituent gluons is neglected. Nevertheless, we hope that it gives a reasonable es- 
timate for the inclusive cross-section. Using this approach we get the inclusive jet production 
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Figure 2.7: Inclusive jet production cross sections for the processes 7*7* at y/s = 
540GeV with Q 2 = 25(GeV/c) 2 and P 2 = l(GeV/c) 2 (a), for j*p at Vs = 
540GeV with Q 2 = 2h{GeV / c) 2 (b) and for pp at the same c.m.s. energy 
(c). 



cross-sections for the "f*p and pp scattering shown in Fig. 2.7.b and 2.7.c respectively. 



A common feature of jet production in all processes is that the cross-section reaches a max- 
imum at k± ~ 1 GeV '/c from which it monotonously goes down both for smaller and larger 
k ± . 

In Fig. 2.5 the cross-section dajdy integrated over k± in the interval 0.5 -j- 20 GeV/c is 



presented for the process pp. The limitations in the numerical calculations of the wave functions 
do not allow us to study higher values of k±, so that we cannot numerically reach the region 
where the asymptotical behaviour ( 2.67|) is strictly valid. 

In Fig. we show jet multiplicities as a function of s for the three studied processes 7*7*, 
7*p and pp. As one observes, their magnitude and behaviour prove to be quite similar. 

We also tried to estimate the average {k±). Unfortunately, although it exists according to 
( |2.67| ), its value results very sensitive to the high momentum tail of the pomeron wave function, 
poorly determined from our numerical calculations. To avoid this difficulty we chose to calculate 
the average (lnfci/A) instead. This average at fixed y depends on s and y only due to the 
existence of two different pomerons, so that this dependence should die out at large enough 
s. In fact the found average (lnfcj_/A) turns our to be practically independent of y and very 
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Figure 2.8: Cross-sections da/dy for the pp process at y^s = 540GeV. 
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Figure 2.9: Multiplicities (n) as a function of the center of mass energy i/s for the 
processes 7*7* (the solid curve), ~fp (the dashed curve) and pp (the dotted 
curve) . 



weakly dependent on s in the whole studied range of s and y, rising from 3.97 at \fs = 20 GeV 
to 4.19 at \/s = 20 TeV. These values imply a rather high average k± rising from 10.6 GeV/c 
at ^ = 20 GeV to 13.2 GeV/c at y/s = 20 Tel/. 

Our calculations show that the introduction of a running coupling constant on the basis of the 
bootstrap condition cures all the diseases of the orthodox BFKL approach for jet production. 
At high k± the cross-section becomes well-behaved and more or less in accordance with the 
expectations based on the quark counting rules. At small k± no singularity occurs, although 
contributions from other states is expected to dominate. 

As to s- and y -dependence, our predictions are even simpler than in the BFKL approach, 
since the running of the coupling converts the branch point in the complex angular momentum 
plane, corresponding to the BFKL pomeron, into poles, of which only two are located to the 
right of unity and contribute at high energies. As a result at super-high energies, when only 
the dominant pomeron survives, the y-dependence completely disappears and the s-dependence 
reduces to the standard s A ° factor. At smaller s some y- and extra s-dependence appears due 
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to the existence of two supercritical pomerons. 

With all these refinements, some basic predictions of the BFKL theory are reproduced. The 
cross-section for mini-jet production rises fast and saturates the total cross-section as s — > oo. 
Jet multiplicities rise logarithmically. 

However (k±) turns out to depend on s weakly and its calculated value results pretty high, 
of the order of 10-12 GeV/c. This should be contrasted to the experimentally observed much 
lower values of (k±) rising with energy. A natural explanation of this discrepancy follows from 
the results |l5|] shown in the first section of this chapter in particular from the study of the 
structure functions and total cross-sections in our model: at present energies the contribution 
from the two supercritical pomerons only covers a part of the observed phenomena because we 
are still rather far from the real asymptotics. The bulk of the contribution comes from other 
states, which produce a much softer spectrum of particles. The observed rise of the (k±) is then 
related to the dying out of these sub-asymptotical states. Our prediction is then that the rise 
of (k±) should saturate at the level of 10-12 GeV/c. 

This circumstance has also to be taken into account when discussing the numerical results 
presented in Figs. 6-9. Probably they also illustrate predictions for considerably higher energies 
than the present ones, at which, according to the estimate made in p5| , they should account for 
~ 20% of the observed spectra. 

To obtain predictions better suited for present energies one should evidently take into ac- 
count all states present in the spectrum of the pomeron equation (2.12) and not only the two 



supercritical pomerons. To realize this program an evolution equation in v following from (2.12) 
seems to be an appropriate tool. As mentioned, it could also effectively take into account the 
non-perturbative effects related to the pomeron coupling to physical hadrons. 



2.3 Evolution of the gluon density. 



For the forward scattering amplitude the pomeron equation ( 2.12| ) for the isotropic wave function 
with an inhomogeneous term reads 

(H - E)ip = Vo (2.70) 
Taking the Mellin transformation of (|2.7C| ) one converts it into an evolution equation in 1/x: 

—^—^(x, k) = -H${x, k) (2.71) 
a in 1/x 

which should be supplemented with an initial condition at some x = xq 

^(x ,k)=Mk) (2.72) 

containing the non-perturbative input about the coupling to the hadronic target. 

The physical interpretation of the pomeron wave function is provided by the fact that in 



the DLL A scheme eq. ( 2.71 ) reduces to an equation for the fully amputated function 4>(x,k) 
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r)(k)ip(x, k): 

d d , , , . 3a s (k 2 ) , . , . „„, 

ai, i -»i.i/, fet » = ^ ( ''" (2J3) 

which coincides with the standard equation for the unintegrated gluon density xg(x,k 2 ) in the 
DLL A limit. In fact, this circumstance lies at the root of our method to introduce a running 
coupling into the scheme. Thus we may identify 

4>(x, k) = cxg(x, k) (2.74) 



The normalizing factor c cannot be determined from the asymptotic equation ( 2.73| ). We shall 



be able to fix it by studying the coupling of the pomeron to the incoming virtual photon in the 
next section. 



2.3.1 Coupling to the virtual photon 

Once the function 4> proportional to the gluon density is determined, one has to couple it to the 
projectile particle to calculate observable quantities. In particular, to find the structure function 
of the target one has to couple the gluons to the incoming virtual photon, that is, to find the 
colour density p(q, k) which connects the photon of momentum q to the gluon of momentum 
k. This problem is trivial within the BFKL approach with a fixed small coupling. Then it is 
sufficient to take the colour density in the lowest order po(q, k), which corresponds to taking for 
it the contribution of a pure quark loop into which the incoming photon goes. 

The problem complicates enormously when one tries to introduce a running coupling into 
p. Then one has to take into account all additional gluon and qq pair emissions which supply 
powers of the logarithms of transverse momenta. Apart from making the coupling run, they will 
evidently change the form of p(q,k). Unfortunately the bootstrap relation can tell us nothing 
about the ultimate form of the colour density with a running coupling, which essentially belongs 
to the t-channel with a vacuum colour quantum number. So we have to find a different way to 
introduce a running coupling into p. 

A possible systematic way to do this consists in applying to the photon-gluon coupling the 
DGLAP evolution equation. One may separate the colour density from the rest of the amplitude 
by restricting its rapidity range to some maximal rapidity yo ~ log Q 2 (which, of course should 
be much smaller than the overall rapidity Y ~ logQ 2 /x). Then the kinematical region of p(q, k) 
will admit the standard DGLAP evolution in Q 2 . Solving this equation one will find the quark 
density at scale Q 2 of the gluon with momentum k (i.e. essentially the structure function of the 
gluon with the virtuality k 2 ). This is exactly the quantity needed to transform the calculated 
gluon density created by the target into the observable structure function of the target. As a 
starting point for the evolution one may take the perturbative colour density po at some low Q 2 
when the logs of the transverse momenta might be thought to be unimportant. 

This ambitious program, combining both evolution in 1/x and Q 2 , does not, however, look 
very simple to realize. As a first step, to clearly see the effects of the introduction of a running 
coupling preserving the gluon reggeization, we adopt a more phenomenological approach here, 
trying to guess a possible correct form for p(q, k) on the basis of simple physical reasoning and 
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Figure 2.10: Forward amplitude for the scattering of a photon with a reggeized gluon.. 



also using the DLLA to fix its final form. 

With a pure perturbative photon colour density one would obtain for the 7*p cross-section 

/ n 2N f d 2 kp (q,k)^(x,k) 
" (X ' 9)= y (2vr)V(fe) ( } 

In fact, the projectile particle should be coupled to the full pomeron wave function (fr/rj 2 . From 
the physical point of view this expression is fully satisfactory for physical particles. However it 
is not for a highly virtual projectile. 

To see this, we first note that for the forward amplitude our method of introducing a running 
coupling reduces to a very simple rule: the scale at which the coupling should be taken is given 



by the momentum of the emitted real gluon ((k — k') 2 in the upper rung in Fig. 2.10 ). Now take 
Q 2 very large and apply the DLLA approximation. Then the momenta in the ladder become 
ordered from top to bottom 

Q 2 » k 2 » k' 2 » 



In this configuration, as can be traced from (|2.70| ) and (2.75), all a s 's acquire the right scale 
(i.e. corresponding to the DGLAP equation) except for the upper rung: a s (k 2 ) appears twice. 



In order to clearly see why this happens let us explicitly write eq. ( 2.70| ) which describes the 



interactions in the ladder of Fig. p.lOj (the kinetic part given by the trajectories of the two 
reggeized gluons is denoted by T) 

(E - T {k))m _ *,(*) + 2 J^ n(k ,^_ kl) m ("6) 
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Constructing by iteration the solution, for example at the second order, one has 

d*)(k) = 2 f ^ V{h) 2 f 1 Mk " ] (2 77) 

y V ; E- T(k) J (2tt) 2 rj(k -k')E- T(k') J (2tt) 2 r](k")rj(k' - k") E - T(k") K ' ' 

and by taking the limit of the DLLA approximations one obtains 

{2)n ^_ 2 f d 2 k' 1 2 f d?k" 1 0oO" 



^ E - T(k) J (2tt) 2 rj(k') E - T(k') J (2tt) 2 V (k") E - T(k") ^' T8 " > 



It is evident that after taking this limit there is no cancellation of one f]{k) term in ( 2.75 ). 



Therefore this defect can be understood if one notices that the upper gluon is, in fact, coupled 



to a virtual particle. If this particle were a gluon, then the interaction K q in ( |2.1l| ) would cancel 
one of the two a(/c 2 )'s and substitute it by an a taken at the scale corresponding to its own 
virtuality. We assume that something similar should take place also for virtual quarks to which 
the gluon chain may couple. The scale of the particle momenta squared which enter the upper 
blob in Fig. 10 should have the order Q 2 (this is the only scale that remains after these momenta 
are integrated out). As a result the lowest order density should be rescaled according to 

Po(q,k)^^^ Po (q,k) (2.79) 

where Q\ has the same order as Q 2 . 



The approximation we assume is that the substitution ( 2.79| ) is sufficient to correctly repre- 



sent the photon colour density with a running coupling. We shall check its validity by studying 
the quark density which results from ( |2.79| ) in the DLLA approximation and comparing it with 
the known result based on the DGLAP equation. 

Explicitly the zeroth order density po has the following forms for the transverse (T) and 



longitudinal (L) photons (see e.g. [44] and do the integration in the quark loop momenta) 

pV\q, *0 = f^E Z f f Q da (V + (1 " «) 2 )((! + 2z2 )9(z) " 1) + Q(1 _ C a)+c (l " *(*))) 

(2.80) 

Here the summation goes over the quark flavours. The dimensionless variables £ and z are 
defined as 

rnj _ k 2 1 
C "g2"' Z ~ Wa(l-a)+( (2 ' 82) 
and m,f and Zf are the mass and charge of the quark of flavour /. The function g(z) is given by 

g{z)= / ln ^H + z (2. 83) 
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The structure function is obtained from the cross-section by the standard relation 



F2{x ,Q2 ) = 9l {a (T) +a (L) ) {2M) 



In the DLLA limit only the transverse cross-section contributes. We can also neglect the 
quark masses in this approximation. Then, with a substitution ( |2.79| ), from ( 2.75| ), ( |2.80| ) and 



(2.84) we obtain an expression for the quark (sea) density of the target 



where g(z) is given by Eq. ( 2.83| ) and we assumed that large values of k 2 < Q 2 contribute in 



accordance with the DLLA approximation. In this approximation the asymptotics of the gluon 
density xg(x, k 2 ) and consequently of (j)(x, k 2 ) is known: 



(x, k 2 ) = cxg(x, k 2 ) ~ cexp \j a\n — In In k 2 (2.86) 



where a = 48/6q. Putting ( 2.86| ) into ( |2.85| ), after simple calculations we find the asymptotical 



expression for the quark density 



, l2 , 4c /In In k 2 f~ 1 n , , " , „. 

xq(x, k J ~ , , i — -r- exp Wain — mm k (2.87) 

7r z 0q y am 1/x V x 

On the other hand, let us consider the DGLAP equation in the small x limit (in such a case it 
is possible to neglect the quark contributions to the sea quarks) 

d , oi s {Q 2 ) f 1 dy x 



x Q 2 ) = ^Rl[ dy Pqg ^ )g{yiQ 2 ) . Pqg{z) = 2 [z 2 + (1 - z) 2 } (2.88) 
ott j x y y 



dlnQ 2 ^ ^ ' 8tt J x y q9K y 
After a double integration we find 



, , n N 4 /lnlnfc 2 / 1~ ~ 

xq(x, k ) ~ —ro\ — \ t exp \/aln — lnlnfc^ (2.89) 

36q y am 1/x V x 



As we observe the approximation ( [2.79 ) for the colour density of the photon projectile leads to 



the correct relation between the quark and gluon densities in the DLLA limit. This justifies the 
use of ( |2.79| ), at least for high enough 1/x and Q 2 . Comparing ( p. 87 ) and (|2.89| ) we also obtain 



the normalization factor c which relates the pomeron wave function to the gluon density 

c = 7r 2 b /3 (2.90) 



2.3.2 The initial distribution 



To start the evolution in 1/x we have to fix the initial gluon density at some small value x = xq . 
Evidently, the smaller is xq, the smaller is the region where we can compare our predictions 



2.3 Evolution of the gluon density. 
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with the experimental data. On the other hand, if xq is not small enough, application of the 
asymptotic hard pomeron theory becomes questionable. Guided by these considerations we 
choose xo = 0.01 as our basic initial x although we also tried x = 0.001 to see the influence of 
possible sub-asymptotic effects. 

The initial wave function 4>{xq, k 2 ) has to be chosen in accordance with the existing data at 
x = xq and all k 2 available. The experimental F2 is a sum of the singlet and non-singlet parts, 
the latter giving a relatively small contribution at x = 0.01. Our theory can give predictions only 
for the singlet part (and one of the criteria for its applicability is precisely the relative smallness 
of the non-singlet contribution). The existing experimental data at x = 0.01 give values for F2 
averaged over rather large intervals of x and Q 2 . For all these reasons, rather than to try to 
adjust our initial <j>(xo, k 2 ) to the pure experimental data, we have preferred to match it with the 
theoretical predictions for the gluon density and the singlet part of F2 given by some standard 
parametrization fitted to the observed F2 in a wide interval of Q 2 and small x. As such we have 



taken the GRV LO parametrization [54|. The choice of LO has been dictated by its comparative 



simplicity and the fact that at x = 0.01 the difference between LO and NLO is insignificant. 
Thus, for the initial distribution we have taken the GRV LO gluon density with an ap- 



propriate scaling factor. Putting this density into eqs. ( 2.74 ),( 2?75|) and ( 2.84 ) one should be 



able to reproduce the sea quark density and thus the singlet part of the structure function. In 
the GRV scheme the relation between the gluon density and the quark density is much more 
complicated and realized through the DGLAP evolution. Since the DGLAP evolution and the 
pomeron theory are not identical, one should not expect that our initial gluon density should 
exactly coincide with the GRV one to give the same singlet structure function. One has also to 
have in mind the approximate character of our colour density p at small Q 2 . In fact, with the 
initial <p given by Q2.74| ) and the gluon density exactly taken from the GRV parametrization at 
x = 0.01 we obtain a 30% smaller values for the singlet part of the structure function as given 
by the same GRV parametrization, the difference growing at low Q 2 . To make the description 
better we used a certain arbitrariness in the scale Q 2 which enters ( 2.79Q and also the scale at 



which the coupling freezes in the density p. A good choice to fit the low Q data is to take 

2 \ 4-7T 1 

a ^ Ql > = Mn((0.17*Q 2 + 0.055 (GeV/c) 2 )/A 2 ) (2 ' 91) 

With this a(Q 2 ) the obtained singlet structure function at x = 0.01 has practically the same 
Q 2 dependence as the GRV one, although it results 30% smaller in magnitude. This mismatch 
can be interpreted in two different ways. Either we may believe that the gluon density given 
by the GRV is the correct one and the deficiency in the singlet part of the structure function 
is caused by our approximate form of the colour density p (which is most probable). Or we 
may think that the colour density to be used in the DGLAP should coincide with ours only for 
large enough Q 2 and 1/x and at finite values they may somewhat differ (our relation ( |2.74[ ) was 
established strictly speaking only in the DLLA limit). Correspondingly we may either take the 



relation (2/74) as it stands and use the GRV LO gluon density at x = 0.01 in it, or introduce 
a correcting scaling factor 1.3 which brings the structure function calculated with the help of 
fl2.75|) - (|2.84 ) into agreement with the GRV predictions. In the following we adopt the second 
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Figure 2.11: The singlet part of the structure function of the proton at x — 0.01. 

The continuous line is the result of our calculation while the dashed line 
correspond to the GRV prediction. 



alternative, that is we assume that our initial gluon distribution at x = 0.01 is 30% higher that 
the one given by the GRV parametrization. Of course to strictly compare the rate of growth 
of the distribution under evolution one should take into account the differences at the starting 
point. 

The singlet part of the structure function at x = 0.01 calculated from ( |2.75| )-( 2.84 ) with this 



choice is shown in Fig. |2.11 together with the GRV predictions. However one can easily pass to 
the first alternative by simply reducing our results by factor 1.3. 



2.3.3 Numerical results 

With the initial wave function (j)(x = 0.01, k 2 ) chosen as indicated in the preceding section we 
solved the evolution equation for 10~ 8 < x < 10~ 2 . 

The adopted computational scheme consists in diagonalizing the Hamiltonian in (|2.70|) , re- 
duced to one dimension in the transverse momentum space after angular averaging, and represent 
the initial wave function as a superposition of its eigenvectors. To discretize k 2 a grid was in- 
troduced, after which the problem is reduced to a standard matrix one. To check the validity of 
the obtained results we have also repeated the evolution using a Runge-Kutta method, resulting 
in a very good agreement. The final results obtained for the gluon distribution xg(x,Q 2 ) as a 
function of x for various Q 2 are shown in Fig. 2.12 where the first two plots correspond to x and 



Q presently available, whereas the last two plots show the behaviour of the calculated gluon 
density in the region up to very small x and very high Q 2 , well beyond the present possibilities. 
For comparison we have also shown the gluon densities for the GRV LO parametrization |54 



for the MRS parametrization [55] and also for the BFKL evolution as calculated in |5q| . 

Putting the found gluon densities into eqs. ( 2.75| )- (|2.84 ) we obtain the (singlet part of) proton 
structure function F2(x,Q 2 ). The results are illustrated in Fig. 2.13| for the x dependence. As 
for the gluon densities, the experimentally investigated region is shown in the first two plots, in 



which the existing experimental data from [ 10 1 are also presented. 

Finally, to see a possible influence of sub-asymptotic effects, we have repeated the procedure 
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taking as a starting point for the evolution a lower value x = 0.001. The resulting gluon 
distributions and structure functions are also presented in the above figures. 

To discuss the results obtained we have to remember that they involve two quantities of a 
different theoretical status. One is the pomeron wave function <p which can be identified with 
the gluon distribution (up to a factor). The other is the quark density (which is equivalent to 
the structure function), for which we actually have no rule for the introduction of a running 
constant and which in the present calculation involves a semi-phenomenological ansatz ( 2.79 ). 
Evidently the results for the latter are much less informative as to the effect of the running 
coupling introduced in our way. Therefore we have to separately discuss our prediction for the 
gluon distribution, on the one hand, and for the structure function, on the other. 

Let us begin with the gluon distribution. Comparing our results with those of GRV, which 
correspond to the standard DGLAP evolution, we observe that at high enough Q 2 and low 
enough x our distributions rise with Q and 1/x faster than those of GRV. This difference is, 
of course, to be expected. The hard pomeron theory in any version predicts a power rise of 
the distribution with 1/x to be compared with ( [2.86 ) for the DGLAP evolution. As to the Q- 
dependence, the fixed coupling (BFKL) hard pomeron model predicts a linear rise, again much 
stronger than ( 2.86 ). Our running coupling model supposedly leads to a somewhat weaker rise. 
From our results it follows that it is still much stronger than for the DGLAP evolution. However 
one can observe that these features of our evolution become clearly visible only at quite high Q 
and 1/x. For moderate Q < WGeV/c and/or x > 10 -4 the difference between our distributions 
and those of GRV is insignificant. As to the DGLAP evolved MRS parametrization, it gives the 
gluon distribution which lies systematically below the GRV one and, correspondingly, below our 
values, the difference growing with Q and 1/x. 



To compare our gluon evolution dynamics with the BFKL one, as presented in 5C], we have 
to remark that we should take a huge correcting factor that depends in Q 2 for the following 
reasons: we must remember of the factor 1.3 discussed in section 4 and we should note that 



the starting values used in [56] for the gluon distribution are lower than the GRV values, by a 
factor 1.92 at Q 2 = 4 GeV 2 /c 2 and by a factor 1.1 at Q 2 = 1000 GeV 2 /c 2 . According to this to 
compare just the evolution in this two models we should multiply our distribution by a factor 
that is 0.4 at Q 2 = 4 GeV 2 /c 2 and that is 0.7 at Q 2 = 1000 GeV 2 /c 2 . Taking into account this 
factors looking at Fig. 2.12 we see that there is a good agreement for not too small x values 
between the two evolutions, while looking at the values at very small x in all the Q 2 interval our 
model predicts a quite slower growth of the gluon distribution. We note also that the spectral 
method we have used shows how the dominant components of the spectrum contribute to the 
evolution of the gluon density; some small changes of the initial gluon distribution may results in 
different projections onto the spectral basis, but cannot lead to a very different behaviour at very 
small x, even if, as it is, the initial gluon distribution is well known in a restricted Q 2 interval. 
This is because the first two dominant states are precisely giving an important contribution in 
this interval. 



Passing to the structure functions we observe in Fig. |2,13| that our results give a somewhat 
too rapid growth with 1/x in the region 10~ 3 < x < 10~ 2 as compared to the experimental 
data (and also to the parametrizations GRV fitted to these data). With the scaling factor 1.3 
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Figure 2.12: The gluon distributions as a function of x evolved from x = 0.01 and 
x = 0.001 for the experimentally accessible kinematical region and for 
asymptotically high values of Q 2 and 1/x. Standard DGLAP evolved 
parametrizations (GRV-LO and MRS) and the BFKL evolved distribu- 
tions from Ji(J (we report only few points connected by lines) are shown 
for comparison. For the latter distribution one should take into account a 
further factor which actually has been chosen in ]56j to reduce the initial 
distribution (see text). 
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Figure 2.13: x dependence of the singlet part of the proton structure function obtained 
by evolution from x — 0.01 and x = 0.001, compared to the GRV prediction 
and the ZEUS 94 data. Values for asymptotically high values of Q 2 and 
1/x are also shown. 



introduced to fit the data at x = 0.01 we overshoot the data at x < 10~ 3 by ~25%. Without 
this factor we get a very good agreement for x < 10~ 3 but are below experiment at x = 0.01 
by the same order. This discrepancy may be attributed either to sub-asymptotic effects or to 
a poor quality of our ansatz ( 2.79[ ). Comparison with the result obtained with a lower starting 
point for the evolution x = 0.001 shows that sub-asymptotic effects together with a correct form 
of coupling to quarks may be the final answer. 



58 Hard pomeron with a running coupling constant: a phenomenological approach 



Chapter 3 



Unit arizat ion in the large N c limit: 
the three pomeron vertex 



In the first chapter we have reviewed and analyzed some aspect and properties of the BFKL 
approximation to the pomeron which is based on the application of the perturbation theory to 
QCD in the Regge kinematical region and corresponds to small in processes such as DIS. 
Specifically the BFKL pomeron is based on the resummation of the leading large logarithms in 
1/x and presents the important feature of a steep increase of the cross sections as x decreases, 
with a power like behavior in 1/x. This is, as we already noted in the first chapter, in contrast 



with the bound set by unitarity. In fact this is the well known Froissart bound [41], which 
sets the limit to the increase behavior with energy of the total hadronic cross section to In 2 s 
and which, for DIS, would lead to the growth bound of In 2 1/x [42, 43 1. Taking the double 
logarithmic limit, where one resums powers of (In 1/x In Q 2 ) and which thus come out also from 
the small x limit of the DGLAP approach, one has a somewhat less dramatic increase, in accord 
to an exp (^1/x) behavior, but still unitarity is strongly violated. 

Therefore, to extend the kinematical region of validity of the approach one has to go beyond 
the leading logarithmic approximation. This is in principle a formidable task and could also 
not be very satisfactory in presence of hadrons which cannot be studied only with perturbative 
methods. Here we shall be interested only in a perturbative analysis, which nevertheless is 
better established for the study of processes like onium-onium scattering, and one hopes that 
any result could give some new insight to the complete theory which surely would need of a 
non-perturbative treatment. 

A systematic approach to take into account a minimal set of corrections necessary to restore 



the unitarity has been studied by Bartels [19]. The higher order amplitudes, which contain any 
number of (reggeized) gluons in the i-channel, are constructed using unitarity and dispersion 
relations starting from the lower order ones and one has a coupled system of equations for the n 
reggeized gluon state amplitudes. The main feature is the presence of transitions vertices which 
change the number of reggeon in the t-channel and thus leads to a reggeon field theory. 

The study of this problem is far from being complete. The first step has been to consider 
the coupled equations for the amplitudes relative to 2, 3 and 4 gluons in the t-channel [pCj ]. The 
most interesting object which appears is an effective vertex for the transition between 2 and 4 
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gluon in the t-channel. This vertex has been found to be conformal invariant (21]] and in the 
appropriate color subspace is related to the triple pomeron vertex. Another element of some 
interest is the 4-gluon interacting state (and its generalizations to n gluons) satisfying the so 



called BKP equation p^, 23]. This amplitude, characterized by a constant number of gluons 



in the f-channel, has been studied in the large iV c limit and the system has been proved to be 
complete integrable and equivalent to an XXX Heisenberg zero spin model |24|, (2(| . 

The BFKL amplitudes and the 2 — > 4 effective reggeon vertex have been interpreted in terms 
of correlation functions of a 2D conformal field theory [27, [2£| hoping to obtain in the future 
such a formulation for the effective QCD in the Regge limit. In this context some work [57] has 
been done to understand the spectrum of the 4 gluon state in terms of the scaling dimension of 
some operators appearing in the twist expansion of the amplitude in the short distance limit. 
The idea is based on the fact that, due to the strong symmetry of a bidimensional conformal 
field theory, all the informations are contained in the scaling dimensions of the operators present 
in the OPE (operator product expansion) and in the expansion coefficients. 

We shall be interested here in the study of the coupled system of equations up to 4 reggeized 
gluons in the i-channel in the large N c limit p9| , |30t| . The analysis will be performed at the 
leading and next to leading order (NLO) in 1/N C since it is at the latter level that the double 
pomeron amplitude will appear. A link to some results coming from the dipole colour approach 
let arise the possibility of a general constructions of the n pomeron amplitudes in the large N c 
limit. In order to better understand this one should go beyond the 4 gluon system. Already 
the 6 gluon system in the large N c limit is capable to show if some equivalence with the A.H. 
Mueller approach to dipoles is maintained, since it seems that the latter leads to a system of fan 
diagrams with only a triple pomeron interactions. 

We shall start introducing the formal general system of equations for the 4 gluons system 
and its main properties. Written the system for the large iV c limit case the solutions in the 
leading order are given and at the NLO in 1/N C system studied. The double pomeron amplitude 
at NLO is found to separate in a direct double pomeron exchange (DPE) contribution and in 
a triple pomeron interaction (TPI) term. The TPI term can be written in terms of 3 point 
function building blocks which are shown to be conformally invariant so that the same property, 
as already known, has to be shared by the triple pomeron vertex. 

It is later shown in detail the interesting fact that the DPE contribution can be substituted 
by a completely equivalent TPI term (the inverse is not true). The absence of DPE terms was 
already noted for N c = 3 in |2(J and, in fact, in this form the final vertex interactions left 
essentially coincides in the two approaches, provided one takes the limit N c — > oo. Coupled 
to pomerons, this vertex effectively reduces to the one found by R. Peschanski in |H|], in a 
colour dipole framework. In particular this fact has been contemporaneously also shown by 
Korchemsky pq| . However this does not mean that the double dipole density in the dipole 
approach coincides with the diffractive amplitude in the s-channel unitarity approach: there are 
certain terms in the latter which are missing in the double dipole density and an asymmetry 
factor which distinguish, in the rapidity evolution, between the original pomeron and the two 
which originate from the split. 

Finally a section is devoted to study some properties of higher order dipole densities in the 
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A.H. Mueller approach which are found to be represented by a set of pomeron fan diagrams 
with only a triple pomeron coupling. 



3.1 General unitarization approach 

A program devoted to considering a minimal subset of next to leading logarithmic contributions 
necessary to construct a unitary amplitude has been carried out in J. Bartels work [pLS[ ] . The 
basic object under consideration is the amplitude D n for the coupling of n reggeons (gluons) 
to two external elementary particles in a colour singlet state. These amplitudes result from a 
multiple discontinuity analysis from which emerges the feature that the number of reggeized 
gluons in the t-channel is not conserved. This fact can be also analyzed in the physical i-channel 
region (t > 0) looking for the discontinuities along the i-cut associated with the threshold for 
the reggeon production. One nice feature due to the gluon reggeization property is the absence 
of contributions, in the adjoint colour representation, to the partial wave amplitude cut related 
to n-particle production with n > 2. 

Let us write the formal system of equations in cascade for the D n amplitudes (seen as vectors 
in the appropriate colour spaces) up to n = 4: 



3 ~ 1 - J2 U ( k i 

1=1 

3 

j - 1 - ^ufc 

i=l 
4 

j - 1 -J^^ih 
i=i 



D 2 {{h}) 

D 4 ({h}) 



D 20 ({h}) + [K 2 ^ 2 
D 30 ({ki}) + [K 2 ^ 3 



D 2 ]({h}) 

D2}({k}) + J2lK2^2®D 3 ]({k i }) 



D m {{ki}) + [K 2 ^ 4 ® D 2 ]({k l }) + J2[K 2 - 

5}#2-+2®£>4]({fci}) 



f3 ®D 3 ]({k i }) + 

(3.1) 



In the above equations D n Q is the lowest order term with the gluons coupled directly to the qq 
loop. K 2 ^ 2 is the BFKL interaction kernel with the gluon propagators included. Remembering 
the expression for its momentum dependence 



V(k 1 ,k 2 ;k' 1 ,k 2 



h^h' 2 + k 2 k' 2 
12 ' "'2 1 



+ k 2 f 



/^l (J^\ ^i ) ^ ^i ^ 



I 2 rj 2 



(3.2) 



it is possible to give the expression for the kernels K 2 ^ 3 and K 2 ^ 4 in a simple way. In fact, 
defining 

W(h,k 2 , k 3 ; k[,k' 3 ) = V(k 2 , k 3 ; k[ - h, k' 3 ) - V(h + k 2 , k 3 ; k[,k' 3 ) (3.3) 
one has, on including the colour structure, 



K 2 -, 3 (h, k 2 , k 3 ; k'x,k 3 ) 
K 2 ^ 4 (h, k 2 , k 3 ,k 4 ; k[,k' 4 ) 



g *f<^f^f^W{k u k 2 , k 3 ; k[, k' 3 ) 



Moreover in the coupled system of equations (|3.1| ) the sums run over all possible pairwise interac- 
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tions and the symbol (g) means an integration over primed variables with a weight (2ir) 3 5^ (J2 k— 

The analysis of the colour structure for the system of 4 gluons was performed by Bartels in 
|Lst l, separating the gluons in two pairs, each in the same colour quantum state belonging to the 
set ( 1, 8,4, 85, 10 + 10, 27). Therefore the amplitude has been characterized by a 5-dimensional 
colour vector for the three different signature assignments (— , — ), (+, +), (+, — ) and it has been 
possible to reduce the general operator K^4, which should appear in principle in the last term 
of the equation for D4, in a sum of operators acting on different colour subspaces. 



3.2 Large N c limit for the system of four reggeized gluons 



The set of equations (3.1) for the system of four reggeized gluons, coupled to a fermion loop in a 



colour singlet state, found by J. Bartels and M. Wuesthoff |2C] for the real case iV c = 3 is quite 
complicated because of the colour structure. 

The large N c limit jl], |2[ is interesting because in the leading order the amplitude for the 
system is based on a single BFKL pomeron exchange, i.e. one has a fully reggeizing contribu- 
tion. Performing a perturbative analysis at the NLO in 1/N C one finds a subleading diffractive 
amplitude which consists of a sum of a DPE and a TPI terms, in full correspondence with the 
Regge-Gribov picture. 



The results presented in this section where studied in [29] and |3(| respectively for the 
forward and for the more general non-forward direction case which is necessary to study also 
the conformal invariance properties of the 2 — * 4 effective gluon vertex. 



3.2.1 Leading order in 1/N C 

The results for a system of two and three gluons are not sensitive to the N c —* 00 limit, which 
instead affects the four gluon case. 

Let us first consider the two gluon system. In the lowest (zero) order approximation the 
2 gluon amplitude discontinuity (Mellin transformed to the complex angular momentum j) is 
given by the qq loop with the two gluons attached to it in all possible ways (4 diagrams in all). 
Their colour indices a\ and ai enter into the colour trace Tr{t ai t a2 } = (l/2)<5 aia2 , where t a is 
the colour of the quark. Separating this trace and the coupling g 2 and projecting on the vacuum 
colour state |0) = (l/N)6 

£tia 2 we write the zero order contribution in the large iV c limit as 

D 2 o(l, 2) = D 20 (2, 1) = g 2 N c (/(l + 2, 0) - /(l, 2)) (3.5) 

where we use the notation in which only the number of the gluon is indicated whose momentum 
enters as a variable and /(l, 2) = /(2, 1) is a contribution of the qq loop with gluon 1 attached 
to q and gluon 2 attached to q. Its explicit form can be easily found (see Appendix A. 2) but 
has no importance for the following. 

The basic quantity is the full amplitude D2. It corresponds to the exchange of two reggeized 
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gluons (the BFKL amplitude) and satisfies the BFKL equation 

S 20 D 2 = D 20 + g 2 N c V l2 D 2 (3.6) 
where S 2 o is the 2 gluon "free" Schrddinger operator for the energy 1 — j 

S 20 =j-1- cj(1) - w(2) (3.7) 

uo{k) is, as usual, the gluon Regge trajectory, V\ 2 is the BFKL interaction and the previous found 
D 2 q acts as an inhomogeneous term. For simplicity we leave here the (g> symbol understood. 

For the number of exchanged gluons n = 3,4 amplitudes D n (j) are defined as integrals of 
D n (ji, ■■■jn-i), depending on n — 1 partial t-channel angular momenta, over all ji subject to 
condition Yn=i(ji — 1) = J — 1 corresponding to conservation of "reggeon energies". Moreover 
we extend the notation for the n reggeized gluon free Schrddinger operator defining S n o = 

The three gluon system is described by an equation which can be solved in terms of the D 2 
amplitudes thanks to the gluon reggeization property. It is worth to study its colour structure. 
The colour factor associated to the lowest order amplitude of three gluons attached to the quark 
loop is given by the trace Tr{t ai t a2 t a3 } = (l/4)/i 0l0203 , where h aia2a3 = d aia2as + if ai a 2 a 3 is 
cyclic symmetric and satisfies the properties 

h *a ia2 a 3 = Kiaiaa , ^ Kcdhcd = S ab 2N c (l - 2/iV c 2 ) , ^ Kcdhbcd = -5 ab (A/N c ) (3.8) 

cd cd 

One may introduce the two colour wave functions 



1 123) — ^——^h aia2a3 , 1 213) — ^——^h a2aia . A (3.9) 



orthonormal in the large iV c limit and project the above lowest order amplitude on to these 
states. The result can be expressed in terms of the function / because the momenta of the 
gluons attached to the same fermion line sum together; there is also an overall sign depending 
on the parity of the number of gluons attached to the quark or to the antiquark lines. One has, 
for N c — > oo, 



L>£ 23) = -Dg 13) = gJ^ \D 20 (2, 1 + 3) - £> 20 (1, 2 + 3) - L> 20 (3, 1 + 2)1 (3.10) 




Due to the global colour singlet state, each pair is in an adjoint representation and the interaction 
cannot change the colour configuration. Remembering the bootstrap property each two reggeon 
state in an adjoint representation gives one reggeon state so that the considered three gluon 
amplitude must reduce to a superposition of BFKL pomeron amplitudes. 

In more detail one can write the equation for the colour state 1 123) 

S 30 D^ 3) = < 3) + + ^Vi2 + V23 + V 31 )D 3 (3.11) 
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( 123) 

The second inhomogeneous term -D2— >3 ^ s related to the contribution which arises for the tran- 
sition from 2 to 3 gluons. In fact, using the vertex -^2^3 recalled in the previous section and 
projecting onto the state 1 123) one finds 



= 9 3 \l ^W(l, 2, 3; 1'3') D 2 {\' , 3') = g\\ ^W 2 (l, 2. 3) 



(3.12) 



In the equation for the orthogonal state 1 213) the inhomogeneous terms have opposite sign. As 
expected the solutions of the equations can simply be obtained promoting the D 2 q in the zero 
order term to a full BFKL pomeron amplitudes D 2 , i.e. 



D 



(123) 



-D 



(213) 




D 2 (2, 1 + 3) - D 2 (l, 2 + 3) - D 2 (3, 1 + 2) 



(3.13) 



This can be easily proved (see appendix A. 4) directly using ( |3.3j ) to transform D 2 ^ an< ^ °y 
applying the bootstrap condition 



„2/U- r j2jL/ 2 AT 

" " (^^i' fc 2;^,A; 2 )^^y 12 ®l = u;(l)+ W (2)- W (l + 2) 



(3.14) 



to obtain contributions coming from BFKL equations for amplitudes where two reggeized gluons 
have fused into one. 



Let us now consider the four gluon system. In the leading order the colour configuration is 
much more conveniently described by the order of gluons along the cylinder surface 1234, 1324, 
etc. If one takes the lowest order amplitude given by the contributions of 16 diagrams with the 
4 gluons attached in all possible ways to the qq loop, one has for each case a color factor given 
by the trace Tr{t ai t a2 t a3 t ai } = (l/8)h aia2b h a3a4b + (l/4N c )5 ai a 2 5 a3a4 . In the N c -> 00 limit only 
the first term, cyclic symmetric, surveys. It is useful to define the colour state 



1 1234) 



2iV c 2 



h. 



aia2bha3a,4b 



(3.15) 



together with others, differing for a permutation, orthonormalized in the large N c limit. The 
zero order amplitude projected on to these states gives the components 



D 



(1234) 
40 



D 



(4321) _ g 2 N c - 



10 



D 20 (l, 2 + 3 + 4) + Z? 20 (4, 1 + 2 + 3) - D 20 (l + 4, 2 + 3) (3.16) 



n (2134) _ n (4312) _ g N c [■ 
-^40 — -^40 — 



Z?2o(2,l + 3+4)+Z? 2 o(3,l + 2+4)-D2o(l + 2,3+4)-D 2 o(l + 3,2+4) 

(3.17) 

The gluons lie on a cylinder attached to the qq loop and all the neighbouring pairs are globally 
in an adjoint representation. Moreover in the leading order in 1/N C all the interactions across 
the cylinder are suppressed and only the neighbouring gluons interact. 
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The equations for the above amplitudes decouple and are given by [^9| 

S m D?^ = D^ + D^ + ^(V 21+ V 13 + V M + V i2 )Dt U) (3.18) 



The inhomogeneous terms related to the transitions which change the number of gluons from 2 
to 4 are 



,(1234) ff 4 iV c 2 



D 2->7' = -Z-j JL W 2 (l,2 + 3,4) (3.19) 



while the transitions between 3 and 4 gluons are due to the terms 



DiZT =9 3 \r~f [W(2, 3, 4; 2', 4') ® D^> (1, 2', 4') + W(l, 2, 3; 1', 3') ® Dg™> (1', 3', 4) 

= -9 3 \[^-[W(h 2, 4; 1', 4') ® Z)f 4) (l', 3, 4') + 3, 4; 1>, 4') ® Z)f 4) (!', 2, 4') 

(3.20) 

It can be shown that also the solutions for the 4 gluon amplitudes can be written in terms of 
the two reggeon amplitudes just by substituting the -D20 terms in ( |3.16|) and ( 3. 171 ) by D 2 ones. 




' a'\x* n( 124 )n <V a'\ 1 wn o q. 1' q'\ ^ n( 134 )n' q' 



Thus in the leading order in 1/N C the 4 reggeized gluon system reduces to a single pomeron. 
3.2.2 Next to leading order in 1/N C 

The corrections at the NLO in 1/N C force to consider the two subsystems of gluons (12) and (34). 
The state which relates to the diffractive amplitude and contributes to the triple pomeron inter- 
action is characterized by both the subsystems in a colour singlet state so that we shall project 
the zero order state of the 4 gluons attached to the quark loop onto |0) = (l/A r2 )J aia2 (5 a3a4 , 
obtaining 

1 4 4 
D fo = ^5 2 (E D ^(h 1 + 2 + 3 + 4- i) -J2 D 20 (l + i, 2 + 3 + 4 - »)) (3.21) 

i=l i=2 

which is explicitly down by a factor 1/N C respect to the 4 gluons amplitudes at leading order 
(its order is g A N c compared to g A N^ of the amplitudes D^ 3 ^ and Z)^ 2 , 134 ^). The full diffractive 
amplitude (directly related to the diffractive cross-section integrated over the diffractive 
mass) results to be the solution of the following equation 



S40 Df = + DH 4 + DSU + D®U + g 2 N c (V 12 + V U )D^ (3.22) 

where there is also an inhomogeneous term which comes from a LO 4 gluon cylindric 

configuration due to the interactions between adjacent gluons as well as across the cylinder. In 



the Bartels approach [19| the corresponding transitions are related to the matrix elements 12 
and 13 (which turns out to be zero) of the matrix Kj_^ 4 . The contributions to the two colour 
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singlet amplitude at the NLO have explicitly been found [p9f : 



D 
D 

D 



(0) 
2^4 

(0) 
3^4 



-<? 4 iW 2 (l,2 + 3,4) 



g 3 y^N c [W(l, 2, 3; l'3 ; ) ® D^ U) (1', 3', 4) - W(l, 2, 4; 1', 4') ® (1', 3, 4') + 

W(2, 3, 4; 2', 4') ® £>< 12 4) (1, 2', 4) - 3, 4; l', 4') ® 4) (l', 2, 4')' 

,(1234) n (2134) 



(134) n / 



(0) 
4^4 



g 2 (V 23 + V 1A -V l3 -V, 



21; 



(3.23) 



From the inspection of these inhomogeneous terms one can observe that each one is given by 
the action of some operator on the BFKL pomeron amplitude D 2 . It is then possible to define 
an operator Z associated to the triple pomeron vertex and acting on the BFKL pomeron 



D 



(o) 

2^4 



(0) 
4^4 



and rewrite eq. ( |3.22 ) as 



5*40^4 



(o) 



D { $ +ZD 2 + g z N c (V 12 + Vu)D% 



(0) 



(3.24) 



(3.25) 



The three-pomeron vertex operator Z can be conveniently reexpressed in terms of a new 
function G(k\, k 2 , k 3 ) (see also appendix A. 3) which is defined as the vertex K 2 -+ 3 for the tran- 
sition of 2 to 3 gluons, integrated with the BFKL pomeron and regularized in the infrared by 
terms proportional to the gluon trajectory in the same manner as in the total BFKL kernel 

G(h,k 2 , k 3 ) = G(k 3 , k 2 , h) = -g 2 N c W 2 (h, k 2 , k 3 ) - D(h, k 2 + k 3 )(uj{k 2 ) - u(k 2 + k 3 )) - 
D(h + k 2 , k 3 )(uj{k 2 ) - u{kx + k 2 )) (3.26) 

It is a generalization of a similar function of two momenta introduced in ^] for the forward 
case and in the graphical notation introduced by Bartels it can be written as illustrated in Fig. 

o. 

If we write ^3^4 using G we obtain (as usual we denote ki by writing just i for brevity), 
■2 



D 



(o) 
3^4 



sr_ 
2 



2 G(2, 3 + 4, 1) + 2 G(3, 1 + 2,4)- G(l + 3, 2, 4) - G(l + 4, 2, 3) 



(7(1 + 4,3,2) -G(2 + 4,3,l 


+ 




D 2 (l,2 + 3 + 4) (2o/(3 + 4) - 


-o/(3) 


-w(2 + 3 + 4)) + 


D 2 (2,l + 3 + 4) (2o/(3 + 4)- 


-o/(3) 


-o/(l + 3 + 4)) + 


L> 2 (3,l + 2 + 4) (2o/(l + 2)- 


-o/(2) 


-o/(l + 2 + 4)) + 


L> 2 (4,l + 2 + 3) (2o/(l + 2)- 


-o/(2) 


-o/(l + 2 + 3)) + 


L> 2 (l + 3,2 + 4) (o/(2 + 4) - 


o/(2)- 


-o/(3)) + 


D 2 (l + 4,2 + 3) (o/(2 + 3) - 


u/(2)- 


- W (3))" 



(3.27) 



(0) 



This expression is not symmetric in gluons (12) nor in gluons (34) and, since the function D\ 
is symmetric in both gluon pairs, one should symmetrize it. The resulting symmetric expression 
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Figure 3.1: Graphical representation of the conformal invariant function G in term of 
which the vertex Z is constructed: for each line one has to write the corre- 
sponding propagator, for each vertex the square of the sum of the momenta 
above or below that vertex and the propagators of the lower lines are omit- 
ted. 



is 



D 



(o) 

3^4 



2 r 



2 G(l, 3 + 4, 2) + 2 G(2, 3 + 4, 1) + 2 G(3, 1 + 2, 4) + 2 G(4, 1 + 2,3) - 
G(l + 3,2,4) -G(l + 3,4,2) - G(2 + 3,1,4) - G(2 + 3, 4,1) - 
G(l + 4,2,3) -G(l + 4,3,2) - G(2 + 4,1,3) - G(2 + 4, 3,1) + 
D 2 (l, 2 + 3 + 4) (2w(2 + 3 + 4)-2 w(2 + 4)-2w(2 + 3)+ w(3) + w(4)) + 
D 2 (2, 1 + 3 + 4) (2 w(l + 3 + 4) - 2 w(l + 4) - 2 w(l + 3) + w(3) + w(4)) + 
D 2 (3, 1 + 2 + 4) (2 w(l + 2 + 4) - 2 w(2 + 4) - 2 w(l + 4) + w(l) + w(2)) + 
£> 2 (4, 1 + 2 + 3) (2w(l + 2 + 3)-2 w(l + 3)-2w(2 + 3)+ w(l) + w(2)) 
£> 2 (1 + 3, 2 + 4) (2 w(l + 3) + 2 w(2 + 4) - w(l) - w(2) - w(3) - w(4)) + 
D 2 (l + 4, 2 + 3) (2 w(l + 4) + 2 w(2 + 3) - w(l) - w(2) - w(3) - w(4))l (3.28) 



The term D2-I4 can be rewritten as 



D. 



(0) 



2^4 



9 



G(l,2+3,4)+I>2(l,2+3+4)(t4;(2+3)-w(2+3+4))+D2(4,l+2+3)(w(2+3)-w(l+2+3)) 

(3.29) 

and, after an analogous symmetrization procedure, it becomes 



D 



{ 2 °i 4 = (l/4) 5 2 [G(l,2 + 3,4) + G(l,2 + 4,3) + G(2,l + 3,4) + G(2,l + 4,3) + 
D 2 (l, 2 + 3 + 4)(w(2 + 3) + w(2 + 4) - 2 w(2 + 3 + 4)) + 
^(2, 1 + 3 + 4)(w(l + 3)+w(l + 4)-2w(l + 3 + 4)) + 
£> 2 (3, 1 + 2 + 4)(w(l +4)+w(2 + 4)-2w(l + 2 + 4)) + 



68 



Unitarization in the large N c limit: the three pomeron vertex 



£> 2 (4, 1 + 2 + 3)(w(l + 3) + w(2 + 3) - 2 u(l + 2 + 3))J (3.30) 

We finally come to the term D^ 4 . The operator V 2 % + V14 — V13 — V24 is antisymmetric under 
the interchange of 1 and 2 and/or 3 and 4. Since the two-pomeron state is symmetric under 
these substitutions, only antisymmetric parts of _d| 1234 ^ and D^ 213 ^ give a nonzero contribution. 
Taking this into account we find 

^ TV 

£>E 4 = ^(^23 + ^14 - lis - V 24 )[D 2 (1 + 3, 2 + 4) - D 2 (l + 4, 2 + 3)] (3.31) 

As we have previously noted these terms are due to the interaction between different gluons (e.g 
Vi3D 2 (l + 4)) or within the same gluon (e.g. V\sD 2 (l + 3)). 



The first type gives the contribution 



D 



(01) _ g A N c 



4^4 



[(V13 + V 2i )D 2 {l + 4) + (V 23 + V U )D 2 (1 + 3)] 



(3.32) 



To reduce it to the expression of the same type as before we use the identity (see appendix A. 4): 
V 13 D 2 (1 + 4) = W 2 (2 + 3, 1, 4) + W 2 {2, 3, 1 + 4) - W 2 (2, 1 + 3, 4) - W 2 (2 + 3, 0, 1 + 4) (3.33) 
Introducing then functions G we find for this part 



D 



(01) _ g 4 N c 



4^4 



[G(l, 2 + 3, 4) + G(l, 2 + 4, 3) + G(2, 1 + 3, 4) + G(2, 1 + 4, 3) - 

G(l, 3,2 + 4)- G(l, 4, 2 + 3) - G(2, 3,1 + 4)- G(2, 4, 1 + 3) - 

G(3, 1,2 + 4)- G(3, 2, 1 + 4) - G(4, 1,2 + 3)- G(4, 2, 1 + 3) + 

2 G(l + 3, 0, 2 + 4) + 2 G(l + 4, 0, 2 + 3) + 

D 2 (l, 2 + 3 + 4)(w(2 + 3) + kj(2 + 4) - w(3) - w(4)) + 

D 2 (2, 1 + 3 + 4)(u;(l + 3) + u(l + 4) - w(3) - w(4)) + 

L> 2 (3, 1 + 2 + 4)(w(2 + 4) + w(l + 4) - w(l) - w(2)) + 

D 2 (4, 1 + 2 + 3)(w(l + 3) + kj(2 + 3) - w(3) - w(4)) - 

4 

(D 2 (l + 3, 2 + 4) + D 2 (l + 4, 2 + 3)) J2 w(i 



(3.34) 



The second part leads to the contribution 

„4 



£>K = "^[(^13 + ^24)^2(1 + 3, 2 + 4)+ 



(V 23 + V 14 )D 2 (l + 4,2 + 3)] 



(3.35) 
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Applying the bootstrap relation one immediately has 



D 



(02) 
4^4 



9 2 N C 



D 2 {1 + 3, 2 + 4) £ - w(l + 4)) - w(2 + 3) + 



i=l 



D 2 (l + 3, 2 + 4) (J] w(») - w(l + 3) - w(2 + 4))' 



i=l 



(3.36) 



The final three-pomeron vertex is obtained after summing all the contributions, explicitly ZD 2 
dRa + D^L + £>4°-S + and is given by: 



ZD 2 = -y[G(l + 3,2,4)+G(l + 3,4,2) + G(l,3,2 + 4) + G(3,l,2 + 4) + 
G(l + 4, 2, 3) + G(l + 4, 3, 2) + G(l, 4, 2 + 3) + G(4, 1,2 + 3)- 
2G(1, 3 + 4, 2) - 2G(3, 1 + 2,4)- G(l, 2 + 3,4)- G(2, 1 + 4, 3) - 
G(l, 2 + 4, 3) - G(2, 1 + 3,4)- G(l + 3, 0, 2 + 4) - G(l + 4, 0, 2 + 3)1 (3.37) 

One may note from the form of the K 2 ^z vertex that, after multiplying by the gluon propagators, 
the G function vanishes if the first or the third argument is zero, but the property in not true 
for the second argument. However at the global vertex level, after having performed the sum 
( |3.37 ) of all the contributions, the validity of the Ward identities is restored. 



Let us end this section looking at the form of the solution of the equation ( 3.25 ). In fact 



it can be easily solved and, evidently, the solution may be constructed as a sum of two terms 

inhomogeneous term 

D (0) =D DPE + D TPI (3>38) 



corresponding to the two parts of the inhomogeneous term D^j and ZD 2 : 



The DPE term D^ PE comes from the inhomogeneous term D^q and its explicit form can be 
conveniently written using the quark loop density in the transverse coordinate space defined by 
the Fourier transform (see Appendix A. 2 for explicit expressions) 

/(1,2)= J d 2 r Pl (r)e^ r (3.39) 

where I = k± + k 2 and r is the dipole transverse dimension. Thus one has a contribution of the 
form 

D° PE = (l/4)g*N c J d 2 r Pl {r)D { : ] (3.40) 

(r) 

Here D\ is a convolution in the "energy" 1 — j of two independent BFKL pomerons 

D { { ] = J dj 12 dj u 5(j + 1 - j 12 - ju)D ( ;] l2 (1, 2)£>g 34 (3, 4) (3.41) 

and to restore the structure of the inhomogeneous term the amplitudes must belong to two 
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equations with a proper inhomogeneous part. jD^(1,2) satisfies the equation 

SmDQ = f[ (e^ r - 1) + 9 2 N c V l2 D { Q (3.42) 
i=i 

and a similar relation regards the second pomeron. In the rapidity space one can directly write 
this contribution as a product of two BFKL amplitudes relative to the two colour singlet pairs 

Dt\y) = D^ l2 \y) D^\y) (3.43) 

The part Dj PI is the TPI contribution. It can be written as a convolution in the rapidity space: 

Dj PI (1,2,3 A;Y) = 

[ Y G 2 (l, 2; 1'2'; Y - y)G 2 (3, 4, ; 3'4'; F - y) ® 2', 3', 4'; l", 2") ® D 2 (l" , 2"; y) (3.44) 
JO 

where G 2 is the BFKL Green function and the symbols <8> mean integrations over intermediate 
momenta. This equation clearly shows that Z is just the three-pomeron vertex. 



3.3 Conformal invariance 

In [^] it was introduced the vertex ^(1234). For a non-zero momentum flow in the i-channel 
it can be written in a way similar to ( |3.37| ), 

V{12U)D 2 = ^g 2 [G(l, 2 + 3, 4) + G(2, 1 + 3,4)+ G(l, 2 + 4, 3) + G(2, 1 + 4, 3) 
- G(l + 2, 3, 4) - G(l + 2, 4, 3) - G(l, 2, 3 + 4) - G(2, 1,3 + 4) + G(l + 2, 0, 3 + 4)1 (3.45) 



This vertex is conformal invariant |2M in the following sense. If one transforms VD 2 to the 
transverse coordinate space and integrates it over the 4 gluon coordinates with a conformal 
invariant function, the resulting integral is invariant under conformal transformation of gluon 
coordinates. Comparing ( 3.45| ) and ( p. 37 ) one immediately see that our vertex Z can be written 



as 

Z = V(1324) + V(1423) (3.46) 

and thus is conformal invariant. 

However a stronger conformal property is satisfied: not only the combination ( 3.45| ) of func- 



tions G is conformal invariant, but each function G(l, 2, 3) is conformal invariant by itself. This 
is related to the fact that the function G defined in ( 3.26| ) represents a natural generalization of 



the BFKL kernel not only in respect to its infrared stability but also in its conformal properties. 

To prove the previous statement one needs a representation of G in the coordinate space. 
We give all the technical details in appendix A. 3 and use here the computed final expressions. 



Denote the integral part of G, given by the W 2 term in ( 3.26 ) as G\ and the remaining terms 
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related to the gluon trajectories as G 2 . In A. 3 the following expression is derived: 

G 1 (r 1 ,r 2 ,r 3 ) = A 1 D 2 (r 1 ,r 3 ) 
where Ai is an operator in the coordinate space 



(3.47) 



9 2 N cl 



8^ 3 



2Tr8 2 (r 23 )dl(c - lnr 13 )<9 3 - 2 + 2v5 2 {r l2 )d 2 (c - lnria)^ 



-2 



F12F23 

2 Jl 



2vr(c - lnr 13 )(<5> 12 ) + 5 z {r 23 )) - 4vr 2 J 2 (r 12 )5 2 (r 23 )(d 1 + d 3 fd^d 3 ' z (3.48) 



r 19 r. 



12' 23 



Here = r « — rj and c = ln(2/m) + ip{X) with m the gluon mass acting as an infrared cutoff. 

The transformation of the part G 2 requires introduction of an ultraviolet cutoff e due to the 
presence of the gluon trajectory terms. Of course the final results do not depend on e. One 
obtains 

G 2 (r 1 ,r 2 ,r 3 )=A 2 D 2 {r 1 ,r 3 ) (3.49) 



where A 2 is another operator in the coordinate space 



Ao 



9 2 N C 
8vr 3 



2vrc,5 2 (r 2 3) + 5 2 {r 23 )u:(-id 3 ) 



23 



9 2 N C [ 1 
8vr 3 



12 



27rc5 2 (r l2 )\ +5 2 (r l2 )oj(-id 1 ) 

(3.50) 

The four terms in A 2 correspond to the four respective terms in G 2 . We remember that the 
singular operators 1/V 2 2 and 1/V 23 are in fact defined with the help of e as 



1 



1 



r 2 + e 2 



+ 27r<5 2 (r)lne, e^O 



(3.51) 



which is the expression in (A.69) and does not depend on e. Summing A\ and A 2 the terms 
containing mm cancel, thus the dependence on the gluon mass desappears and G(ri,r 2 ,r 3 ) is 
infrared stable. 



We want to check the conformal invariance of the integral 

1=1 d 2 r l d 2 r 2 d 2 r 3 ^(r l ,r2,r 3 )G{r 1 ,r 2 ,r 3 ) (3.52) 

where function is conformal invariant. We shall demonstrate that the integral I does not 
change under conformal transformations. In doing so we shall use the fact that the BFKL 
solution ^{r\,r 2 ) = 9 1 ~ 2 5 2 ~ 2 L' 2 (ri, r 2 ) is conformal invariant, as it has been shown in the first 
chapter. 

It is sufficient to study the behaviour of the function G only under the inversion since the 
invariance under translations and rescaling is obvious. The complex notation (see section 3 of 
the first chapter) is more convenient for this aim: under inversion one has 



r — > 1/r, k = —id — > r 2 k 



(3.53) 
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Therefore (for real r) 



and 



d A r — > <i 2 r/r 4 



D 2 (ri,r 2 ) -» r 1 r 2 D 2 (r 1 ,r 2 ) 



(3.54) 



(3.55) 



There are some contributions which are evidently invariant under inversion. For example taking 
the last term from A\ it is easy to check that it leads to the integral 



9 2 N C 
2vr 



d 2 r<fr(r, r,r)d 2 ^(r,r) 



(3.56) 



which is conformal invariant, since both functions $ and \E r are invariant, and the factor r 4 
from d 2 r is cancelled by the factor r 4 from d 2 . 

Terms with the denominators r\ 2 and/or r| 3 from A% + ^2 combine into an integral 

g 2 N c f - 2 
I d 2 rid 2 r 2 d 2 r 3 ^(r 1 ,r 2 ,r 3 )- 

" 12' 23 



a 2 N f r 2 
h = / d 2 r l d 2 r 2 d 2 r^{r l ,r 2 ,n)^ r D 2 {r l ,r z ) (3.57) 

07T J ^19^9*3 



in which the regularization ( |3.51| ) is implied. Under inversion the ultraviolet cutoff e is trans- 
formed into ei = r\r 2 e and £2 = r 2 r%e in the denominators r\ 2 and r| 3 respectively. Consequently 
one finds a change of I 2 : 



AI 2 



9 2 N C 
Att 2 



d rid r 3 <E>(ri, ri, r 3 ) lnr^ + $(ri, r 3 , r 3 ) lnr 3 D 2 (n, ^3 



(3.58) 



All the other terms in Ai+Ajj proportional either to 5 2 (ri 2 ) or to 5 2 (r 2 3) can be divided into 
two parts. On one side one has terms in which the 5-function is multiplied either by a constant 
or by lnri 3 . It will give rise to a part of the integral I3. Terms with a constant are evidently 
invariant under inversion. However those containing lnri 3 are not and lead to the corresponding 
change of J 3 



A/ 3 



9 2 N C 
8vr 2 



d 2 r!d 2 r 3 *(n, r 3 , r 3 ) + $(n, n, r 3 ) ln(rfr 3 )£> 2 (n, ^ 



„22\ 



(3.59) 



The second part contains differential operators acting on D 2 of an Hamiltonian type. It has 
a form 

9 2 N C 



a\5 (n 2 ) + a3<5 (r 23 ) 



where 



ai = d't In r 2 3 df 2 + ln(-a?) = ^ ai 0{ 



2 ~ Q-2 



(3.60) 



(3.61) 



and a 3 is obtained by interchange 1 <-> 3. As shown in the section 1.3 the operator a\ can 
be transformed into a different form more convenient to analyze its properties under inversion. 
Indeed, in the complex notation, one writes the expression: 



di = In ri 3 + In k\ + c.c 



(3.62) 
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A useful relation, between pseudo-differential operators, 

In ri3 + In k\ = ln(r 2 3 k%) — k^ 1 In ri3 k\ 
can be deducted from the eq. (|1.43j ). On applying it one finds 

a\ = k\ ln(rf 3 A;i) ky 1 — lnri3 + c.c 
Since under the inversion we have 

In rf 3 k\ — > lnrf 3 k\ — 21nr 3 

the change in a\ is 

Aoi = —2 In r 13 + ln(rir 3 ) + c.c. = In 



and the variation of the last part of the integral I4, which comes from ( 3.60| ), is 



Ah 



8tt 2 



ci 2 ri(i 2 r3[$(ri,ri,r 3 ) - $(ri, r 3 , r 3 )l £> 2 (n, r 3 ) In 



(3.63) 

(3.64) 

(3.65) 
(3.66) 

(3.67) 



Summing all the changes one has A/2 + A/ 3 + A/4 = and this proves that the total integral 
/ is indeed invariant under inversion and, hence, under a global conformal transformation. 



3.4 Double pomeron exchange and triple pomeron interaction 

We already know the general form of the solution of the 4-gluon equation for the diffractive 
amplitude 

(3.68) 



S m DP = D% + ZD 2 + g 2 N c (V 12 + V U )D)- 



,(o) 



(o) 



It can be constructed as a sum of the DPE and a TPI parts as in eq. ( 3.38 ) from the sum of 
various inhomogeneous terms. 

However a more detailed analysis is instructive. Let us separate from this known exact 
solution some arbitrary function /, which may depend on the angular momentum j: 



DT(j) = f(j) + Df\j) 
This choice leads to an equation for the new 4-gluon function 



n(0) 



S40DP = L>$ + ZD 2 - [S 40 ~ 9 2 N C (V 12 + V34)]/ + 9 2 N C (V 12 + V U )D ( ; 
where the total inhomogeneous part has changed 



D$ + ZD 2 - D%> + ZD 2 -X , X = [540 - 9 Z N C (V 12 + V34)]/ 



>(0) 



(3.69) 



(3.70) 



(3.71) 



This procedure seems quite trivial but, as we shall see, gives some insight into the structure 



of the eq. ( 3.68; ) and of its solution. To this end we choose / to be the BFKL function 
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which depends on some gluon momenta. The idea of separating a reggeization piece which 
correspond to BFKL pomerons has already been used in [2C] and we shall try here to see in 



detail how it works; in this way we shall also be able to compare the results obtained in the two 
approaches. The separation of this simple term / in the amplitude leads to the subtraction of a 
new inhomogeneous term X which acquires the structure of a triple pomeron term. This means 
also that one can calculate some specific triple pomeron contributions expressing them in terms 
of simple functions. 

Let us see how this procedure works in some important cases, differing for the momenta 
which enter in the separated BFKL function. We shall find that, with a particular simple choice 
for /, one is able to completely cancel the D^j which generates the DPE contribution. In the 
following, for simplicity, only the forward case is considered. 

Let / = {l/2)g 2 D 2 (l) = (l/2)g 2 D 2 (2 + 3 + 4). One has 

X = y [S 40 - 9 2 N C (V 12 + V M )]D 2 (1) (3-72) 

Using eq. ( |3.6| ) we can express the j — 1 term in £40 in terms of the forward BFKL interaction 
Vo and u> and have 

2 2 

X = 9 -D w + 9 -[g 2 N c {V Q - V 12 - V u ) + 2w(l) - 5>(i)]D 2 (l) (3.73) 
Making use of the relation 

(V -V 12 )D 2 (1) = -W 2 (l,2,3 + 4) (3.74) 
where W 2 is defined by ( |3.12 ) and Q3.3| ), applying the bootstrap for the (34) gluon system 



g 2 N c V u D 2 (l) = 2(w(3 + 4) - w (3) - u(A))D 2 (l) (3.75) 

and introducing the function G(l,3) with two arguments |2C], which coincides with G(l,2,3) 
for the forward direction case (1 + 2 + 3 = 0), we finally obtain 

q 2 

X = y[ J D 20 (l) + G(l,3 + 4) + J D 2 (l)(^(3)+^(4)-2u;(3 + 4)) + D 2 (l + 2)(c«;(2)-w(l + 2)] (3.76) 

Thus we find that the changed function will satisfy an equation with a new inhomogeneous 
term ^ 

D40 + ZD 2 + |-[-£>2o(l) " G(l, 3 + 4) 

- L> 2 (l)(cu(3) + w(4) - 2w(3 + 4)) - D 2 (l + 2)(w(2) - u(l + 2))] (3.77) 

Note that the additional term — (l/2)g 2 D 2 Q will cancel the identical term in . As a result, 
one has converted the double pomeron exchange contribution coming from (l/2)g 2 D 20 (l) into 
a triple pomeron contribution corresponding essentially to G(l, 3 + 4) and in the same time has 
explicitly separated the term (1/2) g 2 D 2 (l) in the amplitude. In other words, one can calculate 
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the triple pomeron contribution corresponding to a vertex 
Q 2 

y [-G(l, 3 + 4) - D 2 (l)(w(3) + w(4) - 2kj(3 + 4)) - D 2 (l + 2)(w(2) - w(l + 2))] 

as a sum of the double pomeron exchange coupled to — (l/2)g 2 L>2o(l) and a term (1/2)<7 2 L>2(1)- 
Evidently this result is trivially generalized for / = D2(i), i = 2,3,4 by simple permutation 
of indexes 1,2,3 and 4 and one may note that the separated terms found for the amplitude have 
the same momentum dependence of the corresponding zero order cancelled terms. Let us see if 
this is true also for the last three terms in ( 3.21| ). 



Thus we consider the case / = {l/2)g 2 D 2 (l + 2) for which we have 

f a 2 
2° 20 + 1 

In terms of W we have 



X = y -D 2Q + y -[g 2 N c (V - V l2 - V u ) + 2w(l + 2) - £ u(i)]D 2 (l + 2) (3.78) 
V Q D 2 (l + 2) = -W 2 (l + 2,0,3 + 4) (3.79) 



The bootstrap gives 

9 2 N C (V 12 + V U )D 2 (1 + 2) = 2^(1 + 2)(2o;(l + 2) -£>(*)) (3.80) 

so that in terms of G we obtain 
a 2 

X = y [-Z?2o(l + 2) + G(l + 2, 3 + 4) - D{1 + 2)(4w(l + 2) - ^) "(*))] ( 3 - 81 ) 

Again we see that the term with the double pomeron exchange coupled to g 2 (l/2)D 2 o(l + 2) 
can be transformed into a triple pomeron vertex, essentially, into —G(l + 2, 3 + 4) term. 

Finally we study a more complicated case with / = (l/2)g 2 D 2 (l + 3). In this case we find 

Vo£> 2 (l + 3) = -1^(1 + 3,0,2 + 4) (3.82) 

Calculation of V\ 2 or V34 applied to -D2(l + 3) is done using the formula (see appendix A. 4) 

V 12 D 2 (1 + 3) = W 2 (2 + 4, 1, 3) + W 2 (4, 2, 1 + 3) - W 2 (3, 1 + 2, 4) - W 2 {1 + 3, 0, 2 + 4) (3.83) 

and 

^34^2(1 + 3) = W 2 {2 + 4, 3, 1) + W 2 {2, 4, 1 + 3) - W 2 (l, 3 + 4, 2) - W 2 {1 + 3, 0, 2 + 4) (3.84) 

Using these results one has 

f 
2 



X = 9 - 



D 20 (l + 3) + G(l, 2 + 4) + G(2, 1 + 3) + G(3, 2 + 4) + G(4, 1 + 3) 



-G(l, 2) - G(3, 4) - G(l + 3, 2 + 4) - D 2 (l)(u(3 + 4) - w(3)) 
D 2 (2)(o;(3 + 4) - w(4)) - £> 2 (3)(w(l + 2) - cj(1)) - D 2 (4)(w(l + 2) - w(2))l (3.85) 
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The result for / = (l/2)g 2 D 2 (l + 4) is obtained from this after the permutation of 3 and 4. 

Inspecting these results and comparing them with the form of our triple pomeron vertex, we 
see that only four terms 

^ [G(l, 3) + G(l, 4) + G(2, 3) + G(2, 4)] (3.86) 

are not changed under these transformations and thus correspond to a true triple pomeron 
interaction. All the rest can be transformed into terms which are essentially double pomeron 
exchange contribution. Conversely, one can eliminate terms from the double pomeron exchange 
substituting them by equivalent triple pomeron contributions. 
The most radical result follows if one takes 

f = Dg\D 20 ^ D 2 ) (3.87) 

In this case all the double exchange becomes cancelled and the whole amplitude is given by a 
sum of two terms (in an evident symbolic notation) 

Df = D$(D 20 - D 2 ) + [ G 2 (Y - y)G 2 {Y - y)ZD 2 (y) (3.88) 

Jo 



with a new vertex 

ZD 2 = 9 - 
2 



G(l, 3) + G(l, 4) + G(2, 3) + G(2, 4) + G(l + 2, 3 + 4) 



- G(l, 3 + 4) - G(2, 3 + 4) - G(3, 1 + 2) - G(4, 1 + 2) (3.89) 
Comparing to (|3.45|) for the forward case, we observe that it coincides with the part V(1234) of 



the vertex introduced in [20|, being its leading part in the large N c limit. 



3.5 Relation to the dipole approach 

According to the previous section we shall consider the DPE part completely transferred into 
the TPI one and study the triple pomeron vertex in the simpler form ( 3.45| ) for a non- forward 



direction. It is important to note that the coupling to BFKL pomerons of this vertex will lead 
to further simplifications. In fact in the coordinate space of 4 gluons, the dependence on only 
the sum of the momenta of two gluons, say, 1+2, is translated into a factor 5 2 (ri 2 ), so that the 
two gluons have to be taken at the same point. Since the wave functions ^{ri,r 2 ) and ^(r^, r±) 
of the two final pomerons, coupled to the vertex, vanish if r\ = r 2 and r% = respectively, 
all terms in eq. (|3.45| ) which depend either only on the sum 1+2 and/or only on the sum 3+4 
give zero, coupled to the two pomerons. Thus the effective triple pomeron vertex coupled to 
pomerons reads 



ZD 2 = 9 - 
2 



G(l, 2 + 4, 3) + G(l, 2 + 3, 4) + G(2, 1 + 4, 3) + G(2, 1 + 3, 4) (3.90) 
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Since both pomeron functions ^f(ri,r 2 ) and ^f(r3,r^) are symmetric in their respective argu- 
ments, due to the positive signature of the pomeron, all terms in ( 3.9C| ) give identical contribu- 
tions and one has 

ZD 2 = 2 5 2 G(l,2 + 3,4) (3.91) 

Moreover let as consider in the coordinate space the "improper part" of the function G(r\, r 2 , ^3), 
which includes all the terms proportional to H 2 {r\2) or /and <5 2 (r23); since dHOD is also propor- 



tional to 5 2 {r 23 ) we find that in the "improper part" at least three gluons, either 123 or 234, are 
to be taken at the same point in the transverse space. Then these terms will vanish due to the 
mentioned property of the pomeron wave function. 

Thus only the contribution of the "proper part" is left and it is explicitely given by 

n 2 AT 2 

G pr (n,r 2 ,r 3 ) = -L« « ^(n,^) (3.92) 

S7r r 12 r 23 



(with the regularization ( [3.51 ) implied). 



Coupling this triple pomeron vertex to the two final pomerons, one obtains the following 
expression for the triple pomeron contribution to the diffractive (non-forward) amplitude: 



g 1 N f r r 
Dl PI (Y) = -?—£ / dy d 2 r 1 d 2 r 2 d 2 r 3 ^ T D 2 (r l ,r 3 ;y)^ l (r 1 ,r 2 ;Y - y)V 2 (r 2 ,r 3 ;Y - y) 

47T Jo J r 12 r 23 

(3.93) 

If one writes the initial pomeron amplitude D 2 via the non- amputated function ^, which is the 
conformal invariant one, the amplitude reads 

Dl PI (Y) = - 9 -^ f dy I d2r f^ fr l(ri , r2 . y _ y )^ 2 (r 2 , r 3 ; y - y^fffflKn, r 3 ; y) 
4^ 3 Jo J rf 3 rf 2 r^ 

(3.94) 

In this form it is evident that the triple pomeron vertex is not symmetric with respect to the 
initial pomeron and two final ones: there appears an extra operator rf 3 dfd 3 acting on the 
initial pomeron. Note that this operator is essentially a product of the Casimir operators of the 
conformal group for the holomorphic and antiholomorphic parts. Therefore it is convenient to 
expand the functions in terms of the conformal basis, introduced in the first chapter. Using 
( 1.56| ) we do this for all the three pomerons in ( p. 94 ) and compute the action of the Casimir 



operators by means of ( |1.48| ). We then obtain, after the integration over y 
Dl PI (Y) = - 9 -rf E <H*o></.i|*io><H*20> e 



4tt 8 f d 2 nd 2 r 2 d 2 r 3 
/ 2 2 2 — E ^A r i,r 2 )E^ 2 {r 2l r 3 )E^{r u r 3 ) (3.95) 

0"n—l,vO-n+l,v J r 13 7 'l2 ? *23 

In this form the triple pomeron contribution can be compared to the double dipole density 



found by Peschanski [31| in A.H.Mueller's colour dipole approach. One observes that the two 
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expressions differ only in the sign and by a factor 

4tt 8 

which in our approach distinguishes the initial pomeron from the two final ones. The integral over 
the coordinates of the three pomerons is the same. So essentially the three-pomeron contribution 
to the diffractive amplitude found in our approach coincides with the double dipole density in 
the dipole approach. 

However one should not forget that in our s-channel unitarity approach the TPI term ( p. 95 ) 
does not exhaust all the diffractive amplitude. In fact we have seen that 

of = D m (D 20 D 2 ) + Dj PI (3.96) 

At high energies the TPI term behaves essentially as s 2A and the first one as s A where A is the 
BFKL intercept. So one might think that the first term could be neglected. However the correct 
region of the validity of the leading log approximation, implicit in the hard pomeron theory, 
is g 2 In s ~ 1 when the two terms in (80) have the same order of magnitude. Moreover the 
BFKL separated term goes as 1/N C and the TPI as g 2 N c and also these factors are of the same 
order. The dipole approach uses essentially the same leading log approximation but it leads 
to the double dipole density which coincides only with the TPI term in the s-channel unitarity 
approach and shows no trace of the first term. Therefore this fact points to certain differences 
between the two approaches. 



3.6 Higher order densities in the dipole approach 

Keeping in mind this substantial agreement between the results obtained in the s-channel unitar- 
ity approach and in the colour dipole picture for the two pomerons (dipoles) distribution, let us 
investigate the general structure of the higher order dipole densities given by the latter. To see 
if the previous equivalence found is maintained, one should study the multipomeron distribution 
in the former approach. In particular the first step would regard the study of the structure of 
the system of up to six reggeized gluons to look at the transition from 1 to 3 pomerons. In fact, 
as we shall show, this coupling does not directly appear in the standard dipole picture which 
instead leads to a set of fan diagrams with a single "elementary" three-pomeron coupling. If in 
such a context a future analysis will show the presence of pomeron vertices of order higher than 



3 then an extension of the dipole picture, for example as the one proposed by Peshanski [31], 
should be required. 

Following the approach of A.H. Mueller, as briefly reviewed in section 2.4, in the colour 
dipole formalism the k-fold inclusive dipole density is obtained as the fc-th functional derivative 



of the functional D{u(ri, f/)}, taken at it(rj,r/) = 1 [11]. The arguments rj and rf are the 
dipole endpoints in the transverse plane. We remember here the simple equation the functional 
D satisfies 

D(r 1 ,r ,y, U )=u(r 1 ,r )e 2 y^ 
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+ 1HT / dy'e 2 ^-y>^ / d 2 r 2 ^^D( n ,r 2 ,y' \ U )D(r 2 ,ra,y' » (3.97) 
8vr 3 7 7 rf 2 r| 

Here ri and ro are the end points of the qq pair which determine the initial dipole; moreover we 
define y = In z%o where Zxo is the ratio of the scaling factors of the antiquark and of the quark; 
u){r) is not a Fourier transform of the trajectory, but just u(k) with k/m formally substituted 
by r/e, where e is an ultraviolet cutoff. This cutoff is also implied in the singular kernel of 
the integral operator in r 2 . Let as note that eq. ( 3.97| ) is compatible with the normalization 
condition D(u = !) = !. 



The k-th derivative of eq. ( p. 97 ) give an equation for the fc-fold dipole density if computed 
for u = 1. For k > 1 we obtain 



71 k 



„2 AT ny p 

(ri,r ; pi, ...p*; y) = ^ J dy'e 2 ^>^ j d 2 r 2 -^n k {r u r 2 ; Pl , ...p k] y') + ( n m) 



'12' 20 



s 2 ^ r y 



k-l 



8tt 3 



f dy' e 2(^'V(no) f d \ 2 ^_ ^( ni(ri)r2 . pi> y >^(r 2 , r Q ; ...p fc ; J/') 

"'O J r 12 r 20 ; = i 

+ symmetrization terms) (3.i 



where the symmetrization terms (ST) are obtained from the explicitly shown one by taking all 
different divisions of arguments pi,...p k into two groups with I and I — k arguments. In the 
k = 1 case an inhomogeneous term whose form is clear from ( |3.97 ) is present. One should 
note that the operator on the right-hand side acts nontrivially only on the first argument of the 
density n k . Its action on the rapidity variable y, on the contrary, is rather simple and it is better 
treated in a transformed associated equation. Thus it is convenient to multiply the equation by 
e -2yu)(r w ) ^ differentiate then with respect to y and pass to the j-space by the standard Mellin 
transformation. After having performed these operations one obtains 



(j - l)n fc (ri,r ;pi, ...pk'J) 
,9 2 N C 



9 2 N C 



47T 3 

djidj 2 



d 2 r 2 L(r 1 , r 2 , r 20 )n k (n,r 2 ; p 1 , ...p k ;y) 



d 2 r 2 ^\ 

' 12' 20 



8^ 3 J (27Ti) 2 (j + l-Jl -h) 

53(n{(ri,r 2 ;pi, ...pr, ji)n k -i(r 2 ,r ; p l+1 , ...p k ;j 2 ) + ST) 



k-l 



where we introduced the BFKL kernel in the coordinate space 



£(ri2,r 2 o) 



10 



(r 2 2 + e 2 )(r 2 o+e2) 



27 rln^ 



5 2 (ri 2 ) + 5 2 (r 20 ) 



(3.99) 



(3.100) 



Comparing with ( |3.51| ) we see that it does not depend on e and is ultraviolet stable. The solution 
of the equation fl3.99|) is quite easy to obtain if one describes the dependence of the densities on 
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their first two arguments in terms of the conformal basis: 



(3.101) 



Here we have suppressed all other arguments in n k irrelevant for the time being. The densities 
in a given conformal state are obtained from rik(ri,ro) by the inverse transformation which 
follows from orthogonality property ( 1.52 ) and a relation between E n ^ u and E- n - u (see 

d 2 r\d 2 r Q 



n r i 



■K(ri,roK(ri,r ) 



(3.102) 



10 



Therefore, to pass to the conformal basis, we integrate eq. ( |3.99| ) over r\ and rn as indicated 
in ( |3. 102 ) . The first term on the right-hand side can be simplified due to the property of the 
BFKL kernel 

9 2 N C 



J d 2 r 2 L(r 12 ,r 20 )E ll (ri,r 2 ) = u^E^ri, 
Therefore after the integration we obtain 

9 2 N C f dj x dy. 



(3.103) 



(j - 1 -^^)n^(Pi,...p k ,j) 



8vr 3 



(2TTi) 2 (j + l-j 1 -j 2 ) 



d 2 r\d 2 r 2 d 2 TQ 

2 2 2 
'12'20'10 



k-1 



E t( r w) ^2( n i( r i,r 2 ; pi, ...pi;ji)n k _i(r 2 ,r ; pi + i, ...pk]j2) + ST) (3.104) 



i=i 



To find the final form of the equation we have only to present also the densities ni and n^-i as 
3t arguments in the form ( 

(j - l-u lx )n^(pi,...p k ,j) 



functions of their first arguments in the form (3.102). Then one has 

dj\dj 2 



(2TTi) 2 (j + l-j 1 -j 2 ) 

v ^i^^{ n ^ 1 {pi^--Pi\h)n k ^_ l {pi+i,-.-Pk]j2)+ ST) 



k-l 



1=1 



where 



V, 



g 2 N c f d 2 nd 2 r 2 d 2 r 



8vr 3 



2 2 2 
to 1 on' 



^(no)^i (n 2 )E^(r 20 ) 



(3.105) 



(3.106) 



12' 20' 10 

is just one half of the three-pomeron vertex introduced by Peschanski. The general structure 
of this equation allows one immediately to imagine a tree-like expansion of the colour dipole 
densities as we shall show more in detail. At each tree node a vertex ( |3.106] ) will be present. 
The integral of its explicit expression has been calculated |2fi, 58] with the leading quantum 
numbers n = v = for the related three pomerons and it turns out to have a quite huge value 
(7766.679). Of course there will be present some correcting factors but still this means that such 
a perturbative expansion in pomeron vertices may need a complete resummation. 



Eq. ( |3.105| ) allows to obtain successively dipole densities for any number of dipoles starting 
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from the lowest order one-dipole density, for which 

= ■ E f P) \ (3-107) 

(we recall that in this notation p includes two endpoints of the dipole pi and pf, p 2 = pfj). 

Putting this into ( 3.105 ) for k = 2 and integrating over j% and 22 we arrive at the expression 
obtained by Peschanski 

<(P!,p 2 ;i) = E ~ ,/ 7, v ^,KM)KM)4-a (3-ios) 

where w = j — 1. (To compare with [31] one should take into account that factors l/(2a^) are 
included in the definition of sums over /i's in our notation). 

As a further step we study the density for three dipoles. Eq. ( |3.105] ) for k = 3 reads 



djidj 2 



(2vrz) 2 (i + 1 - Ji - J2) 
E V flM (n1 1 (Pi;hH 2 (p2,p 3 ;j2)+ ST) (3.109) 

Ml ,^2 

It is instructive to study explicitely this term inserting the expressions for n^ 1 and n^ 2 obtained 
earlier. Then we get, after integrations over j\ and j 2 '- 

n%(pl,P2,p3-,j) = ,,_,, E V MllMt V WM>i 



LO — U>u 



7 — i E ; i (p 1 )E; 3 (p 2 )E; A ( P3 )^ rT (3.110) 

To this term we have to add terms which symmetrize in the three dipoles. 

Studying ( |3.110j ) we see that it corresponds to the picture when first the initial pomeron 
splits into two pomerons, 1 and 2, and afterwards the pomeron 2 splits into pomerons 3 and 4. 
One does not find here a local vertex for the transition of the initial pomeron into three final 
ones. It is not difficult to see under which condition one would get such a local vertex. If we 
forget about the dependence of the second denominator on then one can sum over Using 
the completeness relation ( |1.50| ) one obtains 

^ 2 \ / r 12 r 10 r 23 r 30 

(3.111) 

which is just the vertex from one to three pomerons introduced by Peschanski. However, the 
described summation is not possible due to the second denominator. It implies that the pomeron 
2 has to evolve in y from the point of its formation from the initial pomeron up to the point of 
its splitting into the final pomerons 3 and 4. 

Thus our conclusion is that the vertex for transition from 1 to k pomerons introduced by 
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Peschanski, in fact, does not appear in the solution of the Mueller equation ( |3.98| ) for the fc-fold 
density, which rather corresponds to a set of all fan diagrams with only the triple pomeron 
coupling. Absence of higher-order couplings can be directly traced to the structure of the 
equation ( |3,97| ) for the generating functional D, quadratic in D, which in fact means that the 
basic vertex changes one dipole into two. 

The general structure of higher order vertices introduced by Peschanski hence refers to an 
extension, which he proposes, to the dipole model. In his work the possibility of a splitting of 
one dipole into many dipoles is considered, going therefore beyond the standard dipole model. 
His idea is related to the fact that in general a contribution to the transition from 1 to p pomeron 
given by the fan diagrams would require a large rapidity interval at each splitting step in order 
to ensure the validity of the Regge approximation. The derived effective field theory is therefore 
more complex but it should be characterized by a strong symmetry since the vertices are related 



each other by integral relations, as for example the one in eq. (3.111). 



A general analysis of the n-reggeon amplitudes approach, based on s-channel unitarity, is 
not known yet, even if some work in this direction is in progress. Any advance in this direction 
will be of great interest to understand the Regge limit effective behaviour of QCD. 

To conclude this section we shall try to analyze in more details the picture one derives from 
the triple pomeron interaction structure so far verified. It has been shown that at asymptotic 
energies the higher-order densities in the dipole approach correspond to the standard Regge- 
Gribov picture, in the tree approximation (fan diagrams), with only the triple pomeron inter- 
action, which however has a highly complicated non-local form. Indeed, the triple pomeron 
interaction present in ( 3.110| ) corresponds to a structure 



g 2 N c f d rid r 2 d r 3 ~ , , ~ , , ~ , , . , 

T = -5-3-/ 2 2 2 — G 3{n,r 2 ;r 1 ,r 2 )G l {r 2 ,r 3 ;r 2 ,rr i )G 2 {r 3 ,r 1 ;r 3 ,r 1 ) (3.112) 



where Gj, i = 1, 2, 3 are Green functions of the interacting pomerons defined as 

G(n,r 2 ; r , r' 2 ) = £ M V ^ lJ 2> (3.113) 

They are not the physical BFKL Green functions. The latter include an extra factor depending 
on fj,: 

G(r 1 ,r 2 ;r l ,r 2 ) = > a n+1 „a n -i v <=■ 3.114) 

However in the limit s — > 00 only the lowest conformal weights contribute n = v = for which 
0"n±l,v = 27T 4 and (3.113) and (3.114) coincide. Then we can forget about tildes in (3.112). 



We transform the Green functions to given total momenta of the pomerons presenting 

r d 2 ] 

G-M^AA) = J j^e il ^-^G l3 (r 12 ,r' 12 ) (3.115) 



where R3 = (1/2) (n + r 2 ) and similarly for the two other Green functions. Introducing R = 



3.6 Higher order densities in the dipole approach 
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r i + r 2 + r 3 we transform the integration over the coordinates as follows 

J d 2 r x d 2 r 2 d 2 r 3 = J d 2 Rd 2 r 12 d 2 r 23 d 2 r 31 5 2 (r 12 + r 23 + r 31 ) (3.116) 
The coordinates themselves are r\ = (1/3) (i? — r 2 i — r 3 i) etc., where from we find 

Ri = (l/6)(2i? - na - ri 3 ), #2 = (l/6)(2i2 - r 21 - r 23 ), fi 3 = (l/6)(222 - r 31 - r 32 ) (3.117) 
and 

3 3 

*E^i = i ( 1 /$) R Y,h -*(l/6)(ri 2 il2 + r23i23 + r 3 il30 (3.118) 

where we denoted l\ 2 = h — l 2 etc. The integral over R gives 9(2vr) 2 5 2 (/i + h + h)- Presenting 
the remaining 5-function in ( |3.116| ) as an integral over an auxiliary momentum q we find an 
expression (for fixed l\,l 2 and Z 3 ) 

g 2 N c f d 2 g d 2 r l2 d 2 r 23 d 2 r 3l ( 1 1 1 \ 

"8^ J W rfr&l 6XP r i2(<? " 6^ + ir23{q ~ 6 h3) + tr3l( * " 6 kl ) 

G h (r 12 ,r[ 2 )G h (r 23 ,r' 23 )G h (r 31 , r' 31 ) (3.119) 

At this point we recall the expression for the BFKL Green function with a fixed total mo- 
mentum: 

GlM = j^j { / + % )2 ^ U) E2\r)E^{r') (3.120) 

where 

mr) = j<fRe,p(UR)( lR + r/m _ rm ) (3.121) 

and where we retained only the dominant isotropic term. At s — ► oo the vicinity of v = gives 
the dominant contribution. If I ^ then the integral in ( p.!21| ) converges at large R and we can 
take the functions E out of the integral over v at v = 0. Taking then the asymptotics of the 
remaining integral, we find 

G iM - i sA (^72^ )(r) ^ )(r,) (3 - 122) 

where A = uj n= o tV= o is the BFKL intercept and a = 7g 2 N c ((3)/(2ir 2 ). As we see, the Green 
function asymptotically factorizes in r and r'. This means that we obtain a quantum field theory 
of pomerons with a propagator 



(not really depending on the momentum /) and an interaction vertex 

xrn . , . 9g 2 N c f d?q d 2 r 12 d 2 r 23 d 2 r 31 
V (.li, 12,13) 



8tt3 J (2tt)8 r 2 2 r 2 3 r 2 



31 
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Unitarization in the large N c limit: the three pomeron vertex 



exp (ir 12 {q - h l2 ) + ir 23 (q - ^ 23 ) + ir 31 (« - jUl) (ri 2 )£#° (r 23 )^ 2) (r 3 i) (3.124) 

The vertex factorizes under the sign of the integration over q: 

Q a 2 N f d 2 a 111 
V(h,h, h) = J jAiAh, Q - g*12) Ah, Q - - Q hz)J(h,q ~ -hi) (3.125) 



d2r Aarr^Q),^ f 



where 

J "-'» = J (2^^ E »"W = J ^- 9 + l /2|'|P-«-</2| < 3 - 126 » 
Note that for I = this derivation is incorrect. Calculations show that in this case eq. 
fl3.125j ) for the vertex remains valid with 

J(Q,q) = ^ (3-127) 
However the Green function ( 3.120| ) at I = has an asymptotics 

G °^'>^ A SM-^) <3 - 128) 

which is evidently no more factorizable in r and r'. 



Chapter 4 

Conclusions 



As stated in the introduction, the thesis has been devoted to the study of some issues relevant for 
the hard pomeron model. Both theoretical and phenomenological aspects have been considered. 

From a phenomenological point of view in the second chapter an extension to include the 
running coupling for the hard pomeron is studied. The leading singularity in this way is softened 
and once the infrared sector is safely included in the model one can avoid also the problems 
due to the diffusion mechanism which spoil out the validity of the pure BFKL pomeron. The 
model, which preserves the gluon reggeization property by means of the bootstrap condition 
and the right gluon distribution behaviour in the DLLA, has been studied both numerically 
and in some aspects analytically. Numerical calculations were addressed to the computation 
of the spectrum and eigenfunctions of the Schrodinger-like equation for the pomeron. The 
spectrum has shown a discrete part and for the wide parameter range two negative "energies" 
states corresponding to two supercritical pomeron states. The slope of each singularity has been 
calculated with a perturbative approach. The sensitivity of the model is mainly with respect to 
only one infrared parameter, a mass which we chose to fix forcing a value for the slope of the 
leading singularity. We have applied such a model to the process of two virtual photon scattering 
and to see some unitarization effects we have considered in the large iV c limit a multi-scattering 
(eikonal) formulation which have shown small unitarity corrections of such a kind to the cross 
sections. Triple pomeron interactions, required by unitarity, are not taken into account in these 
calculations and larger unitarity effects are therefore not excluded. Similar considerations were 
applied to a toy model, with a Gaussian colour dipole distribution, for the proton. 

An application of this model to the inclusive jet production has been considered. A correct 
asymptotic behaviour of the differential cross section for large jet transverse momenta has been 
found. In this case one must note that in the BFKL approach already some unitarity correction, 
which have lead to modify the BFKL equation with a non-linear contribution generating fan 
diagrams, as well as coherence effects were, for example, invoked. Since the leading singularity is 
a pole, the y dependence in the cross section depends on the inclusion of the non-leading states 
of the spectrum. Here we calculated the contribution due to the two supercritical pomerons and 
for a better estimate a much larger part of the spectrum depending on the energy, should be 
considered. 

The contribution from the whole spectrum have been included in the analysis of the evolu- 
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tion of the gluon density and of the proton singlet structure function. The initial proton gluon 
distribution in rapidity has been found by matching the related singlet structure function. The 
main problem has been to find a way to couple the gluon to the photon taking in account the 
introduced running of the coupling. A phenomenological approach has been chosen, based on 
the requirement that for the structure function in the DLLA the analytical asymptotic behaviour 
should match the one calculated using the DGLAP quark distributions. The numerical calcu- 
lations have still shown a steep increase with decreasing x of the gluon distribution, but with a 
better behaviour than for a BFKL pomeron. 

In the theoretical part of the thesis we deal with the unitarity problem. At the basis of this 
contribution are the works of J.Bartels et al. where a program to ensure unitarity requirements 
has been developed and the analysis of the lower order corrections, up to 4 gluons in the t-channel, 
started. The idea of studying these corrections in the large N c limit, a program initiated by 
M.Braun, is pursued in order to avoid the otherwise complicated colour structure. The nice 
property of this limit is that one does not loose many important features of the underlying 
interactions and it is also possible to directly compare the obtained results with some recent 
developments derived in the colour dipole picture. In particular the system of four gluons has 
been studied for a non-zero momentum flow in the t channel in the leading and next-to-leading 
order in l/N c . The solution has been shown to have in general double pomeron exchange and 
triple pomeron interaction contributions. The latter is written in terms of a three point function 
which is by itself infrared stable and conformal invariant. Following the idea of Bartels and 
Wuestoff it is shown the possibility of eliminating the DPE terms in the solution so that the 
unitarity correction to this order may depend only on a particular three pomeron vertex, plus 
some single pomeron terms. There is complete agreement with the Bartels et al. works once 
the large N c limit is taken into account. After coupling this vertex to pomerons one can easily 
compare the obtained two pomeron amplitude with the double dipole colour density analyzed 
by R.Peschanski in the standard dipole picture . One finds a substantial agreement for the triple 
pomeron interaction form, apart from an asymmetry factor and the previous mentioned single 
pomeron terms. The properties of higher order dipole densities, as defined in the A.H. Mueller 
formulation of the colour dipole model, are subsequently studied. A picture is found where in 
the rapidity evolution a pomeron may split into two pomerons, thus leading to a fan diagrams 
description with only three pomeron vertices involved. The resulting pomeron effective field 
theory would hence have a strong signature. Therefore a corresponding analysis is needed in the 
reggeon approach based on s-channel unitarity. The analysis starting from a system of 5 and 6 
reggeized gluons and in general of higher order distributions will show if this picture is correct 
or one has to deal with higher order effective pomeron vertices. In such a case one should also 
consider extension in the colour dipole picture, as for example proposed by R.Peschanski. In any 
case the structure of a possible effective theory of QCD in the Regge limit will then be strongly 
restricted. 



87 



Acknowledgements 

I would like to express my deep gratitude to my supervisor Professor Giovanni Venturi for 
his support, interest and encouragement and for countless stimulating discussions in different 
fields of Physics. 

I am sincerely and infinitely grateful also to Professor Mikhail Braun, for his guidance and 
help in the last years. His ideas have had a strong influence and most of the work present in the 
thesis is inspired and arises from his collaboration. 

It was a great experience to visit the University of Santiago de Compostela and I would like 
to thank Profs. C. Pajares and L. Miramontes for their hospitality. 

Furthermore I thank Prof. R. Peschanski for pleasant and useful discussions. 

Finally a special thank to Fabio Finelli for many discussions and a fruitful collaboration. 



88 Conclusions 



Appendix A 



A.l Asymptotic behaviour of the wave functions. 
A. 1.1 The momentum space behaviour 

Let us consider the asymptotic expression for the gluon Regge trajectory given in ( |2.11| ) 

u(q) ~ rhi- ^ (A.l) 

v V»™ 2 26 lnm 2 v ; 

Hence, for N c = 3, one can write the asymptotic form of the pomeron non-homogeneous equation 
at high momentum for the semi-amputated wave function ip 

lng 2 am 2 lnm 2 If d 2 q'ip(q') 

l n i oviq) H T\ o— = e W\l) H / 77 91 : r , 7775 (A. 2) 

mm z q z mq z tt J [[q — q') z + m z \ ln[{q — q') z + 771/ J 

where we have put E = (3/6)e (we consider the case m = mi for simplicity) Let us study the 
last (integral) term, which we denote by D. Changing the variable q' = \q\n we present it in the 
form 

_ 1 f d 2 Ki;(\q\K) 



7T J [(n — k) 2 + m 2 /q 2 ] ln[q 2 (n — k) 2 + m 2 ] 

= 1 f d 2 ^{\g\{n + K)) 

it J (k 2 + m 2 /q 2 ) \n(q 2 K 2 + m 2 ) 

with n 2 = 1. In the last form it is evident that the leading terms at q — ► 00 come from the 
integration region of small n. Then we split the total k space into two parts: k > kq and k < k 
where kq is a small number 

k << 1 (A.4) 

The contributions from these two parts we denote as D\ and D2, respectively. Our first task is 
to show that D2 cancels the kinetic term in (A. 2) irrespective of the asymptotics of ip(q)- 
With small enough kq and for a "good enough" tp 

ip(\q\ (n + «)) ~ ip(\q\n) = ip(q) (A. 5) 



so that we can take it out of the integral over k in Z?2- Actually ( |A.5| ) is our definition of a 
"good" function, so that we shall have to check if this condition is indeed satisfied for the found 
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asymptotic ip(q). 



With flAg ) D 2 simplifies to 



° 2 = / T7T X 2 m ; ,2 I 57 = ^(9) In 7 2 ( A - 6 ) 

vr Jo (q' +m 2 )ln(q' + m 2 ) mm z 

As q — > 00 we assume that also 

Kog — ► 00 (A. 7) 

Evidently this condition fixes the manner in which kq goes to zero as q turns large. This should 
be taken into account when verifying condition (A. 4). Then (|A.6|) gives 



D 2 =m^ l ^r ( A - 

m m z 



We have also 



But according to ( [A.7| ) 



so that we have 



Then our final result is 



In In K^q 2 = hi (In q 2 + In Kq) 
Inq 2 » I lnKpl 



In In n 2 Q q 2 = In In q 2 + + 



In q 1 



(■ 



In g 



2 



The first term exactly cancels the kinetic energy ( first) term in the pomeron equation. The 
correction term in QA.9| ) evidently is much smaller than eijj(q), since the factor which multiplies 
the function ifi goes to zero in (|A.9|). So we can safely neglect it. 



As a result, the asymptotic equation becomes 

am 2 In m 2 



q 2 In q 2 

with the term D\ given by 



~e^{q)+D 1 (A.10) 



1 f d 2 g'j;(q')e((q - q') 2 - K 2 q 2 ) 
1 nj {q-q'Y\n(q-q'Y 1 J 



where we have omitted the m terms in the denominator because of (A.7). 



We consider the function x(q) defined by 

x{q) = ^{q)q 2 ^q 2 



(A.12) 



A.l Asymptotic behaviour of the wave functions. 
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Multiplying (A. 10) by q 2 lnq 2 we find an equation for x 



am Inm = e X (q) + -J {q _ q , ?Hq _ ql? 



At this point we make a second assumption. Namely we assume that in the integral term 
of QA.13; ) values m « q' « q give the dominant contribution (as typical for logarithmic 
integrals). Of course, this assumption is also to be checked for the final asymptotics. With this 
assumption, we can forget about the function and also put the last factor equal to unity in 
the integral term of ( |A.13 ). We obtain 



am 2 



lnm =e X {q)+ ^ (A.14) 

JO q hiq 

Differentiating with respect to q 2 



I— = X (A.15) 
dq 2 q 2 In q 2 



or 



<,X -X (A.16) 



d In In q 2 
with a solution 

X (q) = Ae W (J^f.) = A(\nq 2 )- 1 " (A.17) 



e 

The initial function ip(q) has the asymptotics 

^q) = ^{\nq 2 f (A.18) 
q 2 

where 

P = -l-l = -l-±. (A.19) 
e bE 

Now we must check that our assumptions are indeed fulfilled for the found asymptotics. 

Let us begin with the second assumption that the values m « q' « q give the bulk of 
the contribution to the integral in ( [A. 13 ). Evidently, in order that the integral be dominated by 



large values of q' , it should diverge as q — > oo. This leads to the condition 

E < 0, (3 > -1 (A.20) 
So our asymptotics can only be valid for negative energies, that is, for bound states. 

Now for the second part of this assumption. To prove that values q' « q dominate we 
shall calculate the contribution from the region q' » q and show that it is smaller. The 
corresponding integral is 

I = Aq 2 lnq 2 f°° dq' 2 (ln q' 2 )^ 1 /^ (A.21) 
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In terms of x = In q' 2 

roo 

I = q 2 lnq 2 / dxx^ 1 exp(-x) (A.22) 

Jhiq 2 

The integral over x can be developed in an asymptotic series in l/lng 2 : 

roo roo (hia 2 )^ ~* f°° 

/ dxx^" 1 exp(-x) = - / x /3 " 1 dexp(-x) = ^— ^ + (3-1) dxx^~ 2 exp(-x) = ... 

Jlnq 2 J\nq 2 T J\nq 2 

(A.23) 

From this we conclude that the integral I has the asymptotics 

I = A(\nq 2 f (A.24) 

to be compared to the contribution from the region q' << q which behaves as (lng 2 )^ +1 . We 
see that we have lost one power of Inq 2 , so that the region q' » q indeed can be neglected. 
Now to the assumption ( A.5| ). We have explicitly 



ip(\q\(n + k)) = (A/ q 2 ){n + k)~ 2 (In q 2 (n + k) 2 ) (3 = 

{A/q 2 )(l - 2n K - k 2 + 4(n K ) 2 )(ln q 2 f(l + (/3(2n« + K 2 )/\nq 2 ) (A.25) 

and it is evident that ( |A11| ) is satisfied with the behaviour of k as indicated in ( |A.7| ). Indeed 
take k ~ q~ s with 5 < 1. Then qn ~ q 1 ^ 5 — > oo and all the correcting terms in ( |A.25| ) have the 
order q~ s . 

Note that this does not mean that ( |A5| ) is quite obvious. It is not valid for, say, the 
exponential function. In fact we have 

exp(ag 2 (n + k) 2 ) = exp(ag 2 ) exp(ag 2 (2nK + k 2 )) 

and since qn is large the second factor cannot be neglected. 

In conclusion, we have verified that our assumptions are fulfilled and therefore the asymp- 
totics ( |A. 18 ) is correct for negative "energies". 



A. 1.2 The coordinate space behaviour 

We are also interested in the coordinate space behaviour of the semi-amputated wave function 

tf>(r)= r^MqrMq) (A.26) 



o 



2vr 

To estimate ip(r) for r — > we make use of the asymptotic momentum behaviour found 



previously (|A.18Q , so we can write 

JO 27T Jan 27T <T 



The first integral for r — > is finite; on the other hand the second integral results to be not 
bounded. Infact using the integration variable y = qr we split the y integration region in two 



A.l Asymptotic behaviour of the wave functions. 
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parts, qor < y < yo <C 1 and y > yo; denning the two contribution l\ and I2 respectively, we 
have 



I qor V 

and 



00 



/ 2 ~ / ^Jo(y)(hiy + lnV ~ (InV (A.29) 

Jyo y r r^O r 



So the asymptotic small r behaviour will be divergent, in particular 



4>(r) rZo (ln±.f +1 (A.30) 
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A. 2 Two gluons coupling to quark loop for non-zero transfer 
momentum 

Consider a qq loop for the scattering of a virtual photon J*(q) + ... — ► 7*(<? + I) + q 2 = —Q 2 - 
The momentum transfer I is taken to be purely transversal. We are interested in a multi- 
discontinuity amplitude so that at the lowest order one retains at each discontinuity a simple qq 
state with the quark and antiquark on shell (due to unitarity) and then with their propagators 
substituted by the corresponding delta distributions. The calculations are quite straightforward. 

For example let us consider the most involved case of a transverse polarized photon and the 
gluons (with colours a\ and a 2 ) coupled one (with momentum q\) to the quark (with momentum 
k\) and the other to the anti-quark (momenta q 2 and k 2 respectively) so that q = k± + k 2 - In this 
case the propagators for the quark with momentum k\ + q\ and the anti-quark with momentum 
k 2 have to be replaced by the mass-shell conditions. 

In the high energy limit only the longitudinal component of the gluons polarization give a 
contribution. It is convenient to work in the light cone coordinates and therefore one has for 
the trace present in the amplitude the following expression 

N± = -^-gapTr J±{m f +j-k 1 )j-(m f +j-(k 1 +q 1 ))^ l l(m f --f-(k2+q2)h-(m f --f-k 2 ) (A.31) 

After some algebra and performing the trace one obtains 

N ± = -86 aia2 q 2 _{m 2 f + [a 2 + (1 - a) 2 ]k l± ■ (h + Ql ) ± - a 2 k l± ■ l ± } (A.32) 

where a is defined such that k\ = ag_ . The propagators contain in the denominator two terms 
which can be written, defining Q 2 = —q 2 , as 

D 1 = (k 2 - ml) = -— *— \Q 2 a(l - a) + (m 2 f + k 2 1± ) 
D 2 = p 2 + q 2 ) 2 ~ m}} = --{a(l-a)(Q 2 + l 2 )+m}+ [{h + qi ) ± - al±\ } (A.33) 

Taking also into account the contribution of the 5 distributions, l/4g_o;(l— a), one easily obtains 
the full expression. Similar expressions can be obtained for different photon polarizations and 
configurations of gluons. 

Thus the calculations give for the function /(1,2) corresponding to the loop with gluon 1 
attached to q and gluon 2 attached to q the following expression 

mM) = e)J G daj—^- (A.34) 

where now we use ki to indicate the transverse momentum vectors of the gluons coupled to the 
photon. 

D = (e 2 + k 2 )(e 2 1 + (k + k 1 -alf) (A.35) 
e 2 = Q 2 a(l - a) +m 2 f , e\ = {Q 2 + l 2 )a(l - a) + m) (A.36) 



A. 2 Two gluons coupling to quark loop for non-zero transfer momentum 
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e f and m / are the quark electric charge and mass and the numerator for a transversal photon is 

N 1 = mj + (a 2 + (1 - a) 2 )k(k + k x ) - a 2 kl (A.37) 



and for a longitudinal photon is 



N L = 4Q 2 a 2 (l - a 



(A.38) 



This expression can be conveniently represented as an integral over the colour dipole density 
p created by the qq pair at a given distance in the transverse space: 



ik\r 



f(h,k 2 ) = / d 2 r Pl (r)e tkl 



(A.39) 



From (]A3l) - ( |A38|) one finds for the transverse and longitudinal photons: 



ft(r) 



(2vr) 



/ ^—ialr 



da 



m 2 f K (er)K (e ir ) + (a 2 + (1 - a) 2 )ee 1 iT 1 (er)A 1 (e 1 r) 



ielr 



a(l - a)(l - 2a)— A (eir)A'i(er) 
(2vr)3 



daa 2 (l — a) 2 K()(er)Ko(eir) 



(A.40) 
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A. 3 Coordinate representation of the non-forward G function 

For the non-forward direction a function G(ki, k 2 , ^3) can be defined by a natural generalization 



of the infrared stable function introduced by J.Bartels [20]: 

G(h, k 2 , k 3 ) = -g 2 N c W 2 (h, k 2 , h) 

- D(k u k 2 + h)(u(k 2 ) - uj(k 2 + k 3 )) - D(ki + k 2 , k 3 )(uj{k 2 ) - uj{h + k 2 )) (A.41) 

where D(ki,k 2 ) is the amputated BFKL function with the gluon momenta k\ and k 2 ; to{k) is 
the gluon trajectory and W 2 (k\, k 2 , k 3 ) is the 2^3 vertex integrated with the function D. A 
graphical representation of the function G, following the notation introduced by Bartels, has 



been given in Fig. 3.1. G can be presented as a sum of proper (integral) G\ and improper (terms 



with u>) G 2 parts which, by definition, are given by 

G 2 (k x , k 2 , fcs) = -D(h, k 2 + k 3 )(uj(k 2 ) -uj(k 2 + k 3 ))-D(k 1 + k 2 , h^uik^-uih+k^) (A.42) 
and 

aJk-, fc„ h,\ = n 2 N„ f - 

(2nY 



(2tt) 3 3 < ^ 2 ^ 1 + ( l3-ki-k 2 - k 3 )D(q u q 3 ) 



(k 2 + k 3 f (fci + fc 2 ) 2 k 2 ( kl + k 2 + k 3 ) 2 \ 

+ 7271 Tv?. T~„ u V2tZ ITvI J2J2 (A.4«5j 



V {qi - h) 2 q 3 q 2 (q3 - h) 2 (qi - h) 2 (q3 - h) 2 q 2 q 3 
For the following we recall that the gluon trajectory is 

where the gluon mass m is a regularization parameter (to be set to zero in the infrared finite 
expressions) . 

We shall consider the Fourier transform 

If 3 
G(n,r 2 ,r 3 ) = -j—^ J d 2 k 1 d 2 k 2 d 2 k 3 exp(iY,k j r j )G{k 1 ,k 2 ,k 3 ) (A.45) 

To simplify the notation we denote 



dr(k) = j^-^d 2 kid 2 k 2 d 2 k 3 , dr(q) = - ' ^ 2 (gi + q3 ~ h - k 2 - k 3 ) (A.46) 



and 



dr{p) = d 2 Pl d 2 P3 (A.47) 



A. 3 Coordinate representation of the non-forward G function 
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A. 3.1 Proper part 



Let us start by considering the terms (four) which appear in the proper part ( A.43| ) and which 
will be denoted as Gu, G\2, ••• etc. For the first term one has 

f (k I k l 2 ^ 

Gu = 9 2 N C / dT(k)dT(q)dT(p)D(p 1 ,p 3 ) 2 , L 2 exp()'^ kjrj - %pi - iq 3 p 3 ) (A.48) 
J [qi — K\) q 3 . =l 

After the substitution (£; 2 + k 3 ) 2 — > — (d 2 + d 3 ) 2 , l/q 2 — > — <9~ 3 2 and gi — &i — > £; 2 (all operators 
d refer to coordinate space), choosing q 3 as independent variable one obtains 

Gu = g 2 N c (d 2 + d 3 ) 2 J dT(k)dT(q)dT( P )^ 

exp(ikm + ik 3 r 3 + i(q 3 - k 3 + £; 2 )r 2 - + fc 2 )pi - iq 3 p 3 )d~ 2 D(pi, p 3 ) (A.49) 

Integrations over ki, k 3 and (73 give {2i:) 2 5 2 {ri — pi), (2ir) 2 5 2 (r2 — r 3 ) and (27r) 2 5 2 (r 2 — p 3 ) 
respectively and the remaining integral over fc 2 has to be regularized: 

/ 7,2 1 2 2 ex P^2( r 2 - Pi) = 27r(c - ln|r 2 - p\\), m -> 0, (A.50) 

where 

c = ln(2/m) + V(l) (A.51) 
Doing the integrations by means of <5-functions, we obtain 

Gu = ^{d2 + d 3 ) 2 5 2 {r2 3 ){c-\nr 1 2)d 3 2 D{r 1 ,r 3 ) (A.52) 

Here and in the following r^- = rj — r^. Note that the operator <9 2 + does not act on the 
5-function. Its action on the In ri 2 and D is evidently reduced to <9 2 and d 3 respectively. However 
in the end one has to put r 2 = r 3 so that we obtain finally 

Gu = 9 ^8 2 (v2 3 )d 2 {c-\nr l3 )d 3 2 D{r ll r 3 ) (A.53) 



The Fourier transform of the second term in ( |A.43 ) is obtained from Gu by interchanging 
1^3: 

Gi 2 = ^5 2 (ri 2 )a i 2 ( C -lnri3)ar 2 J D(ri,r 3 ) (A.54) 
Now we consider to the Fourier transform of the third term in ( A.43j ). Its explicit form is 



f k 2 3 
G13 = ~g 2 N c / dT(k)dr(q)dT(p)D(p 1 , p 3 )- — — ^ exp(i V kjVj - iq x p x - iq 3 p 3 ) 

J {qi - hr(q3 - hr ~[ 

(A.55) 

As usual we substitute k 2 — > — d 2 , k\ — > k\ + q\, k 3 — > k 3 + q 3 and choose & 2 as a dependent 
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variable, obtaining 

G13 = 9 2 N c d 2 2 I dr(k)dr(q)dr(p)D( Pl ,p 3 ) 



exp(f(fci + qi)n + i(k 3 + q 3 )r 3 - i(ki + k 3 )r 2 -iq\p\- iq 3 p 3 ) (A.56) 



k 2 k 2 

K 1 K 3 

The integrals over q\ and q 3 give (27r) 2 5 2 (ri — pi) and (27r) 2 <5 2 (r3 — p 3 ) respectively. The two 
integrals over k\ and k 3 give 2n(c — lnri 2 ) and 2ir(c — lnr23), so that one has 

G i3 = 9 ^dl{c-\nr l2 ){c-\nr 23 )D{ n ,r 3 ) (A.57) 
and, after performing the derivatives in r 2 , 

G ^ = -^TT f 2 ¥? + M<= ~ lnr 13 )(5 2 (r 23 ) + S 2 (r l2 ))) D(n,r 3 ) (A.58) 

We are left with the last proper term 

G\a = -g 2 N c ( dT(k)dT(q)dT(p)D(pi , p 3 ) ^— + ^ ^ ^ exp(i ^ fyr.,- - %pi - iq 3 p 3 ) (A.59) 

Substituting all fc's by the derivatives w.r.t. r's and q 7 s by the derivatives w.r.t. p's we can do 
all the integrals to obtaining 

Gu = ^(di + d 2 + d 3 ) 2 S 2 (r 12 )5 2 (r 23 )d^ 2 d 3 2 D(r 1 , r 3 ) (A.60) 

Since the sum X)f=i does not act on the two S- functions, we can transfer the first operator 
to the right and note that, acting on D, d 2 gives zero. So finally one has 

Gu = 4^Vi2)5 2 (r 23 )(3i + d 3 ) 2 d^d 3 2 D( n ,r 3 ) (AM) 

Z7T 



A. 3. 2 Improper part 

A similar notation is used for the improper part: G 2 i, G 22 , ... correspond to the 4 terms in 
( |A.42|) . Take the first term 

3 

G 2 \ = - / dT(k)dT(p)uj(k 2 )D(pi, p 3 ) ex.p(i^2kjrj - ik x p x - i{k 2 + k 3 )p 3 ) (A.62) 

i=l 

We transform this expression in two possible ways. One is evident: we just take uj out of the 
integral as an operator to obtain immediately 

G 21 = -uj(-id 2 )5 2 (r 23 )D{ n ,r 3 ) (A.63) 



A. 3 Coordinate representation of the non-forward G function 
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In this form it is clear that the contribution can be split into a sum of two independent part, 
a holomorphic one acting on r's and an antiholomorphic one acting on r*'s. However this 
expression is difficult to analyze as to its conformal properties, due to the fact that the operator 
u is applied to the <5-function, that is, to both D and a function which is to be integrated with 
G 



A convenient alternative form for G21 is obtained by means of integration over k% and ^3 in 

G21 = -D(n,r 3 ) J ^e ik ^co(k 2 ) = -u(r 23 )D(r u r 3 ) (A.64) 
where u>(r) is just the Fourier transform of uj{k). 



A useful relation, similar to ( |A.50| ) 

d 2 r 



Akr 



If* ^ j ^ 2 



2vr(ci-ln/c) , ci = ln(2/e) + ij>(l) 



allows us to write 



to(k) 



g 2 N c , k 2 g 2 N c r f d 2 r 



In 



ikr 



On noting that 
from eq. ( |A.66 ) we conclude that 



8tt 2 m 2 87r 3 

ci — In m = c — In e 



27r(ci — lnm) 



tjj\r) 



9 2 N C 
8vr 3 



1 



+ 2vr(lne-c)(5 2 (r^ 



(A.65) 

(A.66) 
(A.67) 

(A.68) 



It also follows that the right-hand side of ( A.68 ) does not, in fact, depend on e. Evidently the 
term with In e serves to regularize the singularity of the first term at r = 0. In the following this 
regularization will always be implied, so that 



1 



r 2 + e 2 



+ 2vrlne5 2 (r), e -» 0, 



With this convention we finally obtain 



u>(r) 



9 2 N Cr 1 



8 ^3 V 



(-5 -27TC(5») 



(A.69) 



(A.70) 



The same result could be obtained directly by Fourier transforming the integral representation 
of u(k) but in this way we obtain (|A.69| ) byproduct. Using this we get our final expression 



G 



g 2 N c , 1 



21 



8^ 3 



( 2%c8\r 2 z))D{r 1 ,r 3 ) 



(A.71) 



23 
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The term G23 is obtained from this by the permutation 1 and 3 in the operator acting on D: 



G23 = -^r(4" " ^c5 2 (r 12 ))D(n,r 3 ) 
Improper terms 2 and 4 are calculated in a simpler manner. In fact we have 



(A.72) 



G22 = ~ I dT(k)d,T(p)u}(k 2 + k 3 )D(p 1 , p 3 ) exp(i kjrj - ik\p\ - i(k 2 + k 3 )p 3 ) (A.73) 
J 3=1 

In this case we can substitute to(k2 + k 3 ) by an operator uj(—id P3 ) acting on D. After that all 
integrations are trivially done and we obtain 



and in a similar manner 



G22 = 6 2 (r 23 )uj(-id 3 )D(r 1 ,r 3 ) 
G 2 4 = 5 2 (ri2)uj{-idi)D(ri 1 r 3 ) 



(A.74) 
(A.75) 



A. 3. 3 Infrared stable pieces 

Let us see how to combine different proper and improper terms in order to obtain expressions 
which do not depend on m and so are infrared finite. There are three infrared finite combinations. 
Summing G\\ and G22 we find 

G« = G n + G22 = S 2 (r 23 )^ [2 In 2 + 2^(1) - k 2 3 \nr 2 13 k 3 2 - \nk 2 ]D( ri ,r 3 ) (A.76) 

Here we have used the explicit expression for oj in the momentum space and also used the 
notation k 3 for the operator — id 3 acting in coordinate space. The important point is that the 
terms with mm have cancelled in the sum, so that the expression is infrared finite. Note that 
the expression in the brackets is a half of the Hamiltonian for the two gluons (with a minus 
sign) . It represents the Hamiltonian of the gluon 3 in interaction with the gluon 1 considered as 
an external source. 

The other half of the Hamiltonian appears in the sum Gu and G24, for which we also find 
an infrared finite expression 



G^ = G12 + G 24 = 5 2 (ri 2 )^J [2 In 2 + 2^(1) - k 2 lnr 2 3 k^ 2 - In k\] D(r u r 3 ) 
The rest of the terms give an expression 

G^ = G\ 3 + G14 + G21 + G23 



(A.77) 



9 2 N C 
8tt 3 



-ATT 2 5 2 (r 1 2)5 2 (r 23 )(k 1 + k 3 ) 2 k^ 2 kf + 27^(^3) + 5 2 (r l2 ) lnr 13 - 



rur. 



12' 23 



D(n,r 3 ) 
(A.78) 

It is also infrared finite since all the mm terms cancel. One has only to remember that the 



A. 3 Coordinate representation of the non-forward G function 
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singularities of the last term at r\i 
according to the prescription ( A.69Q . 



and r23 = have to be regularized, if necessary, 
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A. 4 Some useful relations in graphical notation. 



We introduce in the following a graphical notation in order to easily derive some relations 
involving the elementary reggeized gluons (see also in |2jJ for the analytical counterpart). A 
line will denote a reggeized gluon momentum, a line split a dependence in the sum of the 
momenta corresponding to the lines, the squares the interaction vertices defined in (|3.2j) -( |3~4"|) 
and a small circle the contribution of a Regge trajectory (if more than one is present a sum of 
all the possible one-circle graphs is meant). In Fig. [A.l| we show the vertices with their explicit 
definition, related to the 2 —* 2, 2 — ► 3 and 2^4 transitions in (a), (6) and (c) respectively, 
and in Fig. [A. 2 a possible representation of the bootstrap condition. 

It is simple to check the validity of the relation ( |3.33 ) which was needed to calculate the triple 
pomeron vertex Z. The steps are shown in Fig. |A,3 and depend exclusively on the repeated 
application of the relation (b) of Fig. [A.l . The relation in line (a) is directly derived by applying 
on W(4, 1, 3) (first term in the rhs) this relation with the indices (123) — * (413). In (b) and (c) 
in a similar manner other useful relations are derived and, applying them to (a), one finds in 
(d) the identity flOj) . 

Let us show in this graphical representation how to apply the bootstrap condition, giving 
the proof that the solution of eq. ( 3.11 ) is in fact the expression ( |3.13 ). The proof is based on 
the reduction, by means of the bootstrap, of the equation for the 3 gluon system to the equation 
for two gluon in the singlet state, i.e. a BFKL non-homogeneous equation. Precisely one is able 
to show this for an equation with an inhomogeneous term which is a part of the one present 
in ( |3.11 ) so one completes the solution by applying the linearity of the equation. Therefore we 
decompose the second inhomogeneous term L*2^3 in a sum of three terms as in D30 and study 
one of the three similar equations, corresponding to the case with Z?2o(l) and the third and 
fourth terms in the decomposition of Z?2— >3 given in Fig. A. 4. a. 



We can now easily check that the solution of the resulting equation (Fig. A. 4 .6) is given in 



fact by = — Z?2(l ; 2 + 3). With this substitution the equation is represented in Fig. A. 4 
for the Z?2- On applying the bootstrap condition one has the relation in Fig. kA.d and in Fig 



A.4.e two relations are derived from Fig. A.l.fo. Substituting these identities in the relation of 



Fig. A.4.c one finally obtains for D2 the equation in Fig. A. 4./ which is the expected BFKL 
equation. 



A. 4 Some useful relations in graphical notation. 
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1' 2' 



(a) 



1 2 



r 3' 



r 3' 



r 3' 



(b) 



1 2 3 



1 2 3 



1 2 3 



1' 4' 



1' 4' 



(c) 



1 2 3 4 12 34 

Figure A.l: Graphical representation of the elementary vertices. 



Figure A. 2: Graphical representation of the bootstrap condition. 
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(a) 



7V 



41 32 4132 4132 



(b) 



4 1 3 2 4 1 3 2 



T~7\ 

4 13 2 



= ( 



TTA 

4 13 2 



a^ 



)- ( 



4 13 2 41 32 



T^ 
4 13 2 



TTA 

4 13 2 



w 

41 32 



(c) 



TV 

4 1 



3 2 



AAT 

4 13 2 



ATA 

4 13 2 



4 13 2 



4 103 2 



(d) 



4 1 



3 2 



TA 

13 2 



41 32 



A^ 



4 13 2 



4 103 2 



Figure A. 3: Proof of the relation ( ^3^ ). 



A. 4 Some useful relations in graphical notation. 
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(a) 



1 2 3 3 2 1 



"S" "55" "Q" 



(b) 



213 231 123 132 312 321 




+ i ( 



> + ± < 



1 2 3 1 2 3 1 2 3 



12 3 13 2 



1 2 3 2 3 1 3 1 2 



(c) (j-D 




(en 



2 ( 




12 3 13 2 



1 2 3 1 2 3 1 3 2 




12 3 12 3 



(e) 



A 
1 2 3 




1 3 2 




13 2 13 2 




A 



123 123 123 123 



Figure A. 4: Proof of the relation (3.13) 
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